Predicted mutation
evidence seq id position mutation annotation gene description
RA CP028101 880,795 A→T K1889N (AAA→AAT)  C7Y58_04335 → autotransporter domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP028101880,7950AT100.0% 147.9 / NA 45K1889N (AAA→AAT) C7Y58_04335autotransporter domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (24/21);  total (24/21)

AGATAAAGAATTTAAATATCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAAGATATAGGAAAATCAAAAGAACATATGTTACAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAATATCAGGA  >  CP028101/880555‑880885
                                                                                                                                                                                                                                                |                                                                                          
aGATAAAGAATTTAAATATCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGaa                                                                                   >  2:422121/1‑250 (MQ=255)
 gATAAAGAATTTAAATATCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACa                                                                                 >  2:707373/1‑251 (MQ=255)
               atatCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTgcaagcaa                                                                   >  1:334700/1‑251 (MQ=255)
               atatCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTgcaagcaa                                                                   >  1:365157/1‑251 (MQ=255)
                 atCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGGGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACCTGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCTAGCAAAg                                                                 >  1:655540/1‑251 (MQ=11)
                  tCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAgg                                                                >  2:315752/1‑251 (MQ=255)
                  tCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAgg                                                                >  1:69262/1‑251 (MQ=255)
                  tCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAgg                                                                >  2:541400/1‑251 (MQ=255)
                         gaTGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGatt                                                         >  2:161111/1‑251 (MQ=255)
                         gaTGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGatt                                                         >  1:330903/1‑251 (MQ=255)
                          aTGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGatt                                                         <  2:315894/250‑1 (MQ=255)
                            gAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGattatt                                                      >  1:118083/1‑251 (MQ=255)
                             aaTGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTTATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGAttatta                                                     <  2:326725/251‑1 (MQ=255)
                             aaTGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGAttatta                                                     <  2:270424/251‑1 (MQ=255)
                              aTGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATGCAAAACAGATACAGCTGGAGTGATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTaa                                                    <  1:536106/251‑1 (MQ=14)
                              aTGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCACGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTaa                                                    >  1:510779/1‑251 (MQ=255)
                              aTGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTaa                                                    >  2:654334/1‑251 (MQ=255)
                                ggAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGCGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGATCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTaaaa                                                  >  2:637048/1‑251 (MQ=255)
                                    aaCAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTATAGGAGAATACAAAACAGATACAGCAGGAGTTATACACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATc                                                <  1:571004/249‑1 (MQ=255)
                                    aaCAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATc                                                <  1:161111/249‑1 (MQ=255)
                                    aaCAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATc                                                <  1:110031/249‑1 (MQ=255)
                                    aaCAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATc                                                <  2:69262/249‑1 (MQ=255)
                                         cTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGatt                                                         <  1:37894/235‑1 (MQ=255)
                                         cTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGatt                                                         >  2:37894/1‑235 (MQ=255)
                                         cTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCa                                         >  2:565891/1‑251 (MQ=255)
                                           tCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGatt                                                         >  1:364634/1‑233 (MQ=255)
                                           tCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGatt                                                         <  2:364634/233‑1 (MQ=255)
                                                  aTTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACCGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGACATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTgatgat                                <  1:637048/251‑1 (MQ=255)
                                                     cAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAAc                             <  1:707373/251‑1 (MQ=255)
                                                          aaaataaaaGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAg                        <  2:94689/251‑1 (MQ=255)
                                                           taataaaaGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAg                        >  2:379378/2‑250 (MQ=255)
                                                            aataaaaGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTACGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCt                      <  1:654334/251‑1 (MQ=255)
                                                             ataaaaGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATCGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCtt                     <  1:565891/251‑1 (MQ=255)
                                                                    cATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTgg              >  1:270202/1‑251 (MQ=16)
                                                                     aTTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGCATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCAATTGATGATAACAATAGCTTGAATTGGa             >  2:315464/1‑251 (MQ=255)
                                                                          gAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACaaaa        >  2:33491/1‑249 (MQ=255)
                                                                           aaCTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAatat       >  1:449642/1‑251 (MQ=12)
                                                                           aaCTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAatat       >  1:361939/1‑251 (MQ=12)
                                                                            aCTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAATATc      <  2:261946/251‑1 (MQ=12)
                                                                               agagGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGAc            <  1:401561/242‑1 (MQ=9)
                                                                               agagGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGATGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGAc            >  2:401561/1‑242 (MQ=9)
                                                                                 aggagAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAATATCAgg   <  1:422121/249‑1 (MQ=2)
                                                                                 aggagAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAATATCAgg   <  2:90810/249‑1 (MQ=2)
                                                                                 aggagAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAATATCAGga  <  2:237929/250‑1 (MQ=2)
                                                                                 aggagAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAACAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAATATCAGga  <  1:121896/250‑1 (MQ=2)
                                                                                                                                                                                                                                                |                                                                                          
AGATAAAGAATTTAAATATCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAAGATATAGGAAAATCAAAAGAACATATGTTACAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAATATCAGGA  >  CP028101/880555‑880885

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑