Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,706,607 A→C 10.2% L173R (CTC→CGC) ‡ lepA ← back‑translocating elongation factor EF4, GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,706,6070AC10.2% 281.7 / 23.0 88L173R (CTC→CGC) ‡lepAback‑translocating elongation factor EF4, GTPase
Reads supporting (aligned to +/- strand):  ref base A (46/33);  new base C (4/5);  total (50/38)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

TTTGTCGCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGCACAACTTCGA  >  NC_000913/2706459‑2706754
                                                                                                                                                    |                                                                                                                                                   
tttGTCGCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCga                                                                                                                                                     <  3:666070/149‑1 (MQ=255)
tttGTCGCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCga                                                                                                                                                     <  3:486402/149‑1 (MQ=255)
      gCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCTCGAACGt                                                                                                                                               >  4:336829/1‑149 (MQ=255)
       cccTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGt                                                                                                                                               >  1:825086/1‑148 (MQ=255)
        ccTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTcc                                                                                                                                             <  3:700256/149‑1 (MQ=255)
         cTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCTCGAACGTctg                                                                                                                                            <  4:625281/149‑3 (MQ=255)
         cTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCTCGAACGTctg                                                                                                                                            <  3:590175/149‑3 (MQ=255)
          ttACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGTTAGTTGTCGAACTATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCTGAATGTCGCGCACCAGACGTTCTCGAACGtatg                                                                                                                                            <  2:599934/148‑4 (MQ=255)
          ttACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCTCGAACGTctgg                                                                                                                                           <  3:405017/149‑4 (MQ=255)
          ttACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCTCGAACGTctgg                                                                                                                                           <  3:336829/149‑4 (MQ=255)
           tACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTg                                                                                                                                           >  1:590205/1‑148 (MQ=255)
                       cGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGcc                                                                                                                               >  2:602762/1‑148 (MQ=255)
                       cGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGcc                                                                                                                               >  1:781249/1‑148 (MQ=255)
                        gTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCg                                                                                                                              >  2:797569/1‑148 (MQ=255)
                         tttttAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCgg                                                                                                                             <  2:825086/148‑1 (MQ=255)
                          ttttAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGt                                                                                                                            >  2:729207/1‑148 (MQ=255)
                            ttAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGtttt                                                                                                                         <  4:802978/149‑1 (MQ=255)
                             tAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGtttt                                                                                                                         >  1:171737/1‑148 (MQ=255)
                                 aCGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCt                                                                                                                     >  1:719081/1‑148 (MQ=255)
                                   ggATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGa                                                                                                                   <  1:192624/148‑1 (MQ=255)
                                     aTAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACa                                                                                                                >  4:415534/1‑149 (MQ=255)
                                         gTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAgcgc                                                                                                            >  3:904975/1‑149 (MQ=255)
                                          tGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCa                                                                                                           >  4:706610/1‑149 (MQ=255)
                                            aaaCAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCAc                                                                                                          >  1:531632/1‑148 (MQ=255)
                                                aacGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGt                                                                                                     >  4:163440/1‑149 (MQ=255)
                                                   gCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCgg                                                                                                  >  3:609231/1‑149 (MQ=255)
                                                      cAGGTAGTTGTCGAACCATGAGTCGATAATTAGTTCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGtt                                                                                                                           >  1:226635/1‑121 (MQ=255)
                                                      cAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGtt                                                                                                                           <  2:226635/121‑1 (MQ=255)
                                                       aGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGc                                                                                              <  3:427485/149‑1 (MQ=255)
                                                        ggTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACCTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGc                                                                                              <  2:781249/148‑1 (MQ=255)
                                                            gTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTc                                                                                                                                              <  4:371301/96‑1 (MQ=255)
                                                            gTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTc                                                                                                                                              >  3:371301/1‑96 (MQ=255)
                                                             ttGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGa                                                                                         <  2:171737/148‑1 (MQ=255)
                                                              tGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGAt                                                                                        >  2:81467/1‑148 (MQ=255)
                                                                 cGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGtt                                                                                                                           <  1:492301/110‑1 (MQ=255)
                                                                 cGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGtt                                                                                                                           >  2:492301/1‑110 (MQ=255)
                                                                   aaCCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCcgac                                                                                  <  3:163440/149‑1 (MQ=255)
                                                                     ccATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCcgacg                                                                                 >  2:589702/1‑148 (MQ=255)
                                                                     ccATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCcgacg                                                                                 <  2:380382/148‑1 (MQ=255)
                                                                         gAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGatat                                                                             <  1:164935/148‑1 (MQ=255)
                                                                         gAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATc                                                                            >  3:351182/1‑149 (MQ=255)
                                                                            tCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTc                                                                         >  3:74120/1‑149 (MQ=255)
                                                                             cGAGAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTc                                                                         >  2:821822/1‑148 (MQ=255)
                                                                              gATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCaa                                                                       >  3:893387/1‑149 (MQ=255)
                                                                                tAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAAtt                                                                     >  4:842145/1‑149 (MQ=255)
                                                                                     aGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATttcttc                                                                <  3:892618/149‑1 (MQ=255)
                                                                                      gTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATttcttc                                                                >  1:750483/1‑148 (MQ=255)
                                                                                         ccTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGc                                                             <  1:602762/148‑1 (MQ=255)
                                                                                           tGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCAc                                                          <  3:589736/149‑1 (MQ=255)
                                                                                            gCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGATGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCAcg                                                         >  4:738785/1‑149 (MQ=255)
                                                                                                cGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGtt                                                     <  4:74120/149‑1 (MQ=255)
                                                                                                 gggCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTc                                                    >  3:605187/1‑149 (MQ=255)
                                                                                                  ggCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTc                                                    <  1:81467/148‑1 (MQ=255)
                                                                                                   gCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTACGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCa                                                   <  1:267836/148‑1 (MQ=255)
                                                                                                       tCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCATCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGAt                                               <  1:359502/148‑1 (MQ=255)
                                                                                                       tCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGAt                                               >  2:161300/1‑148 (MQ=255)
                                                                                                       tCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGAt                                               >  1:343854/1‑148 (MQ=255)
                                                                                                       tCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGAt                                               >  1:333936/1‑148 (MQ=255)
                                                                                                          ggATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGc                                           >  4:487045/1‑149 (MQ=255)
                                                                                                           gATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCt                                          >  4:770177/1‑149 (MQ=255)
                                                                                                            aTCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTg                                         >  4:716220/1‑149 (MQ=255)
                                                                                                               gCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCg                                      >  3:722296/1‑149 (MQ=255)
                                                                                                                ccTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCgg                                     >  4:482688/1‑149 (MQ=255)
                                                                                                                ccTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCg                                      >  2:582984/1‑148 (MQ=255)
                                                                                                                  ttCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGc                                    <  1:301687/148‑1 (MQ=255)
                                                                                                                       gcggcggAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCa                              >  4:129932/1‑149 (MQ=255)
                                                                                                                          gcggAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATc                           <  4:609231/149‑1 (MQ=255)
                                                                                                                             gAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCtt                         <  1:797569/148‑1 (MQ=255)
                                                                                                                               aTGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGtt                      <  4:904975/149‑1 (MQ=255)
                                                                                                                               aTGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGtt                      <  3:706610/149‑1 (MQ=255)
                                                                                                                               aTGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGtt                      <  3:487045/149‑1 (MQ=255)
                                                                                                                               aTGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGtt                      <  3:312405/149‑1 (MQ=255)
                                                                                                                                 gTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCa                    <  3:770177/149‑1 (MQ=255)
                                                                                                                                 gTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCa                    <  3:482688/149‑1 (MQ=255)
                                                                                                                                 gTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCa                    <  3:415534/149‑1 (MQ=255)
                                                                                                                                   cgcgCACCAGACGTTCGAGAACGTCTTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTAGTTCAGt                  >  3:778751/1‑149 (MQ=255)
                                                                                                                                    gcgcACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTa                 >  4:357560/1‑149 (MQ=255)
                                                                                                                                     cgcACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCGATCTTGtt                      <  3:833609/143‑1 (MQ=255)
                                                                                                                                     cgcACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGtt                      >  4:833609/1‑143 (MQ=255)
                                                                                                                                      gcACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGt                       >  2:908431/1‑141 (MQ=255)
                                                                                                                                      gcACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGt                       <  1:908431/141‑1 (MQ=255)
                                                                                                                                          aggACGTTCTCGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGca           >  4:77292/3‑149 (MQ=255)
                                                                                                                                          aggACGTTCTCGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGc            >  1:146170/3‑148 (MQ=255)
                                                                                                                                           aGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGca           >  2:589277/1‑148 (MQ=255)
                                                                                                                                            gACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGcaca         >  3:249995/1‑149 (MQ=255)
                                                                                                                                            gACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGcac          >  2:668771/1‑148 (MQ=255)
                                                                                                                                             aCGTTCTCGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGTGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGcaca         >  2:167984/1‑148 (MQ=255)
                                                                                                                                             aCGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATttc                                                                   <  2:145460/90‑1 (MQ=255)
                                                                                                                                             aCGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATttc                                                                   >  1:145460/1‑90 (MQ=255)
                                                                                                                                                ttCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGCACAACtt     >  3:89942/1‑149 (MQ=255)
                                                                                                                                                   gagaACGTCCTGCACACCAACGCCGGTTTTCGCTGGACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGCACAAATTcta  <  4:89942/149‑3 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTTGTCGCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGCACAACTTCGA  >  NC_000913/2706459‑2706754

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: