Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,706,134 Δ1 bp 5.2% coding (991/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,1340A. 5.2% 329.1 / 10.5 97coding (991/1800 nt)lepAback‑translocating elongation factor EF4, GTPase
Reads supporting (aligned to +/- strand):  ref base A (42/50);  new base . (2/3);  total (44/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

TCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGCCAGGCAGCGCCTTTTCT  >  NC_000913/2705986‑2706282
                                                                                                                                                     |                                                                                                                                                    
tctcTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGa                                                                                                                                                      <  1:74404/149‑1 (MQ=255)
  tctTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGc                                                                                                                                                    >  4:491913/1‑149 (MQ=255)
     tGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTGG‑GCttt                                                                                                                                                 <  3:426674/149‑1 (MQ=255)
          tGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGc                                                                                                                                            <  2:591869/149‑1 (MQ=255)
          tGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGc                                                                                                                                            >  2:97285/1‑149 (MQ=255)
            gTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTc                                                                                                                                          <  3:714604/149‑1 (MQ=255)
             tttCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCa                                                                                                                                         >  4:859188/1‑149 (MQ=255)
              ttCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCAt                                                                                                                                        <  2:836132/149‑1 (MQ=255)
               tCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATa                                                                                                                                       <  2:92637/149‑1 (MQ=255)
                 aaCTTCATCCACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATa                                                                                                                                       <  3:203058/147‑1 (MQ=255)
                 aaCTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATa                                                                                                                                       >  4:203058/1‑147 (MQ=255)
                  aCTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGaa                                                                                                                                    <  4:462030/149‑1 (MQ=255)
                    ttCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACa                                                                                                                                  <  1:229352/149‑1 (MQ=255)
                     tCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAg                                                                                                                                 <  1:355017/149‑1 (MQ=255)
                      cATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGt                                                                                                                                <  4:732791/149‑1 (MQ=255)
                       aTACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTg                                                                                                                               <  2:666069/149‑1 (MQ=255)
                         acacTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTTCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGTCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCAAAGAACAGTGAg                                                                                                                             <  4:80644/149‑1 (MQ=255)
                         acacTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAg                                                                                                                             <  3:158382/149‑1 (MQ=255)
                         acacTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAg                                                                                                                             <  2:307936/149‑1 (MQ=255)
                         acacTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAg                                                                                                                             <  2:260883/149‑1 (MQ=255)
                             tACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCAt                                                                                                                         >  2:430954/1‑149 (MQ=255)
                             tACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGG‑GCtttg                                                                                                                                                >  4:423592/1‑124 (MQ=255)
                             tACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGG‑GCtttg                                                                                                                                                <  3:423592/125‑2 (MQ=255)
                                   cGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCa                                                                                                                   <  3:555721/149‑1 (MQ=255)
                                         aGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGa                                                                                                             <  3:391589/149‑1 (MQ=255)
                                         aGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGa                                                                                                             <  3:376807/149‑1 (MQ=255)
                                               gATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGt                                                                                                           >  4:518438/1‑145 (MQ=255)
                                               gATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGt                                                                                                           <  3:518438/145‑1 (MQ=255)
                                               gATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTAc                                                                                                       >  4:170734/1‑149 (MQ=255)
                                                    gatccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAgc                                                                                                  <  3:165594/149‑1 (MQ=255)
                                                     atccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAgcg                                                                                                 >  4:724313/1‑149 (MQ=255)
                                                      tccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGG‑GCttt                                                                                                                                                 <  4:358166/99‑1 (MQ=255)
                                                      tccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGG‑GCttt                                                                                                                                                 >  3:358166/1‑99 (MQ=255)
                                                      tccaaatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAgcgc                                                                                                >  2:27981/1‑149 (MQ=255)
                                                          gatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCa                                                                                            >  2:719345/1‑149 (MQ=255)
                                                          gatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCa                                                                                            >  1:613285/1‑149 (MQ=255)
                                                            tcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCtt                                                                                                                                                  >  4:163078/1‑93 (MQ=255)
                                                            tcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCtt                                                                                                                                                  <  3:163078/93‑1 (MQ=255)
                                                               taatCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGACGAGGAACCCGCAGCGGAAACCAAAGT‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAAAAGTTAGGCATCGTTCAGGCTGAGTTTAGTCAGCGCGTCACGGaa                                                                                       <  4:828383/146‑1 (MQ=255)
                                                                 ttCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAg                                                                                                                             <  1:662752/109‑1 (MQ=255)
                                                                 ttCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAg                                                                                                                             >  2:662752/1‑109 (MQ=255)
                                                                 ttCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAgg                                                                                     >  1:390490/1‑149 (MQ=255)
                                                                  tCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGaa                                                                                       <  1:121499/146‑1 (MQ=255)
                                                                  tCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGaa                                                                                       >  2:121499/1‑146 (MQ=255)
                                                                    aCGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCtt                                                                                  >  3:10492/1‑149 (MQ=255)
                                                                           aGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAgtc                                                                           <  1:430954/149‑1 (MQ=255)
                                                                           aGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAgac                                                                           >  4:554088/1‑147 (MQ=255)
                                                                              cGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAgtcgtc                                                                        <  1:126946/149‑1 (MQ=255)
                                                                                  ccTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGaa                                                                    <  2:384275/149‑1 (MQ=255)
                                                                                   cTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAAc                                                                   >  4:82646/1‑149 (MQ=255)
                                                                                   cTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAAc                                                                   <  4:351095/149‑1 (MQ=255)
                                                                                    tGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACt                                                                  >  4:841543/1‑149 (MQ=255)
                                                                                    tGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACt                                                                  >  4:634691/1‑149 (MQ=255)
                                                                                     ggATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACtt                                                                 >  2:322388/1‑149 (MQ=255)
                                                                                           tctcCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCggg                                                           >  1:7095/1‑149 (MQ=255)
                                                                                            ctcCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGa                                                          >  1:185221/1‑149 (MQ=255)
                                                                                             tcCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGaa                                                         <  2:822982/149‑1 (MQ=255)
                                                                                             tcCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGaa                                                         >  4:844438/1‑149 (MQ=255)
                                                                                             tcCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGaa                                                         <  1:27981/149‑1 (MQ=255)
                                                                                             tcCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGATGGATCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTCCCCAGCGCGTCACGGAACGCTTCATAGTCGTCGGAACTTACCGGGat                                                         <  2:158347/149‑2 (MQ=255)
                                                                                               cATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACa                                                       >  2:123839/1‑149 (MQ=255)
                                                                                               cATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACa                                                       <  1:443335/149‑1 (MQ=255)
                                                                                               cATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACa                                                       >  1:484567/1‑149 (MQ=255)
                                                                                                aTGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAg                                                      <  1:123839/149‑1 (MQ=255)
                                                                                                 tgtgCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGa                                                     >  2:359073/1‑149 (MQ=255)
                                                                                                 tgtgCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGa                                                     <  2:2807/149‑1 (MQ=255)
                                                                                                   tgCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGAcc                                                   <  3:389785/149‑1 (MQ=255)
                                                                                                       gcagGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCg                                               >  4:813376/1‑149 (MQ=255)
                                                                                                         agGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGa                                                     >  2:351484/1‑141 (MQ=255)
                                                                                                         agGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGa                                                     <  1:351484/141‑1 (MQ=255)
                                                                                                             cGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATAcc                                         >  3:80354/1‑149 (MQ=255)
                                                                                                                ggAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGc                                      <  1:97285/149‑1 (MQ=255)
                                                                                                                 gAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCg                                     >  3:404680/1‑149 (MQ=255)
                                                                                                                  aaGCCGGAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACGGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGAATCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCgg                                    <  4:635207/149‑1 (MQ=255)
                                                                                                                  aaGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCgg                                    <  4:541503/149‑1 (MQ=255)
                                                                                                                     ccGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGttt                                 <  2:613285/149‑1 (MQ=255)
                                                                                                                           cGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACTGACCGGCGTATACCTGCGGTTTGACttt                           <  3:554088/149‑1 (MQ=255)
                                                                                                                            ggAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTc                          >  4:792008/1‑149 (MQ=255)
                                                                                                                            ggAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTc                          <  4:859081/149‑1 (MQ=255)
                                                                                                                            ggAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTc                          >  1:101433/1‑149 (MQ=255)
                                                                                                                                aCCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTaa                      >  1:795730/1‑149 (MQ=255)
                                                                                                                                 ccAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTaaa                     >  3:812937/1‑149 (MQ=255)
                                                                                                                                 ccAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTaaa                     <  3:535647/149‑1 (MQ=255)
                                                                                                                                 ccAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTaaa                     <  2:505618/149‑1 (MQ=255)
                                                                                                                                 ccAAAGC‑CCAGCGCCGTGGAGCTTTCCGGCGCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTGACCTAGCGCGTCACGGAAGGCGTCATAGTCGGCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACATTCTTaaa                     <  2:161571/149‑1 (MQ=255)
                                                                                                                                  cAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGa                               >  4:859170/1‑138 (MQ=255)
                                                                                                                                  cAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGa                               >  3:508609/1‑138 (MQ=255)
                                                                                                                                  cAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGa                               <  4:508609/138‑1 (MQ=255)
                                                                                                                                  cAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGa                               <  3:859170/138‑1 (MQ=255)
                                                                                                                                  cAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTc                          <  1:68738/143‑1 (MQ=255)
                                                                                                                                  cAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAg                    >  4:410506/1‑149 (MQ=255)
                                                                                                                                  cAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACAGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTc                          >  2:68738/1‑143 (MQ=255)
                                                                                                                                        tccAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGCCAgg               >  2:32472/2‑149 (MQ=255)
                                                                                                                                         ccAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGCCAGGc              >  4:255849/1‑149 (MQ=255)
                                                                                                                                          cAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGCCAGGCa             >  3:611010/1‑149 (MQ=255)
                                                                                                                                            gcgcgcTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGCCAGGCAgc           <  1:71279/149‑1 (MQ=255)
                                                                                                                                              gcgcTGGAGCTTTCCGGCTCATAGTACAGTGAGGCATCGTTCAGGATGAGTTTACCCAGCGCGTCACGGAAGGCTTTATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGCCAGGCAgcgc         <  1:371723/149‑1 (MQ=255)
                                                                                                                                                    gAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGCCAGGCAGCGCCTTTTc   <  3:813376/149‑1 (MQ=255)
                                                                                                                                                     aGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGCCAGGCAGCGCCTTTTCt  <  1:502533/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGC‑CCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGCCAGGCAGCGCCTTTTCT  >  NC_000913/2705986‑2706282

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