Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,940,095 C→T M33I (ATG→ATA lpxM ← Lipid A biosynthesis myristoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,940,0950CT97.6% 139.0 / ‑3.9 42M33I (ATG→ATAlpxMLipid A biosynthesis myristoyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (20/21);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.84e-01

CGCGCGGCTGAATTTTCTCCGGCCCGCGTATTGCCAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCCATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCT  >  NC_000913/1939851‑1940341
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
cgcgcgGCTGAATTTTCTCCGGCCCGCGTATTGCCAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGACAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCgct                                                                                                                                                                                                                                                  <  1:154263/251‑1 (MQ=255)
                                        aGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGCTACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGACGCTCGGGTGGCgc                                                                                                                                                                                                           >  2:468958/1‑249 (MQ=255)
                                         gCCACCATTGCCATCGCGTGCGGCGCGGTGGCAAACATCTCAACAACAACCGCTCCGCGTCCAGCTTCCCTACGTCCTGGAAAGCAGAGCGCCAGATTGATTAACGCACGACGGCGTGAGCTGGTTCCCAGTCGTCCGGCAAAACTGCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTgtagaaa                                                                                                                                                                                                         <  2:483497/251‑7 (MQ=255)
                                               atTGCCATCGCTTGCGGCGCGGTGGCAAACATCGCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGAttt                                                                                                                                                                                                   <  1:476930/251‑1 (MQ=255)
                                                               gcgcGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATg                                                                                                                                                                                   <  1:571529/251‑1 (MQ=255)
                                                                      ggCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCg                                                                                                                                                                             <  1:35663/250‑1 (MQ=255)
                                                                        cAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCtt                                                                                                                                                                          >  1:270797/1‑250 (MQ=255)
                                                                                tcatcaACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTCttttt                                                                                                                                                                  >  1:398824/1‑251 (MQ=255)
                                                                                  atcaACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTAttttttt                                                                                                                                                                <  1:246707/251‑1 (MQ=255)
                                                                                      acaaTCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTTGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGt                                                                                                                                                            <  1:62007/251‑1 (MQ=255)
                                                                                        aaTCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGttt                                                                                                                                                          <  1:443714/251‑1 (MQ=255)
                                                                                         aTCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGCGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGAGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTc                                                                                                                                                         >  1:205802/1‑251 (MQ=255)
                                                                                                    ttCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAgg                                                                                                                                                                                                       <  2:71909/194‑1 (MQ=255)
                                                                                                    ttCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAgg                                                                                                                                                                                                       >  1:71909/1‑194 (MQ=255)
                                                                                                      cAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTcc                                                                                                                                              >  1:141573/1‑249 (MQ=255)
                                                                                                            cccatggTTCTGGAAAGCAGTGCTACAGATTGAATAACGCACGACGGCGTGTGTATTTCCCCATTCGTCCGGCAAAACGCCCCAGCCGTGCGAGAATGGGATCACGGAACTGTGGCGGCGTTAAAGCGATACCCGGTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCTGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATGCGCTATTATGTGTTTGCGTGTCCATGCTTTTGCAGTTTCgg                                                                                                                                      <  2:291611/245‑1 (MQ=255)
                                                                                                             aCTACGATCTGGAAAGCAGAGCGTCAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTACAGCACAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTCAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCAACAGAAGCGCGGGTGGCGAACGGCTTTATCAAACTCAGGACAGTATTCGATATTATTTTTTTTAGTTTCCATGCTTTTCCAAtttagtt                                                                                                                                     >  2:375483/1‑247 (MQ=255)
                                                                                                               tACGTTCTGGAAAGAGGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGCTTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGAt                                                                                                                                    >  1:424286/1‑250 (MQ=255)
                                                                                                                 cGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGACCTTTGGCGTCGTTAAAGCGTTACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTCTTATTTTTTTTCGATTCCATGCTTTTCCCGTTTCGGATacg                                                                                                                                 >  2:392888/1‑249 (MQ=255)
                                                                                                                 cGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAg                                                                                                                                 <  1:232594/251‑1 (MQ=255)
                                                                                                                      tGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATg                                                                                                                                                                                   <  2:365095/196‑1 (MQ=255)
                                                                                                                      tGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATg                                                                                                                                                                                   >  1:365095/1‑196 (MQ=255)
                                                                                                                         aaaGCAGAGCGACAGATTGATTAACGCACGACGGCGTGGGCTTTTTCCCAGCCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATCCCCGCCATCGCTGCTACGCCCAGCCATGCTCCACAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTGTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAaat                                                                                                                         <  1:144351/251‑1 (MQ=255)
                                                                                                                           aGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAaatc                                                                                                                        >  1:527583/1‑250 (MQ=255)
                                                                                                                                    aCAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGa                                                                                                              <  1:17702/251‑1 (MQ=255)
                                                                                                                                     atgattgagtAACGCACGACGGGGGTAGCTGTCTCCCAGTCGTTCGGCAAAACGGCCCAGCCGTGCCAGAATGGGATCACGGAACTGTGGCGGCGTTAAAGCGATACCCGCTAGCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTAGCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTGTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGAt                                                                                                             <  2:205802/249‑1 (MQ=255)
                                                                                                                                                aCGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAgcggc                                                                                                  <  1:265952/251‑1 (MQ=255)
                                                                                                                                                  gCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTa                                                                                                                                                                                                                        >  1:128719/1‑131 (MQ=255)
                                                                                                                                                  gCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTa                                                                                                                                                                                                                        <  2:128719/131‑1 (MQ=255)
                                                                                                                                                    acgacgGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGACTGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCa                                                                                               >  1:478085/1‑250 (MQ=255)
                                                                                                                                                            gTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGAGTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGTTAAGGCAGAAATCCATCTGGTGATAGTGCAACGGCGCAACTTGcccc                                                                                      >  2:479936/1‑251 (MQ=255)
                                                                                                                                                                     ttCCCAGTCGTCCGGCAAAACGTCCCAGACGTGCCAGAATGGGATCACGGAAATTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGAGCGGGTGGCGAAAGGCTTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCATTTCCCTGCTTTTCAAGTTTCGGAAAAGACAAAAATCAATATGGTGATCGTGTAGCGGCGCCACTTTcccacccccaaat                                                                             >  2:478873/1‑241 (MQ=255)
                                                                                                                                                                                  ggCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGt                                                                <  1:348451/251‑1 (MQ=255)
                                                                                                                                                                                        aCGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTactgac                                                          <  1:255998/251‑1 (MQ=255)
                                                                                                                                                                                                         gAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGa                                          >  1:114058/1‑250 (MQ=255)
                                                                                                                                                                                                                     ggAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGCATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGCCTGCGTACCGGCTGCGAATGGATGttaatt                             >  1:324929/1‑251 (MQ=255)
                                                                                                                                                                                                                        aCTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTtaattaat                          <  1:392888/251‑1 (MQ=255)
                                                                                                                                                                                                                           ttGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCa                        <  1:490478/250‑1 (MQ=255)
                                                                                                                                                                                                                           ttGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCACCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCCAATAAAAAAAGCCGGTACTGACTGCGAACCGGCTGCGAAGGGATGTTAATTAATCa                        >  2:490478/1‑250 (MQ=255)
                                                                                                                                                                                                                             ggcggcGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCGCCCCATTAGCGCGGTTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTCTTCGTTTCCATGCTTTTCCAGTTTCGGATAACGCAAAAATCAATCTGGTGCTAGTGTAGCGGCGCAACTTGCCCCGCATCAAATAAAAAAAGCCGGTACTGACTGTGTACCGGCTGCGAATGGATGTTAACTAATCaaa                      <  2:415390/250‑1 (MQ=255)
                                                                                                                                                                                                                                              aCCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAAt    >  1:382495/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ccGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAAtct  >  1:157974/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
CGCGCGGCTGAATTTTCTCCGGCCCGCGTATTGCCAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCCATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCT  >  NC_000913/1939851‑1940341

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