Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,796,254 C→T D92D (GAC→GAT 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,2540CT98.4% 211.3 / ‑4.5 64D92D (GAC→GATwaaCADP‑heptose:LPS heptosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (24/39);  total (24/40)
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.42e-01

ATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCC  >  NC_000913/3796016‑3796499
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     
aTGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTg                                                                                                                                                                                                                                           >  2:363267/1‑251 (MQ=255)
  gTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGt                                                                                                                                                                                                                                         >  1:311138/1‑251 (MQ=255)
       ccATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTaaaaa                                                                                                                                                                                                                                    <  1:217794/251‑1 (MQ=255)
           aCGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAgcg                                                                                                                                                                                                                                 <  2:123215/250‑1 (MQ=255)
             gTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCgcgg                                                                                                                                                                                                                              <  2:313254/251‑1 (MQ=255)
               tGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCgcggc                                                                                                                                                                                                                             <  1:201402/250‑1 (MQ=255)
                     cactcactGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTg                                                                                                                                                                                                                      >  1:157048/1‑251 (MQ=255)
                      actcactGATGCCCCGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGATCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTg                                                                                                                                                                                                                      >  1:69861/1‑250 (MQ=255)
                         cactGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTGTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGAcg                                                                                                                                                                                                                   <  1:66701/250‑1 (MQ=255)
                            tGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTc                                                                                                                                                                                                               >  1:416528/1‑251 (MQ=255)
                            ggATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGGTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTc                                                                                                                                                                                                               <  2:311138/250‑1 (MQ=255)
                             gATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCt                                                                                                                                                                                                              <  1:348449/251‑1 (MQ=255)
                              aTGCCCAGCAGGCAAGCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTGGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTg                                                                                                                                                                                                             <  1:390603/251‑1 (MQ=255)
                                          cAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGCTTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGGGACGCGTCTGGCGCAAGGCGt                                                                                                                                                                                                  >  2:388341/1‑250 (MQ=255)
                                              cccAGGGATTAAGTTTGACTGGTTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATGCGTCGCTTTCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAGCTATGACGCAGTTAGCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAgc                                                                                                                                                                                             <  1:416469/251‑1 (MQ=255)
                                                      ttAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGAACAGATTCCTTGCTGGCACGCTGCCGTTGAGAGAGTTATTCCTGTGGCAATAAGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGACGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAGgcatggattg                                                                                                                                                                                     >  1:418621/1‑247 (MQ=255)
                                                       tAAGTTTGATTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAgcatggattg                                                                                                                                                                                     <  2:78550/250‑5 (MQ=255)
                                                       tAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAgcatggattg                                                                                                                                                                                     <  1:90498/250‑5 (MQ=255)
                                                        aaGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGCGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGACCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATAGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAgcatggattgg                                                                                                                                                                                    >  1:427394/1‑245 (MQ=255)
                                                             tGACTGGGTGGTGTAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCGGTGGCAATACGTCGCTGGCGTAAATCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCGCAGGGGCTGTTAAAAAGCGCGGCGCGGGTGACGCGTCTGGCGCAGGGCGTAAAGCATTGATTGGACTgg                                                                                                                                                                               <  2:434070/250‑1 (MQ=255)
                                                                  gggTGGTGGAAGAAGGGTTCGCACAGATTCCGTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAAcc                                                                                                                                                                         <  2:190962/251‑1 (MQ=255)
                                                                    gtggtgGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGc                                                                                                                                                                       >  1:128292/1‑251 (MQ=255)
                                                                      ggtggAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGGCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTc                                                                                                                                                                     <  2:326800/251‑1 (MQ=255)
                                                                                gggTGCGCACAGAGTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCGGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCTGCCCCCATAAAAGCTGAACGCATAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGTTGACGCGTCGGGCGCATGGCGTAAAGCATGGATTGGACTGGCAATCCGCTCGCGAAccgtt                                                                                                                                                           <  2:37955/251‑4 (MQ=255)
                                                                                      gCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGCTGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTagccac                                                                                                                                                     >  2:330714/1‑250 (MQ=255)
                                                                                      gCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTagcc                                                                                                                                                       >  1:66211/1‑249 (MQ=255)
                                                                                       cacaGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTagccagc                                                                                                                                                    >  2:258296/1‑251 (MQ=255)
                                                                                         caGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAATCCGCTCGCTAACCTTTagccagc                                                                                                                                                    <  2:398355/249‑1 (MQ=255)
                                                                                          aGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTg                                                                                                                                                 <  1:207815/251‑1 (MQ=255)
                                                                                            aTTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGtt                                                                                                                                               >  1:222317/1‑251 (MQ=255)
                                                                                             ttccttccTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGttt                                                                                                                                              >  1:93760/1‑251 (MQ=255)
                                                                                                  tcctGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACa                                                                                                                                          <  2:66211/250‑1 (MQ=255)
                                                                                                   cctGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACaa                                                                                                                                         <  1:258296/250‑1 (MQ=255)
                                                                                                     tGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAgcatgg                                                                                                                                                                                         <  1:201291/200‑1 (MQ=255)
                                                                                                     tGGCACGCTGCCGTTGAGAGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAgcatgg                                                                                                                                                                                         >  2:201291/1‑200 (MQ=255)
                                                                                                      ggCACGCTGCCGTTGATCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTTGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTGCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTTGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGTCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGt                                                                                                                                     <  1:173339/251‑1 (MQ=255)
                                                                                                      ggCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGt                                                                                                                                     >  1:90213/1‑251 (MQ=255)
                                                                                                                   ggAGCGAGTTATTCCTGTGGCAATACGTCGCTGGTGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAATCATCATATTg                                                                                                                        <  1:102405/250‑1 (MQ=255)
                                                                                                                      gtgAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCaa                                                                                                                     <  2:112926/249‑1 (MQ=255)
                                                                                                                          gTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAAca                                                                                                                 <  2:222317/251‑1 (MQ=255)
                                                                                                                                         aTACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGa                                                                                                  <  2:128292/251‑1 (MQ=255)
                                                                                                                                           aCGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAAc                                                                                                <  2:22249/251‑1 (MQ=255)
                                                                                                                                                gCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCAcc                                                                                           <  1:363267/251‑1 (MQ=255)
                                                                                                                                                   ggCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTTGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATAGTGCAAAACAGCAGCACGCCGTAGAGCGCAccc                                                                                          <  2:254669/249‑1 (MQ=255)
                                                                                                                                                   ggCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCcg                                                                                         <  1:205151/250‑1 (MQ=255)
                                                                                                                                                   ggCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACGATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCAccc                                                                                          <  1:43721/249‑1 (MQ=255)
                                                                                                                                                      gTAAAGCGTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCATTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGGCTGTCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCTCTCTCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGaa                                                                                     <  1:267040/251‑1 (MQ=255)
                                                                                                                                                                       gCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGtct                                                                    >  1:415183/1‑249 (MQ=255)
                                                                                                                                                                           ccATAAACGCTGAACGCAACGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGc                                                                >  1:318111/1‑251 (MQ=255)
                                                                                                                                                                                  agctgAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAgc                                                          >  1:75312/5‑250 (MQ=255)
                                                                                                                                                                                  agctgAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAgc                                                          >  1:174758/5‑250 (MQ=255)
                                                                                                                                                                                        aCGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAgcaaaccgc                                                   >  2:347232/1‑251 (MQ=255)
                                                                                                                                                                                           cAAAGCGAGTCGTAAAGCGCTACAAGCAGAGAACTATGACGCAGTTAGCGATGCTCAGGGGCTGGTAGAAAGCGGGGCGCTGGTGACGCGGCTGGCGCATGGCGTAACGCATGGATTGTACTGGCAAACCGCACGGGAGCCGTGAGCCAGCCTGTGTGGTAATCGGAAGCATCATGTTGCAAATCAGGAGCACGCAGGAGAACGCACGCGCGAACGGTTGGCCAAAAGTGTGGGCTATTgcaatcggcaaa                                                <  2:23009/251‑1 (MQ=255)
                                                                                                                                                                                                 gTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCACACCCAgg                                          <  1:401801/251‑1 (MQ=255)
                                                                                                                                                                                                       ggAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCGCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGAtt                                    <  2:157048/250‑1 (MQ=255)
                                                                                                                                                                                                             gcTACAAGCGGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAACAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCCAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCt                               >  2:181613/1‑250 (MQ=255)
                                                                                                                                                                                                              cTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTAt                             >  1:335658/1‑251 (MQ=255)
                                                                                                                                                                                                                aCAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTAttg                           <  1:180347/251‑3 (MQ=255)
                                                                                                                                                                                                                   aGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGcac                        <  2:415183/251‑1 (MQ=255)
                                                                                                                                                                                                                    gCAGAGAACTATGACGCAGTTAGCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGcaca                       <  2:376643/251‑1 (MQ=255)
                                                                                                                                                                                                                        agaACTATGACGCAGTTATCGATGCACAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGGAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAGCCCAGGGCGATTATGCTATTGCACAGCAt                   <  2:416528/251‑1 (MQ=255)
                                                                                                                                                                                                                              aTGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCAtttcct              <  2:287234/250‑4 (MQ=255)
                                                                                                                                                                                                                                gACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGGGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACTCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGa            <  1:187813/250‑1 (MQ=255)
                                                                                                                                                                                                                                         aTCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGAAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGcc  >  1:51684/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     
ATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCC  >  NC_000913/3796016‑3796499

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: