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.7% 143.5 / NA 44M33I (ATG→ATA) lpxMLipid A biosynthesis myristoyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (12/31);  total (12/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCGCGGCTGAATTTTCTCCGGCCCGCGTATTGCCAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCCATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGC  >  NC_000913/1939849‑1940331
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
aaCGCGCGGCTGAATTTTCTCCGGCCCGCGTATTGCCAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCg                                                                                                                                                                                                                                          >  2:438137/1‑251 (MQ=255)
                ttCTCCGGCCCGCGTATTGCCAACTCAGCCATCATTGCCATCGCTTGCGTCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCTACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTGTGGCGGCGTTAAAGCGATACCCGCTAGCGCTGCTACGCCCAGcc                                                                                                                                                                                                                           <  1:353971/250‑1 (MQ=255)
                         ccGCGTATTGCCAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCCACTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCAAGACGGCGGGAGCTTTTTCCCAGTCGTCCGGCAACACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCTCTATCGCTGCAACGCCCAGCCATGCTcccc                                                                                                                                                                                                                  >  2:564571/1‑250 (MQ=255)
                          cgcgTATTGCCAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGTTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCTGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGTCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAg                                                                                                                                                                                                                <  2:491438/251‑1 (MQ=255)
                            cgTATTGCCAACCCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCGTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGAGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTAGCGCTGCTACGACCAGTCTTGCTCCCCAGTa                                                                                                                                                                                                              <  2:414191/251‑1 (MQ=255)
                             gTATTGCCAACTCAGCCATCATTGCCATCGCTGGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAg                                                                                                                                                                                                             <  2:188369/251‑1 (MQ=255)
                                  gCCAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCgg                                                                                                                                                                                                        <  1:424502/251‑1 (MQ=255)
                                      aCTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGAGTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCGGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTgg                                                                                                                                                                                                    <  2:99630/251‑1 (MQ=255)
                                         cAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGc                                                                                                                                                                                                   <  1:530973/249‑1 (MQ=255)
                                           gCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGaaa                                                                                                                                                                                               <  1:310500/251‑1 (MQ=255)
                                                catTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGAt                                                                                                                                                                                           <  1:108431/250‑1 (MQ=255)
                                                   tGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACCACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTa                                                                                                                                                                                        <  1:7104/250‑1 (MQ=255)
                                                     ccATCGCTGGCGGCGCGGTGGCAAACATCTCATCAACAAGCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCGTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCa                                                                                                                                                                                     <  2:5463/251‑1 (MQ=255)
                                                           cgtGCGGCGCGGGGGCAAACATCTCAGCAAGATTCGCGTCGCGTTCAGCTGCACTACGTGCTGGAAAGCATAGCGACAGATTGATTAACGCACGACGGCGTTAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCa                                                                                                                                                                               <  2:414449/249‑1 (MQ=255)
                                                              gcggcgCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGa                                                                                                                                                                            <  1:181240/251‑1 (MQ=255)
                                                              gcggcgCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCATGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGa                                                                                                                                                                            <  1:231507/251‑1 (MQ=255)
                                                                    cgGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTAt                                                                                                                                                                      >  1:346182/1‑251 (MQ=255)
                                                                           aaaCATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTCTCAAACTCAGGAATGTATTCGCtat                                                                                                                                                               >  2:478495/1‑251 (MQ=255)
                                                                           aaaCATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCtat                                                                                                                                                               >  2:406949/1‑251 (MQ=255)
                                                                            aaCATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCtatt                                                                                                                                                              <  2:346182/251‑1 (MQ=255)
                                                                            aaCATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCtatt                                                                                                                                                              <  1:10851/251‑1 (MQ=255)
                                                                                      caacaaTCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGATAGGATTTAGCAAACTCAGGAATGTATTCGCTATTAtgtttttt                                                                                                                                                     <  2:61391/250‑1 (MQ=255)
                                                                                          aaTCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGGAGCGCGGGTGGCGAAAGGAGTTAGCAAACTCAGGAAAGAAGTCGCTAGTATGTTTTTTCGttt                                                                                                                                                >  1:20398/1‑251 (MQ=255)
                                                                                           aTCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTc                                                                                                                                               <  1:263254/251‑1 (MQ=255)
                                                                                                       tCAGCTTCACTACGTTCTGGTAAGCAGAGCGCCAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGCTACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGCCAGGATTTATCAAACTCAGGAATGTATACGCTATTTTTTTTTTTCGTTTCCATGCTTTTcc                                                                                                                                    >  2:40171/1‑250 (MQ=255)
                                                                                                          gCTTCACGACGTTCTGGAAAGCAGAGCGACAGATTGCTTAACGCACGACGGCGTGAGCTTTTTACCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGACATGCTCCCCAGTAGCGCGGGTGGCGATAGGCTTTATCATACTCAGGAATGTATTAGCTATTATTTTTTTTCGTTTCCATGCTTTTcc                                                                                                                                    >  2:257924/1‑247 (MQ=255)
                                                                                                             tCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTc                                                                                                                              <  1:569935/250‑1 (MQ=255)
                                                                                                                    gTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAgg                                                                                                                      <  1:276222/251‑1 (MQ=255)
                                                                                                                      tCTGGAAAGCTGAGCTACAGATTGATTAACGCCCGACGGCGTGATCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGc                                                                                                                     <  2:487153/250‑1 (MQ=255)
                                                                                                                         ggAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTT‑CGTGTCCATGCTTTTCCAGTTTCGGATAAGGCaaaa                                                                                                                 <  2:190056/250‑1 (MQ=255)
                                                                                                                                             gattAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCTAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTTGTGATAGTGTa                                                                                             <  2:345830/251‑1 (MQ=255)
                                                                                                                                                                     ttttCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCAcc                                                                       <  2:131123/249‑1 (MQ=255)
                                                                                                                                                                     ttttCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTGTCGGAGAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCAcc                                                                       <  2:385466/249‑1 (MQ=255)
                                                                                                                                                                                gtcCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGc                                                          >  1:309954/1‑251 (MQ=255)
                                                                                                                                                                                      aaaaccggcccAGACGTGCCAAATTGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCTCTGCTACGCCCAACCGTGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACCCAGTAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGTCGCAACTTTCCCCGCACCAAATAAAAAAAGCCGGTac                                                    <  2:153353/243‑1 (MQ=255)
                                                                                                                                                                                               ccAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGt                                           >  1:263695/1‑251 (MQ=255)
                                                                                                                                                                                                 aGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTAc                                         <  1:580911/251‑1 (MQ=255)
                                                                                                                                                                                                               gggATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGa                           <  1:406949/251‑1 (MQ=255)
                                                                                                                                                                                                                    cACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGtta                      <  1:176930/251‑1 (MQ=255)
                                                                                                                                                                                                                      cGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATGCGGCGCCGCTCCCTCTCGTGTTAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGttaat                    <  2:20398/251‑1 (MQ=255)
                                                                                                                                                                                                                          aCTTTGGCGGCGTTAAAGCGATTCCCGCTATCGCTGCTACGCCCAGCCATGCTCCTCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATATTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTtaattaat                <  1:32197/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cgcgTTAAAGCGATACCCGCCAGCGCTGCTTCGCCCCGCCATCTTACCCAGTAGCGCGGGTGGCGAACGTCTTTATCCAACGCAGGAATGTATTCGCTATTAATTTTTTTCGTTCCCAATCTTTTCCAGTGTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGTTACTGGCTGCGTACCGGCTGCGAATGTATGTTAATTAATCAAACCGt        <  2:397396/250‑1 (MQ=255)
                                                                                                                                                                                                                                    cgTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGt        >  1:299913/1‑249 (MQ=255)
                                                                                                                                                                                                                                        aaaGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGAAGCTGc  >  1:199791/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
AACGCGCGGCTGAATTTTCTCCGGCCCGCGTATTGCCAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCCATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGC  >  NC_000913/1939849‑1940331

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