Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,567,347 C→A 6.6% L391L (CTG→CTT mdtM ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,567,3470CA 6.6% 275.7 / 5.3 91L391L (CTG→CTTmdtMmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (39/45);  new base A (4/3);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGA  >  NC_000913/4567210‑4567486
                                                                                                                                         |                                                                                                                                           
aCCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTACTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAg                                                                                                                                           >  1:900210/1‑140 (MQ=255)
     ttGGCGCGTCTGAGTTCGATGATTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATAACAGCTCCACCACTAGCTCGGCTGCCTGATTCTGGCGCACGCGATTGAGCAATCCCCCCAGGGTGa                                                                                                                                      >  1:1392725/1‑140 (MQ=255)
     ttGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGa                                                                                                                                      >  1:1112423/1‑140 (MQ=255)
      tGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGaa                                                                                                                                     >  2:368115/1‑140 (MQ=255)
          gcgTCTGATTTCGATTATTTTATACATAACGTTAAGAATAGAGCCTGACGGCCGGAACGCGCAAAAATCACTGCTCCACCACTAGCTCGGCTGCCTTATGCTGGCGCACGCGATTGACAAATCCCGCCAGGGTGAAAacg                                                                                                                                 <  1:1244082/140‑1 (MQ=255)
          gcgTCTGAGTTCGATGAGTTTATACATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGTTGAAAacg                                                                                                                                 <  1:1244061/140‑1 (MQ=255)
          gcgTATTAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGGGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAacg                                                                                                                                 <  2:710806/140‑1 (MQ=255)
            gTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAacgac                                                                                                                               >  2:1395128/1‑140 (MQ=255)
            gTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAacgac                                                                                                                               >  2:518121/1‑140 (MQ=255)
             tCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAacgacg                                                                                                                              >  1:608294/1‑140 (MQ=255)
                   ttCGATGTGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACCAGCTCGGCTGCCTGTTGCTGGCGCACGCGATTGAGCAAGCCCGCCAGGGTGAAAACGACGATAAcg                                                                                                                        <  1:1472617/140‑1 (MQ=255)
                   ttCGATGAGTTTATTCATACCGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAAcg                                                                                                                        <  1:1473283/140‑1 (MQ=255)
                   ttCGATGAGTTTATTCAGAACGTTAAGATTAGTGCCTGACGGCCGTATCGCGCATAATACACTGCTCCTCCACTCGCTCGGCTGCGTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAACACGACGATAAcg                                                                                                                        <  2:1334556/140‑1 (MQ=255)
                      gATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCAAGGGTGAAAACGACGATAAcgccc                                                                                                                     >  1:5457/1‑140 (MQ=255)
                      gATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCAAGGGTGAAAACGACGATAAcgccc                                                                                                                     >  1:5423/1‑140 (MQ=255)
                      gATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCAAGGGTGAAAACGACGATAAcgccc                                                                                                                     >  1:5401/1‑140 (MQ=255)
                      gATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCAAGGGTGAAAACGACGATAAcgccc                                                                                                                     >  1:5126/1‑140 (MQ=255)
                        tGAGTTTATTCATAACTTTAAGAATAGTGCCTGACGGCCGTATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGTTTGAGCAATCCCGCCAGGGTGAAAACGACGATAAcgcccgc                                                                                                                   <  1:1012809/140‑1 (MQ=255)
                            tttATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCcaca                                                                                                               <  2:263606/140‑1 (MQ=255)
                            tttATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCcaca                                                                                                               <  2:569659/140‑1 (MQ=255)
                            tttATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCcaca                                                                                                               <  2:569630/140‑1 (MQ=255)
                                   ttAACGTTTCGAATCGAGGCATACGGCCGGATGGCGAAAAAATCACTGCACCTCCACTAGGTCGGCTGCCTGATGATGGCGCACGCGATTGAGCAATCCCTCCAGGGTTAATACGACGATATCGCCCGCCACAACGGCTa                                                                                                        <  2:1656880/139‑1 (MQ=255)
                                   aaaaCGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTa                                                                                                        <  2:1656887/138‑1 (MQ=255)
                                   aTAAGGTTAAGAATAGTGCCTGAGGGCCGGATCGCGCTAAAATGAGTGCTCCTCCTCTTGCTCGGCTGCCTGATGCTGGCGCCCGGGATTGAGCAATCCCGCCTGGGTTAAAACGACGATAACGCCCGCCACAACGGCTa                                                                                                        <  2:458255/140‑1 (MQ=255)
                                   aTAACGTTATGATTCGAGCCTGACGGCCGGATCGCGCAAAAAGCACTGCTCCTCCACTTGCTCGGCTGCCTGATGCGGGCGCACGCGAGTGAGCAATCCCGCCAGGGTTAAAACGACGAGAACGCCCGCCACAACGGCTa                                                                                                        <  2:172095/140‑1 (MQ=255)
                                   aTAACGTTAAGAATTGTGCCTGACGGCCGGTTCGCGGATAAATCACTGCTCCTCCAGTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTa                                                                                                        <  2:458336/140‑1 (MQ=255)
                                   aTAACGTTAAGAATAGAGCCTGACGGCCGGATCGCTCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTa                                                                                                        <  2:136387/140‑1 (MQ=255)
                                   aTAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTa                                                                                                        <  2:40090/140‑1 (MQ=255)
                                     aaCGTTAAGATTAGAGCCTGACGGCCGGATCGCGCATAAATCACTGCTCCTCCACTCGCTCGGCTGCCTGTTGCTGGCGCAGGCGATGGAGCAATCCCGCCTGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAAc                                                                                                      <  2:211068/140‑1 (MQ=255)
                                     aaCGTTAAGAATAGTGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAAc                                                                                                      <  2:1185279/140‑1 (MQ=255)
                                     aaCGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTCGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAAc                                                                                                      <  2:212688/140‑1 (MQ=255)
                                     aaCGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAAc                                                                                                      <  2:1185164/140‑1 (MQ=255)
                                          tAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATg                                                                                                 >  2:555963/1‑140 (MQ=255)
                                            aGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGaa                                                                                               >  2:142851/1‑140 (MQ=255)
                                                 tgagCCTGACGGCCGTATCGCGAAAAAATCACTGCTCCTCAACTAGCTCGGCTGCCTGTTGCTGGCGCACGCGTTTGAGCAATCCCGCATGGTTTAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGc                                                                                          <  2:971867/139‑1 (MQ=255)
                                                   agCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCaa                                                                                        >  1:1464393/1‑140 (MQ=255)
                                                   agCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCaa                                                                                        >  1:214951/1‑140 (MQ=255)
                                                   agCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCaa                                                                                        >  1:1464464/1‑140 (MQ=255)
                                                   agCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCaa                                                                                        >  1:721459/1‑140 (MQ=255)
                                                   agCATGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGATACGGCaa                                                                                        >  2:786556/1‑140 (MQ=255)
                                                     ccTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGc                                                                                      >  1:151408/1‑140 (MQ=255)
                                                       tGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGa                                                                                    >  1:768333/1‑140 (MQ=255)
                                                        gcgggCCGGATCGCGCAAAATTCACTGCTCCTCCACTCTCCCGGCGGCCTTATGCTGGCGCCCGCGTTTGAGAAATCCCGCAAGGGTGTAAACGACGATAACGCCCGCCACAACGGATAACTGTTGACACGGCAAGCGAc                                                                                   <  1:1616205/137‑1 (MQ=255)
                                                        gccGGCCGTATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGTTGCTGGCGCACGCGTTTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGAc                                                                                   <  1:1389073/138‑1 (MQ=255)
                                                        gACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGAc                                                                                   >  1:1199326/1‑140 (MQ=255)
                                                        gACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGAc                                                                                   >  1:347128/1‑140 (MQ=255)
                                                        gACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGAc                                                                                   >  1:51217/1‑140 (MQ=255)
                                                        gACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGAc                                                                                   >  1:1198905/1‑140 (MQ=255)
                                                        gACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGAc                                                                                   >  1:1557262/1‑140 (MQ=255)
                                                        gACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGAc                                                                                   >  1:140048/1‑140 (MQ=255)
                                                        gACGGCCGGAACGAGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCTCCACAACGGCTAACAGATGAAACGGCAAGCGAc                                                                                   <  2:932342/140‑1 (MQ=255)
                                                         ccggcggTATCCCGCAAAAATCACTGCTCCTCCACTCGCCCGGCTGCCTGTTGCTGGCGCCCGCGATTTAGAAATCCCGCCTGGGTTAACACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGAcc                                                                                  <  1:981968/139‑1 (MQ=255)
                                                          cGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGAccg                                                                                 <  2:506451/140‑1 (MQ=255)
                                                          cGGCCGGATCGCGCAAAAATCACTGCTCCTCCACGAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGAccg                                                                                 <  2:1651347/140‑1 (MQ=255)
                                                          cGGCCGGATCGCGCAAAAATCACTGCTCCTCCACCAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGAccg                                                                                 <  2:789472/140‑1 (MQ=255)
                                                          cGGCCGGATCGCGCAAAAATCACGGCTCCTCCACGATCTCGGCTGCATTATGGTGGCGCACGCGATTGAGCTTTCCGGCCAGGTTGAAAACGACTATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGAccg                                                                                 <  2:1082992/140‑1 (MQ=255)
                                                           ggCCGGATCGCGAAAAAATCACTGCCCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCTTTTGAGAAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGAccgc                                                                                <  1:1083790/140‑1 (MQ=255)
                                                            gCCGTATCGCGCAAAAATAACTGCTCCCCCACTTGCTCGGCTGCCTGTTGCGGGCGCACGCGATTGACAAATCCCGCCTGGGTTAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGAccgcc                                                                               <  1:123228/140‑1 (MQ=255)
                                                               gtaTCGCGAAAAAATCACTGCACCTCCACTAGCTCGGCGGCCTTATGGTGGCGCCCGCGTTTGAGCATTCCCGCCCGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGtt                                                                            <  1:175493/138‑1 (MQ=255)
                                                               ggATCGCGAAAAAATCACGGCTCCTCCACTCGCTCGGCTGCCTGATGCTGGCGCACGCTATTGAGCAACCCCGCCCGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGtt                                                                            <  1:437867/140‑1 (MQ=255)
                                                                     cgcAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCa                                                                      >  1:697978/1‑140 (MQ=255)
                                                                     cgcAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCa                                                                      >  2:266815/1‑140 (MQ=255)
                                                                      gaaaaaaTCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCTCATGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCAt                                                                     <  2:222920/138‑1 (MQ=255)
                                                                                 tGCTCCTCCACTCGCTGGGCTGCCTGATGCTGGCGCACGCTATTTAGAAATCCCGCCTGGGTCAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGcc                                                          <  1:959337/140‑1 (MQ=255)
                                                                                         cACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCga                                                  >  2:1124625/1‑140 (MQ=255)
                                                                                          cctcgCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCgac                                                 <  1:434224/136‑1 (MQ=255)
                                                                                          aCTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCTGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCgac                                                 <  1:435432/140‑1 (MQ=255)
                                                                                           cTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCCGCGGCCGATTTCgacc                                                >  2:1012684/1‑140 (MQ=255)
                                                                                           cTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCgacc                                                >  2:1074696/1‑140 (MQ=255)
                                                                                               cTCGGCTGCCTTTTGATGGCGAACGCTTTTGAGCAATCCCGCCTGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCgaccgaga                                            <  1:550229/140‑1 (MQ=255)
                                                                                               cTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCgaccgaga                                            >  1:649946/1‑140 (MQ=255)
                                                                                               cTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCgaccgaga                                            >  1:580254/1‑140 (MQ=255)
                                                                                               cTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCgaccgaga                                            >  2:1193804/1‑140 (MQ=255)
                                                                                                   gCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCgaccgagaccga                                        >  2:1172683/1‑140 (MQ=255)
                                                                                                    cTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGAc                                       >  2:90364/1‑140 (MQ=255)
                                                                                                      gCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGTCCGATTACGACCGAGACCGACAt                                     >  2:1101133/1‑140 (MQ=255)
                                                                                                            tGCTGGCGCACGCGTTTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCAt                               <  1:1233506/140‑1 (MQ=255)
                                                                                                            tGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCTCCAt                               >  2:643510/1‑140 (MQ=255)
                                                                                                               tGGCGCACGCGTTTGAGCAATCCCGCCTGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAg                            <  2:137450/140‑1 (MQ=255)
                                                                                                               tGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAg                            <  2:1112423/140‑1 (MQ=255)
                                                                                                               tGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAg                            <  2:563005/140‑1 (MQ=255)
                                                                                                                 gcgcACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGa                          <  2:1496835/140‑1 (MQ=255)
                                                                                                                 gcgcACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGa                          <  2:608819/140‑1 (MQ=255)
                                                                                                                      cgcgTTTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCAcc                     <  1:1199194/140‑1 (MQ=255)
                                                                                                                      cgcgATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCAcc                     >  1:621050/1‑140 (MQ=255)
                                                                                                                              aGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAg             <  2:235307/140‑1 (MQ=255)
                                                                                                                                cAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCg           <  2:685331/140‑1 (MQ=255)
                                                                                                                                   tCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATg        >  2:257052/1‑140 (MQ=255)
                                                                                                                                   tCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATg        >  2:907758/1‑140 (MQ=255)
                                                                                                                                   tCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATg        >  2:1499578/1‑140 (MQ=255)
                                                                                                                                     ccGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCg      >  1:1443460/1‑140 (MQ=255)
                                                                                                                                         cAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGgaga  <  2:41702/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                           
ACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGA  >  NC_000913/4567210‑4567486

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: