Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,706,121:1 +T 10.5% coding (1004/1800 nt) 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,1211.T10.5% 387.6 / 48.2 105coding (1004/1800 nt)lepAback‑translocating elongation factor EF4, GTPase
Reads supporting (aligned to +/- strand):  ref base . (45/49);  new base T (4/7);  total (49/56)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01

Only 100 of 111 total aligned reads displayed.

ACGTAGATAACTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGCCAG  >  NC_000913/2705974‑2706268
                                                                                                                                                    |                                                                                                                                                   
aCGTAGATAACTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGc‑c                                                                                                                                                    >  4:887893/1‑149 (MQ=255)
aCGTAGATAACTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGc                                                                                                                                                      >  1:252945/1‑148 (MQ=255)
 cGTAGATAACTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGc‑c                                                                                                                                                    >  2:108408/1‑148 (MQ=255)
     gATAACTTCTCTTGACGTGGTTTCAACTTCATACACTCCGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAgcg                                                                                                                                               >  3:305781/1‑149 (MQ=255)
     gATAACTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAgc                                                                                                                                                >  1:147090/1‑148 (MQ=255)
      aTAACTTCTCTTGACGTGGTTTCAACTTCATAGACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAGAGC‑CCAgcgc                                                                                                                                              >  4:774167/1‑149 (MQ=255)
        aaCTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAgcgcg                                                                                                                                             <  2:252945/148‑1 (MQ=255)
        aaCTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAgcgcg                                                                                                                                             >  2:135560/1‑148 (MQ=255)
           ttCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAgcgc                                                                                                                                              >  4:801989/1‑144 (MQ=255)
           ttCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAgcgc                                                                                                                                              <  3:801989/144‑1 (MQ=255)
           ttCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTgg                                                                                                                                         >  3:841709/1‑149 (MQ=255)
           ttCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTg                                                                                                                                          <  1:135560/148‑1 (MQ=255)
            tctcTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTgg                                                                                                                                         >  1:51213/1‑148 (MQ=255)
            tctcTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGa                                                                                                                                        <  3:90767/149‑1 (MQ=255)
                ttGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCtt                                                                                                                                    >  4:299921/1‑149 (MQ=255)
                 tGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTgggctt                                                                                                                                     <  2:402398/148‑5 (MQ=255)
                 tGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTgggctt                                                                                                                                     <  2:280112/148‑5 (MQ=255)
                 tGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTGG‑GCttt                                                                                                                                   <  4:306783/149‑1 (MQ=255)
                  gACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCttt                                                                                                                                   <  2:136997/148‑1 (MQ=255)
                     gTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCgg                                                                                                                               >  4:269979/1‑149 (MQ=255)
                     gTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCg                                                                                                                                <  1:610/148‑1 (MQ=255)
                      tGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCgg                                                                                                                               <  2:813616/148‑1 (MQ=255)
                        gTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTc                                                                                                                            >  4:546370/1‑149 (MQ=255)
                        gTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTc                                                                                                                            <  4:762382/149‑1 (MQ=255)
                          ttCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTGG‑GCttt                                                                                                                                   <  3:844329/140‑1 (MQ=255)
                          ttCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTGG‑GCttt                                                                                                                                   >  4:844329/1‑140 (MQ=255)
                          ttCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAgcgc                                                                                                                                              <  4:186611/129‑1 (MQ=255)
                          ttCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAgcgc                                                                                                                                              >  3:186611/1‑129 (MQ=255)
                           tCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTgggctt                                                                                                                                     >  2:621872/1‑134 (MQ=255)
                           tCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTgggctt                                                                                                                                     <  1:621872/138‑5 (MQ=255)
                           tCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATa                                                                                                                         >  4:731923/1‑149 (MQ=255)
                            cAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGa                                                                                                                                        <  3:902653/133‑1 (MQ=255)
                              aCTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGa                                                                                                                       >  2:281940/1‑148 (MQ=255)
                                ttCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTgggctt                                                                                                                                     <  1:733045/133‑5 (MQ=255)
                                ttCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTgggctt                                                                                                                                     >  2:733045/1‑129 (MQ=255)
                                ttCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAAc                                                                                                                     >  1:750900/1‑148 (MQ=255)
                                ttCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACa                                                                                                                    >  3:174333/1‑149 (MQ=255)
                                ttCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACa                                                                                                                    >  4:312625/1‑149 (MQ=255)
                                 tCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACa                                                                                                                    >  2:254945/1‑148 (MQ=255)
                                  cATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGt                                                                                                                  <  3:335761/149‑1 (MQ=255)
                                  cATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGt                                                                                                                  <  4:417392/149‑1 (MQ=255)
                                  cATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGt                                                                                                                  <  3:259509/149‑1 (MQ=255)
                                   aNACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTg                                                                                                                 >  3:148849/1‑149 (MQ=255)
                                    tACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTg                                                                                                                 <  2:750900/148‑1 (MQ=255)
                                     acacTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAg                                                                                                                   >  1:568208/1‑145 (MQ=255)
                                     acacTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAg                                                                                                                   <  2:568208/145‑1 (MQ=255)
                                     acacTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGa                                                                                                                <  1:281940/148‑1 (MQ=255)
                                     acacTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAg                                                                                                               <  3:887893/149‑1 (MQ=255)
                                          aCGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGc‑cc                                                                                                                                                   >  2:82966/1‑108 (MQ=255)
                                          aCGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGc‑cc                                                                                                                                                   <  1:82966/108‑1 (MQ=255)
                                              tCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGtt                                                                                                       <  1:254945/148‑1 (MQ=255)
                                               cGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCa                                                                                                     <  4:841709/149‑1 (MQ=255)
                                                 gcgcAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTAGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAgg                                                                                                   <  4:604672/149‑1 (MQ=255)
                                                  cgcAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAgg                                                                                                   >  2:578441/1‑148 (MQ=255)
                                                   gcAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCAt                                                                                                                          >  1:329726/1‑124 (MQ=255)
                                                   gcAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCAt                                                                                                                          <  2:329726/124‑1 (MQ=255)
                                                     aGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGa                                                                                               <  4:305781/149‑1 (MQ=255)
                                                      gtggtgATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGa                                                                                               <  2:147090/148‑1 (MQ=255)
                                                           gATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTAc                                                                                         <  4:149169/149‑1 (MQ=255)
                                                            aTCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTgggctt                                                                                                                                     >  1:510675/1‑101 (MQ=255)
                                                            aTCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTgggctt                                                                                                                                     <  2:510675/105‑5 (MQ=255)
                                                            aTCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTgggctt                                                                                                                                     <  2:203448/104‑5 (MQ=255)
                                                            aTCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTgggctt                                                                                                                                     >  1:203448/1‑100 (MQ=255)
                                                             tCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTAccc                                                                                       >  3:354018/1‑149 (MQ=255)
                                                               agatccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCa                                                                                      <  1:203774/148‑1 (MQ=255)
                                                                gatccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCa                                                                                                                           >  3:444074/1‑110 (MQ=255)
                                                                gatccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCa                                                                                                                           <  4:444074/110‑1 (MQ=255)
                                                                 atccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCTATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAgcg                                                                                   >  3:805032/1‑149 (MQ=255)
                                                                  tccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAgcgc                                                                                  <  3:269979/149‑1 (MQ=255)
                                                                  tccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAgcgc                                                                                  <  4:191212/149‑1 (MQ=255)
                                                                    cagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGt                                                                                <  4:320813/149‑1 (MQ=255)
                                                                       atcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCAc                                                                             <  3:312625/149‑1 (MQ=255)
                                                                                   ttCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCa                                                                  <  1:286068/148‑1 (MQ=255)
                                                                                    tCCAGACGTTCCTGGATGATCTGCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGATGGCTTCAt                                                                 <  1:783480/148‑1 (MQ=255)
                                                                                    tCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCAt                                                                 >  2:615816/1‑148 (MQ=255)
                                                                                      cAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAgt                                                              >  3:503070/1‑149 (MQ=255)
                                                                                          cGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAgtcgt                                                           <  2:646151/148‑1 (MQ=255)
                                                                                          cGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACC‑AAGC‑TCAGCGCGCTGCAGCTTTCCGGCTCATCGAGCAGTGAGGCATCGTTCATGCTGAGTTTACCCAGCGGGTCACGGAAGGCCTCACAgacgtc                                                          <  2:869133/148‑1 (MQ=255)
                                                                                                tGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACt                                                    >  3:716588/1‑149 (MQ=255)
                                                                                                  gatgatCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACtt                                                   >  2:424844/1‑148 (MQ=255)
                                                                                                      atCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCgg                                              >  3:314850/1‑149 (MQ=255)
                                                                                                        ctcCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAAGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCggg                                             >  2:397060/1‑148 (MQ=255)
                                                                                                                       gCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGtat                              <  1:108408/148‑1 (MQ=255)
                                                                                                                           aGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCgg                                   <  2:629080/139‑1 (MQ=255)
                                                                                                                           aGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCgg                                   >  1:629080/1‑139 (MQ=255)
                                                                                                                             gAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGCCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCg                       <  3:774167/149‑1 (MQ=255)
                                                                                                                                gCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGt                     <  1:450805/148‑1 (MQ=255)
                                                                                                                                gCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGt                     <  1:578441/148‑1 (MQ=255)
                                                                                                                                 ccGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGttt                   <  4:815440/149‑1 (MQ=255)
                                                                                                                                 ccGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGtt                    <  2:51213/148‑1 (MQ=255)
                                                                                                                                    cAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGa                 >  2:906736/1‑148 (MQ=255)
                                                                                                                                         gAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCt           >  4:27051/1‑149 (MQ=255)
                                                                                                                                             ccAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTaa        <  1:906736/148‑1 (MQ=255)
                                                                                                                                             ccAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTaa        >  2:611235/1‑148 (MQ=255)
                                                                                                                                              cAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAg      >  3:395716/1‑149 (MQ=255)
                                                                                                                                              cAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAg      >  4:154461/1‑149 (MQ=255)
                                                                                                                                                aaGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGcc    <  3:627072/149‑1 (MQ=255)
                                                                                                                                                aaGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGcc    <  4:805032/149‑1 (MQ=255)
                                                                                                                                                aaGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGcc    >  3:894258/1‑149 (MQ=255)
                                                                                                                                                   c‑ccAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGCCAg  <  1:615816/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACGTAGATAACTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGCCAG  >  NC_000913/2705974‑2706268

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