Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,337,144 C→T G38E (GGG→GAG)  mlaB ← intermembrane phospholipid transport system protein MlaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,337,1440CT97.4% 121.3 / NA 40G38E (GGG→GAG) mlaBintermembrane phospholipid transport system protein MlaB
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (20/18);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TGCCGGACAGGCGGTTAACGCCATATCCGGCCTGAAAAAATTTAACGAGGCAGAACATCAGCAGGCAAATTATACAATTTCGCCAGGGTATACACTTTGTCGTTTACCCCCTGAAGCGTCACATTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCCCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTCAGTTGCGCAGTCAGGCCGTCGATACCTTTGGTACGCAGCAGCGTTCCCCACTCGTTTTGTTTGGTGGTGATCATACTGACGCCTTCAG  >  NC_000913/3336922‑3337389
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     
tGCCGGACAGGCGGTTAAAGCCCTATCCCGCCTGAAAAAATTTAACGAGGCAGAAAAACAGCAGGCCAATTATACAATTTTGCCAGGGGATACACCTTGTCGTTTACCCCCTGAAGCGTCACATTGTTTCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCCGCAGGGCCAGTCCCCCCCTATCCCCGCGGGGGACACGGCTAAGATCGAAGCAGGTAATACCCTTCACCCCCTCCTCACGCATTTcccc                                                                                                                                                                                                                           >  2:393806/1‑250 (MQ=255)
            ggTTAACGCCATATCCGGCCTGAAAAAATTTAACGAGGCAGAACATCAGCAGGCAAATTATACAATTTCGCCAGGGTATACACTTTGTCGTTTACCCCCTGAAGCGTCACATTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCaa                                                                                                                                                                                                               <  1:280320/251‑1 (MQ=255)
             gTTAACGCCATATCCGGCCTGAAAAAATTTAACGAGGCAGAACATCAGCAGGCAAATTATACAATTTCGCCAGGGTATACACTTTGTCGTTTACCCCCTGAAGCGTCACATTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGTGAGACACGGCTAAGATCGCTGCAGGCAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGAAGCaa                                                                                                                                                                                                               >  2:19614/1‑250 (MQ=255)
                                       aTTTAACGAGGCAGAACATCAGCAGGCAAATTATACAATTTCGCCAGGGTATACACTTTGTCGTTTACCCCCTGAAGCGTCACATTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATaa                                                                                                                                                                                    <  1:187826/251‑1 (MQ=255)
                                           aaCGAGGCAGAACATCAGCAGGCAAATTATACAATTTCGCCAGGGTATACACTTTGTCGTTTACCCCCTGAAGCGTCACATTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGcc                                                                                                                                                                                >  2:273898/1‑251 (MQ=255)
                                             cGAGGCAGAACATCAGCAGGCAAATTATACAATTTCGCCAGGGTATACACTTTGTCGTTTACCCCCTGAAGCGTCACATTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAg                                                                                                                                                                              <  1:159902/251‑1 (MQ=255)
                                                 gCAGAACATCAGCAGGCAAATTATACAATTTCGCCAGGGTATACACTTTGTCGTTTACCCCCTGAAGCGTCACATTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCgtg                                                                                                                                                                          >  1:379397/1‑251 (MQ=255)
                                                              agGCAAATTATACAATTTCGCCAAGGTATACACTTTGTCGTTTACCCCATGAAGCGTCACCTTGTTGCACTGCTTTTTAGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGAATAACGCCAGCGTGTCACCCGGCTGc                                                                                                                                                              >  2:251455/1‑250 (MQ=255)
                                                                     atataCAATTTCGCCAGGGTATACACTTTGTCGTTTACCCCCTGAAGCGTCACATTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCTCCGCTTCCTCACGCATTTCCCAAAGCGGGAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCt                                                                                                                                                      <  2:253298/250‑1 (MQ=255)
                                                                       ataNAATTTCGCCAGGGTATACACTTTGTCGTTTACCCCCTGAAGCGTCACATTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCa                                                                                                                                                    <  2:379397/251‑1 (MQ=255)
                                                                            aTTTCGCCAGGGTATACACTTTGTCGTTTACCCCCTGAAGCGTCACATTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGAc                                                                                                                                               >  1:1524/1‑251 (MQ=255)
                                                                               tCGCCAGGGTATACACTTTGTCGTTTACCCCCTGAAGCGTCACATTGTTGCCCTGCTTTGTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTc                                                                                                                                             <  1:367000/250‑1 (MQ=255)
                                                                                 gCCAGGGTATACACTTTGTCGTTTACCCACTGAAGCGTCACATTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCt                                                                                                                                          <  2:120806/251‑1 (MQ=255)
                                                                                        tataCACTTTGTCGTTTACCCCCTGAAGCGTCACATTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCattatt                                                                                                                                   >  1:59849/1‑251 (MQ=255)
                                                                                                ttGTCGTTTACCCCCTGACGCGTCACATTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTctct                                                                                                                           >  2:35338/1‑251 (MQ=255)
                                                                                                ttGTCGTTTACCCCCTGAAGCGTCACAGTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAAAGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTtct                                                                                                                             <  1:173019/249‑1 (MQ=255)
                                                                                                           ccccTGAAGCGTCACCTTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGCTGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGGGTGAt                                                                                                                >  2:26790/1‑251 (MQ=255)
                                                                                                              cTGAAGCGTCACATTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGAtttt                                                                                                             >  1:150153/1‑251 (MQ=255)
                                                                                                                        ccaTTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGaa                                                                                                   <  1:16987/250‑1 (MQ=255)
                                                                                                                          aTTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCatt                                                                                                                                      >  1:141569/1‑214 (MQ=255)
                                                                                                                          aTTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCatt                                                                                                                                      <  2:141569/214‑1 (MQ=255)
                                                                                                                             gttgCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCg                                                                                               <  1:393806/250‑1 (MQ=255)
                                                                                                                                 cccTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTc                                                                                          >  1:367557/1‑251 (MQ=255)
                                                                                                                                         tttCGCAAGATCGATCAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTACGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTCAGTTgcgc                                                                                  >  1:110796/1‑251 (MQ=255)
                                                                                                                                          ttCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCTCGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTCAGtt                                                                                      <  1:312867/246‑1 (MQ=255)
                                                                                                                                          ttCGCCAAATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGCAATCGATTTCAGtt                                                                                      >  2:312867/1‑246 (MQ=255)
                                                                                                                                                 gATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTCAGTTGCGCAGTCAGGc                                                                          <  1:347489/251‑1 (MQ=255)
                                                                                                                                                  aTCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTCAGTTGCGCAGTCAGGcc                                                                         >  1:148768/1‑251 (MQ=255)
                                                                                                                                                                         ccAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTCAGTTGCGCAGTCAGGCCGTCGATACCTTTGGTACgcagca                                                  <  1:19614/251‑1 (MQ=255)
                                                                                                                                                                           aGTCCACCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGCAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGCGAAATCGATTTCAGTTGCGCAGTCAGGCCGTCGATACCTTTGGTACGCagcagc                                                >  2:141194/1‑251 (MQ=255)
                                                                                                                                                                                 cccGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTCAGTTGCGCAGTCAGGCCGTCGATACCTTTGGTACGCAGCAGCGTTccc                                          >  2:258906/1‑251 (MQ=255)
                                                                                                                                                                                 cccGTATCCACGCGGGAGAAACGGCCCAGATCGATGAAAGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGATATACGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGCAATCGATTTCAGTTGCGCAGTCAGGCCGTCGATACCTTTGGTACGCAGAAGCGTTccc                                          >  2:258207/1‑251 (MQ=255)
                                                                                                                                                                                       tCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTCAGTTGCGCAGTCAGGCCGTCGATACCTTTGGTACGCAGCAGCGTTCACCACTCg                                    >  2:46645/1‑251 (MQ=255)
                                                                                                                                                                                          aCGCGGGAGACACGGCTAAGATCGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTCAGTTGCGCAGTCAGGCCGTCGATACCTTTGGTACGCAGCAGCGTTCCCCACTCGttt                                 <  1:413860/251‑1 (MQ=255)
                                                                                                                                                                                                                cGATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTCAGTTGCGCAGTCAGGCCGTCGATACCTTTGGTACGCAGCAGCGTTCCCCACTCGTTTTGTTTGGTGGTGATCATACTg            >  1:354801/1‑250 (MQ=255)
                                                                                                                                                                                                                 gATGCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAAAGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTCAGTTGCGCAGTCAGGCCGTCGATACCTTTGGTACGCAGCAGCGTTCCCCACTCGTTTTGTTTGGTGGTGATCATACTGAc          >  2:333245/1‑251 (MQ=255)
                                                                                                                                                                                                                    gCAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTCAGTTGCGCAGTCAGGCCGTCGATACCTTTGGTACGCAGCAGCGTTCCCCACTCGTTTTGTTTGGTGGTGATCATACTGACGc        <  2:148768/250‑1 (MQ=255)
                                                                                                                                                                                                                     cAGGTAATCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTCAGTTGCGCAGTCAGGCCGTCGATACCTTTGGTACGCAGCAGCGTTCCCCACTCGTTTTGTTTGGTGGTGATCATACTGACGCCt      <  1:35338/251‑1 (MQ=255)
                                                                                                                                                                                                                        gTAATCGCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTCAGTTGCGCAGTCAGGCCGTCGATACCTTTGGTACGCAGCAGCGTTCCCCACTCGTTTTGTTTGGTGGTGATCATACTGACGCCTTCa   <  2:150153/251‑1 (MQ=255)
                                                                                                                                                                                                                          aaTCTCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGTATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTCAGTTGCGCAGTCAGGCCGTCGATACCTTTGGTACGCAGCAGCGTTCCCCACTCGTTTTGTTTGGTGGTGATCATACTGACGCCTTCAg  <  2:367557/250‑1 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     
TGCCGGACAGGCGGTTAACGCCATATCCGGCCTGAAAAAATTTAACGAGGCAGAACATCAGCAGGCAAATTATACAATTTCGCCAGGGTATACACTTTGTCGTTTACCCCCTGAAGCGTCACATTGTTGCCCTGCTTTTTCGCCAGATCGATAAGATGGAGCAGCAGTGCCAGTCCCCCCGTATCCACGCGGGAGACACGGCTAAGATCGATGCAGGTAATCCCCTTCACCGCTTCCTCACGCATTTCCCAAAGCGGTAGCAAAACGTCCTGATCCAGCTCTCCGGATAACGCCAGCGTGTCACCCGTCTGCATCCAGCTCAGTGACTCGCTCATTATTTTTTCTCTTCCAGAGTGATTTTCTGTTGAGAAATCGATTTCAGTTGCGCAGTCAGGCCGTCGATACCTTTGGTACGCAGCAGCGTTCCCCACTCGTTTTGTTTGGTGGTGATCATACTGACGCCTTCAG  >  NC_000913/3336922‑3337389

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