Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,796,011 G→C 5.9% M11I (ATG→ATC waaC → ADP‑heptose:LPS heptosyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,796,0110GC 5.9% 804.7 / 27.3 236M11I (ATG→ATCwaaCADP‑heptose:LPS heptosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (104/118);  new base C (5/9);  total (109/127)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

Only 100 of 244 total aligned reads displayed.

CACAAAGTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGT  >  NC_000913/3795862‑3796158
                                                                                                                                                     |                                                                                                                                                   
cacaAAGTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATg                                                                                                                                                     >  2:1811495/1‑150 (MQ=255)
 acaAAGTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATgg                                                                                                                                                    >  2:1183029/1‑150 (MQ=255)
   aaaGTTCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCATGACGGATGCGTGTTTTGATCGTTAAAACATCGTCGATGGGc                                                                                                                                                  <  2:812142/150‑1 (MQ=255)
    aaGTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCg                                                                                                                                                 <  1:869442/150‑1 (MQ=255)
      gTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGTAATAACTCAACGCGCTATTGTTAAAAGAGGAAGCCTGACGGATGCTTGTTTTGATGGTTAAAATATCGTCCATCGACGAt                                                                                                                                               <  2:883628/150‑1 (MQ=255)
         cGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGtt                                                                                                                                            >  2:1567644/1‑150 (MQ=255)
           tAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTtct                                                                                                                                          >  1:301162/1‑150 (MQ=255)
             aaGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTcc                                                                                                                                        <  2:1258722/150‑1 (MQ=255)
              aGGTGACGCTGTGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGCCCATCGACGATGttttaac                                                                                                                                       <  1:1099090/150‑6 (MQ=255)
              aGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCa                                                                                                                                       <  1:1463797/150‑1 (MQ=255)
              aGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGCCCATCGACGATGttttaac                                                                                                                                       <  2:722135/150‑6 (MQ=255)
               ggTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCAt                                                                                                                                      <  1:198066/150‑1 (MQ=255)
                gTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGCCCATCGACGATGttttaacga                                                                                                                                     <  1:538455/150‑8 (MQ=255)
                  gACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACg                                                                                                                                   >  1:1929038/1‑150 (MQ=255)
                       tGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGcc                                                                                                                              >  2:1362407/1‑150 (MQ=255)
                       tGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGcc                                                                                                                              <  1:1200512/150‑1 (MQ=255)
                           gAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGca                                                                                                                          <  2:1743083/150‑1 (MQ=255)
                           gAGGGTTATCACCAGAGATTAATCGACATTACTCCCCATCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTTACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGca                                                                                                                          <  2:1545577/150‑1 (MQ=255)
                             ggtttATCACCAGAGCTTATTCTACATTACTCCCCAGCGCGTACTGGAAGAACTCATCGCGCTATTGTTAGAAGAGGAAGCCTGACGTATTCGGGTTTTGATCTTTAAAAGATCGTGGATGGGCGATGTTCTCCATACGTTGCCCGcact                                                                                                                        <  2:1602293/147‑1 (MQ=255)
                              ggTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGCACTGGAAGAAATCAACGCGCTATTGTTACAAGACGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACAACGTCGATGGGCGATCTTATCCATACGTTGCCCGccctc                                                                                                                       >  2:1671416/1‑150 (MQ=255)
                                   tCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGa                                                                                                                  <  1:1650134/150‑1 (MQ=255)
                                   tCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGa                                                                                                                  <  1:594992/150‑1 (MQ=255)
                                   tCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGa                                                                                                                  <  1:338067/150‑1 (MQ=255)
                                     aCCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATg                                                                                                                <  2:264636/150‑1 (MQ=255)
                                     aCCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATg                                                                                                                >  2:572195/1‑150 (MQ=255)
                                     aCCAGAGCTGAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGACGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATg                                                                                                                >  2:932265/1‑150 (MQ=255)
                                      ccAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGc                                                                                                               <  1:1459724/150‑1 (MQ=255)
                                        agagCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGccc                                                                                                             <  2:1483204/150‑1 (MQ=255)
                                           gCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCcagc                                                                                                          <  2:1395960/150‑1 (MQ=255)
                                           gCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCcagc                                                                                                          >  1:267833/1‑150 (MQ=255)
                                            cTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCcagca                                                                                                         >  2:1947945/1‑150 (MQ=255)
                                            cTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCcagca                                                                                                         >  1:1299456/1‑150 (MQ=255)
                                             ttAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCcagcag                                                                                                        <  1:1533752/150‑1 (MQ=255)
                                              tAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAgg                                                                                                       >  2:1697792/1‑150 (MQ=255)
                                                aTCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCa                                                                                                     >  2:468826/1‑150 (MQ=255)
                                                    aCATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATcc                                                                                                 >  1:377423/1‑150 (MQ=255)
                                                        tACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAgg                                                                                             >  1:622101/1‑150 (MQ=255)
                                                            ccccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGAtt                                                                                         >  1:781895/1‑150 (MQ=255)
                                                            ccccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGAtt                                                                                         >  1:1773026/1‑150 (MQ=255)
                                                            ccccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGAtt                                                                                         >  2:1996066/1‑150 (MQ=255)
                                                            ccccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGAtt                                                                                         >  1:752748/1‑150 (MQ=255)
                                                                aGCGCCTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGGTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGt                                                                                     <  1:166036/150‑1 (MQ=255)
                                                                  cgcgTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGttt                                                                                   <  1:1364384/150‑1 (MQ=255)
                                                                     gTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGAc                                                                                <  1:1736531/150‑1 (MQ=255)
                                                                        cTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGtt                                                                                    >  2:322661/1‑143 (MQ=255)
                                                                           gaagaaCTCAACGCGCTATTATTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggtg                                                                          >  1:1588574/1‑150 (MQ=255)
                                                                            aagaaCTCAACGCGCTATTGTTACACGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTTTCTATAGGTTTCCGGCACTCACTGATTCTAATCAGGCATTCCCTGTGATTAATTTTGACTGtgtgg                                                                         <  2:1775365/150‑1 (MQ=255)
                                                                            aagaaCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggtgg                                                                         <  1:114739/150‑1 (MQ=255)
                                                                                aCTCAACGCGCTTTTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGTTGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGga                                                                     >  2:1343033/1‑150 (MQ=255)
                                                                                aCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTGACTGATGcc                                                                                                              <  2:1894872/109‑1 (MQ=255)
                                                                                aCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTGACTGATGcc                                                                                                              >  1:1894872/1‑109 (MQ=255)
                                                                                 cTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaa                                                                    <  2:1724815/150‑1 (MQ=255)
                                                                                 cTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTCGTGgaa                                                                    >  1:1009925/1‑150 (MQ=255)
                                                                                  tCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaag                                                                   >  1:980647/1‑150 (MQ=255)
                                                                                          cTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAgg                                                               <  2:216213/146‑1 (MQ=255)
                                                                                          cTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAgg                                                               >  1:209079/1‑146 (MQ=255)
                                                                                            aTTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGc                                                         >  2:1294088/1‑150 (MQ=255)
                                                                                             ttGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGca                                                        <  2:1053156/150‑1 (MQ=255)
                                                                                              tGTTACAAGTGGAAGCCTGACGGATGCTGGTTTTGATCGTTAAAACTTCGTCGATGGGCGATGTTCTCCATTCGTTGCGCGCACTCACTGATGCCCAGCAGTCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAATGGTTCGcac                                                       <  2:781895/150‑1 (MQ=255)
                                                                                              tGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaa                                                                    <  2:424498/137‑1 (MQ=255)
                                                                                              tGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaa                                                                    >  1:424498/1‑137 (MQ=255)
                                                                                              tGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGcac                                                       <  2:301162/150‑1 (MQ=255)
                                                                                               gTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGcaca                                                      <  2:840302/150‑1 (MQ=255)
                                                                                               gTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGcaca                                                      >  2:1877443/1‑150 (MQ=255)
                                                                                                 tACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggt                                                                           <  1:218565/127‑1 (MQ=255)
                                                                                                 tACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggt                                                                           >  2:218565/1‑127 (MQ=255)
                                                                                                 tACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGa                                                    >  1:1573984/1‑150 (MQ=255)
                                                                                                  aCAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGc                                                                                                      <  1:1719399/99‑1 (MQ=255)
                                                                                                  aCAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATccc                                                                                                <  1:832686/105‑1 (MQ=255)
                                                                                                  aCAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAt                                                   >  1:1402269/1‑150 (MQ=255)
                                                                                                  aCAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAt                                                   <  1:1513204/150‑1 (MQ=255)
                                                                                                    aaGAGGAAGCCTGACGGATGCGGGGTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttc                                                 <  2:1445946/150‑1 (MQ=255)
                                                                                                     agagGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttcc                                                >  2:6186/1‑150 (MQ=255)
                                                                                                     agagGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttcc                                                <  1:1783115/150‑1 (MQ=255)
                                                                                                         gAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttccttcc                                            <  1:879055/150‑1 (MQ=255)
                                                                                                          aaGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATtccttcct                                           <  1:105971/150‑1 (MQ=255)
                                                                                                           aGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTg                                          >  1:556404/1‑150 (MQ=255)
                                                                                                              cTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCa                                       >  1:1969042/1‑150 (MQ=255)
                                                                                                              cTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCa                                       <  1:1490997/150‑1 (MQ=255)
                                                                                                                        cGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGtt                             <  2:1117382/150‑1 (MQ=255)
                                                                                                                        cGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGtt                             >  1:653241/1‑150 (MQ=255)
                                                                                                                           gTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAg                          >  2:1158727/1‑150 (MQ=255)
                                                                                                                           gTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAg                          <  2:1622013/150‑1 (MQ=255)
                                                                                                                             tttGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCg                        <  2:1781229/150‑1 (MQ=255)
                                                                                                                               tGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTg                                          >  1:203308/1‑130 (MQ=255)
                                                                                                                               tGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTg                                          <  2:203308/130‑1 (MQ=255)
                                                                                                                               tGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAg                      >  1:1045140/1‑150 (MQ=255)
                                                                                                                                gATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGt                     <  1:6186/150‑1 (MQ=255)
                                                                                                                                gATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGt                     <  2:1045140/150‑1 (MQ=255)
                                                                                                                                gATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGt                     <  1:1294088/150‑1 (MQ=255)
                                                                                                                                gATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGt                     >  1:443918/1‑150 (MQ=255)
                                                                                                                                gATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGt                     >  1:1848487/1‑150 (MQ=255)
                                                                                                                                   cGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTAt                  >  2:62446/1‑150 (MQ=255)
                                                                                                                                     ggagaaCATCGCCCATCGACGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTc                >  2:1278889/5‑150 (MQ=255)
                                                                                                                                      tAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTcc               <  1:1275212/150‑1 (MQ=255)
                                                                                                                                      gagaaCATCGCCCATCGACGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTcc               >  2:1135360/4‑150 (MQ=255)
                                                                                                                                            aTCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGc         <  2:549837/150‑1 (MQ=255)
                                                                                                                                            aTCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGc         >  1:1085090/1‑150 (MQ=255)
                                                                                                                                                 cgATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATAc    >  2:758687/1‑150 (MQ=255)
                                                                                                                                                 cgATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGGGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATAc    <  2:1929038/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CACAAAGTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGT  >  NC_000913/3795862‑3796158

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: