Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,012,478 T→G 5.1% intergenic (+85/+136) yqeA → / ← yqeB putative amino acid kinase YqeA/XdhC‑CoxI family protein YqeB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,012,4780TG 5.1% 588.8 / 17.3 193intergenic (+85/+136)yqeA/yqeBputative amino acid kinase YqeA/XdhC‑CoxI family protein YqeB
Reads supporting (aligned to +/- strand):  ref base T (94/89);  new base G (5/5);  total (99/94)
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.77e-01

Only 100 of 197 total aligned reads displayed.

TTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCA  >  NC_000913/3012329‑3012627
                                                                                                                                                     |                                                                                                                                                     
   ccGTATTGAGCAAAGCGAAAGAAAGGATTGAAGGCAAAACAGGAACGGTGATTTTTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGAGCTTTGTTGATTCCGGATGAGGCGCGACCGCCTTAt                                                                                                                                                    >  2:1049899/1‑150 (MQ=255)
         tGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTcc                                                                                                                                              >  1:1056045/1‑150 (MQ=255)
          gAGCAAAGCGAAAGAAGGGATTTAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCt                                                                                                                                             <  2:1623504/150‑1 (MQ=255)
          gAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCt                                                                                                                                             >  2:1266657/1‑150 (MQ=255)
          gAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATAAGGCGCTCGCGCCGCATCCggcat                                                                                                                                             <  1:1435077/150‑5 (MQ=255)
           aGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTa                                                                                                                                            >  2:703293/1‑150 (MQ=255)
           aGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTa                                                                                                                                            >  2:323536/1‑150 (MQ=255)
            gCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTAc                                                                                                                                           >  1:1939587/1‑150 (MQ=255)
            gCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTAc                                                                                                                                           <  1:1136282/150‑1 (MQ=255)
             cAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGTGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACa                                                                                                                                          <  2:1253795/150‑1 (MQ=255)
              aaaGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACaa                                                                                                                                         >  2:1755024/1‑150 (MQ=255)
                 gCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAAc                                                                                                                                      >  1:2133478/1‑150 (MQ=255)
                   gAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACAt                                                                                                                                    >  2:1131805/1‑150 (MQ=255)
                    aaaGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATc                                                                                                                                   <  2:1216439/150‑1 (MQ=255)
                     aagaagGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCg                                                                                                                                  >  1:1410535/1‑150 (MQ=255)
                       gaagGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCa                                                                                                                                >  1:543090/1‑150 (MQ=255)
                             aTTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCa                                                                                                                          <  2:2036806/150‑1 (MQ=255)
                               tGAAGGCAGAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAAt                                                                                                                        <  2:529807/150‑1 (MQ=255)
                                gAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGAGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATa                                                                                                                       >  2:853883/1‑150 (MQ=255)
                                  aGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGa                                                                                                                     >  2:1232876/1‑150 (MQ=255)
                                     cAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTg                                                                                                                  >  2:921166/1‑150 (MQ=255)
                                      aaaaCAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGc                                                                                                                 >  1:411393/1‑150 (MQ=255)
                                        aaCAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCaa                                                                                                               <  1:1450778/150‑1 (MQ=255)
                                        aaCAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCaa                                                                                                               >  2:606502/1‑150 (MQ=255)
                                          cAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATg                                                                                                             >  2:1839236/1‑150 (MQ=255)
                                           aGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGc                                                                                                            <  1:1963716/150‑1 (MQ=255)
                                            ggAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCa                                                                                                           >  1:1043450/1‑150 (MQ=255)
                                            ggAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCa                                                                                                           <  1:1999542/150‑1 (MQ=255)
                                                cGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGGGCTTg                                                                                                       >  1:1102854/1‑150 (MQ=255)
                                                 ggTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGt                                                                                                      >  2:1549974/1‑150 (MQ=255)
                                                   tGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAg                                                                                                    <  2:2299239/150‑1 (MQ=255)
                                                      tttGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGcc                                                                                                 >  1:1803240/1‑150 (MQ=255)
                                                      tttGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGcc                                                                                                 >  1:1794232/1‑150 (MQ=255)
                                                           cAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATa                                                                                            >  1:1321404/1‑150 (MQ=255)
                                                           cAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATa                                                                                            <  2:1056045/150‑1 (MQ=255)
                                                           cAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATa                                                                                            <  1:606502/150‑1 (MQ=255)
                                                             gTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAg                                                                                          <  2:1501948/150‑1 (MQ=255)
                                                             gTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAg                                                                                          >  1:1041230/1‑150 (MQ=255)
                                                             gTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAg                                                                                          >  1:710176/1‑150 (MQ=255)
                                                             gTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAg                                                                                          >  2:977624/1‑150 (MQ=255)
                                                             gTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAg                                                                                          >  1:1352969/1‑150 (MQ=255)
                                                              tAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGc                                                                                         >  1:463329/1‑150 (MQ=255)
                                                              tAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGc                                                                                         <  2:543090/150‑1 (MQ=255)
                                                               aaTTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCg                                                                                        <  2:1429270/150‑1 (MQ=255)
                                                               aaTTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCg                                                                                        <  2:240236/150‑1 (MQ=255)
                                                                aTTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGt                                                                                       <  1:230258/150‑1 (MQ=255)
                                                                aTTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGt                                                                                       >  1:1332651/1‑150 (MQ=255)
                                                                  ttAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAg                                                                                     <  1:1687412/150‑1 (MQ=255)
                                                                    aaaGGCCGGGGGGGTAATACCCTCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGAGGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAgcg                                                                                   <  2:2250232/150‑1 (MQ=255)
                                                                          cGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                             <  2:2372522/150‑1 (MQ=255)
                                                                            ggtgggtAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                           <  1:402712/150‑1 (MQ=255)
                                                                             gtgggtAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                          <  2:2075461/150‑1 (MQ=255)
                                                                                 gtAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAAttt                                                                      <  2:2303881/150‑1 (MQ=255)
                                                                                   aaTACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTgg                                                                    <  1:1755024/150‑1 (MQ=255)
                                                                                        cATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGttt                                                               >  2:2365654/1‑150 (MQ=255)
                                                                                        cATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGttt                                                               >  1:317574/1‑150 (MQ=255)
                                                                                            cGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGtca                                                           >  1:1900654/1‑150 (MQ=255)
                                                                                             ggCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGtcat                                                          >  2:1155206/1‑150 (MQ=255)
                                                                                              gCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGtcatc                                                         >  2:800460/1‑150 (MQ=255)
                                                                                                cAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAg                                                       >  2:1645552/1‑150 (MQ=255)
                                                                                                cAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAg                                                       <  2:1457100/150‑1 (MQ=255)
                                                                                                 aaaGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGt                                                      <  1:2071133/150‑1 (MQ=255)
                                                                                                   aGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGtct                                                    >  2:1400342/1‑150 (MQ=255)
                                                                                                    gATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGtctc                                                   <  1:1141301/150‑1 (MQ=255)
                                                                                                      ttGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAg                                                 >  2:504835/1‑150 (MQ=255)
                                                                                                       tGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTAATCAGTCTCAgg                                                <  2:1706792/150‑1 (MQ=255)
                                                                                                        gTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGc                                               >  2:879009/1‑150 (MQ=255)
                                                                                                        gTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGc                                               >  2:172462/1‑150 (MQ=255)
                                                                                                        gTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGAGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGc                                               >  2:1870597/1‑150 (MQ=255)
                                                                                                          cAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCg                                             >  1:2318708/1‑150 (MQ=255)
                                                                                                           aagtaagtGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                            >  1:2256987/1‑150 (MQ=255)
                                                                                                            agtaagtGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCggg                                           >  2:375173/1‑150 (MQ=255)
                                                                                                              taagtGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtg                                         <  1:853883/150‑1 (MQ=255)
                                                                                                                agtGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtggg                                       <  1:963334/150‑1 (MQ=255)
                                                                                                                 gggcgcTTTGTCGACGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtgggt                                      <  2:1011716/148‑1 (MQ=255)
                                                                                                                   gcgcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTaa                                    >  2:2315894/1‑150 (MQ=255)
                                                                                                                    cgcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAAt                                   <  1:502192/150‑1 (MQ=255)
                                                                                                                     gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATa                                  >  1:2475772/1‑150 (MQ=255)
                                                                                                                     gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATa                                  >  1:1282980/1‑150 (MQ=255)
                                                                                                                           ttgATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATcc                            >  1:1866821/1‑150 (MQ=255)
                                                                                                                            tgATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCg                           <  1:576252/150‑1 (MQ=255)
                                                                                                                             gATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCgg                          >  1:523347/1‑150 (MQ=255)
                                                                                                                              aggacGGATAAGGCGCTCGCGCCGCATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGc                         >  1:1463596/5‑150 (MQ=25)
                                                                                                                              aTGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGc                         <  1:1264908/150‑1 (MQ=255)
                                                                                                                              aTGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGCCGTTTGTCATCAGTCTCAGGCCGGGAGCGTAATACCATCCGGc                         >  2:1329889/1‑150 (MQ=255)
                                                                                                                                 ccGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATTGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCAt                      >  2:422406/1‑150 (MQ=255)
                                                                                                                                    gataaggcgCTCGCGCCGCATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCAtttt                   >  1:1971483/6‑150 (MQ=37)
                                                                                                                                        cggcgCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAtt               <  2:1939587/150‑1 (MQ=255)
                                                                                                                                        cggcgCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAtt               >  1:2019515/1‑150 (MQ=255)
                                                                                                                                         ggcgCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAttt              <  1:2090409/150‑1 (MQ=255)
                                                                                                                                          gcgcGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTa             >  2:1342187/1‑150 (MQ=255)
                                                                                                                                            gcgAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATATATTGCAATTCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAg           <  1:2390645/150‑1 (MQ=255)
                                                                                                                                             cgAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGc          >  1:1638463/1‑150 (MQ=255)
                                                                                                                                             cgAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGc          <  1:2365654/150‑1 (MQ=255)
                                                                                                                                             cgAGCGCATTATCCGTCCTACAAAACATCCCAAATTCAATAGATTGCAATGCACTTGTAGGCATCATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGc          >  2:1521151/1‑150 (MQ=255)
                                                                                                                                                gcgcCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCAc       >  1:2461050/1‑150 (MQ=255)
                                                                                                                                                gcgcCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCAc       >  1:1250916/1‑150 (MQ=255)
                                                                                                                                                 cgcCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACt      >  1:361854/1‑150 (MQ=255)
                                                                                                                                                 cgcCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACt      >  1:926540/1‑150 (MQ=255)
                                                                                                                                                     ttATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCa  <  1:1159297/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCA  >  NC_000913/3012329‑3012627

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: