Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,012,479 T→C 5.2% intergenic (+86/+135) 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,4790TC 5.2% 606.6 / 14.9 193intergenic (+86/+135)yqeA/yqeBputative amino acid kinase YqeA/XdhC‑CoxI family protein YqeB
Reads supporting (aligned to +/- strand):  ref base T (94/89);  new base C (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 = 8.15e-01

Only 100 of 197 total aligned reads displayed.

CCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAG  >  NC_000913/3012332‑3012628
                                                                                                                                                   |                                                                                                                                                     
 cGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                                                    <  1:43620/150‑1 (MQ=255)
   tATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCTCGCGCCGCATCCg                                                                                                                                                  <  2:1055181/150‑1 (MQ=255)
    aTTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGt                                                                                                                                                 <  1:919401/150‑1 (MQ=255)
      tGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTcc                                                                                                                                               >  1:1056045/1‑150 (MQ=255)
       gAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCt                                                                                                                                              >  2:1911837/1‑150 (MQ=255)
       gAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATAAGGCGCTCGCGCCGCATCCggcat                                                                                                                                              <  1:1435077/150‑5 (MQ=255)
        aGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTa                                                                                                                                             >  2:703293/1‑150 (MQ=255)
         gCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTAc                                                                                                                                            >  1:1939587/1‑150 (MQ=255)
         gCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTAc                                                                                                                                            <  1:1136282/150‑1 (MQ=255)
          cAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACa                                                                                                                                           >  2:172161/1‑150 (MQ=255)
                gAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACAt                                                                                                                                     >  1:488663/1‑150 (MQ=255)
                gAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACAt                                                                                                                                     <  1:779204/150‑1 (MQ=255)
                  aagaagGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCg                                                                                                                                   >  1:1410535/1‑150 (MQ=255)
                    gaagGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCa                                                                                                                                 >  2:1502023/1‑150 (MQ=255)
                         gATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTc                                                                                                                            <  2:1155104/150‑1 (MQ=255)
                            tGAAGGCAGAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAAt                                                                                                                         <  2:529807/150‑1 (MQ=255)
                             gAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGAGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATa                                                                                                                        >  2:853883/1‑150 (MQ=255)
                                 gCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGAtt                                                                                                                    >  1:583427/1‑150 (MQ=255)
                                  cAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTg                                                                                                                   >  2:1667492/1‑150 (MQ=255)
                                  cAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTg                                                                                                                   >  2:921166/1‑150 (MQ=255)
                                     aaCAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCaa                                                                                                                <  1:1450778/150‑1 (MQ=255)
                                     aaCAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCaa                                                                                                                >  2:606502/1‑150 (MQ=255)
                                        aGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGc                                                                                                             <  2:1104581/150‑1 (MQ=255)
                                        aGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGc                                                                                                             <  1:1963716/150‑1 (MQ=255)
                                         ggAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCa                                                                                                            >  1:1043450/1‑150 (MQ=255)
                                         ggAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCa                                                                                                            <  2:1406338/150‑1 (MQ=255)
                                         ggAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCa                                                                                                            <  2:186005/150‑1 (MQ=255)
                                         ggAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCa                                                                                                            <  1:1999542/150‑1 (MQ=255)
                                             cGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGGGCTTg                                                                                                        >  1:1102854/1‑150 (MQ=255)
                                                   tttGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGcc                                                                                                  <  1:294615/150‑1 (MQ=255)
                                                   tttGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGcc                                                                                                  >  1:1803240/1‑150 (MQ=255)
                                                   tttGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGcc                                                                                                  >  1:1794232/1‑150 (MQ=255)
                                                     tGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTg                                                                                                <  2:2324324/150‑1 (MQ=255)
                                                      gTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGa                                                                                               >  1:687535/1‑150 (MQ=255)
                                                        cAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATa                                                                                             >  1:1321404/1‑150 (MQ=255)
                                                          gTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGGACTTGTAGGCCTGATAAg                                                                                           >  2:1734097/1‑150 (MQ=255)
                                                          gTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAg                                                                                           >  1:1352969/1‑150 (MQ=255)
                                                          gTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAg                                                                                           >  2:977624/1‑150 (MQ=255)
                                                          gTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAg                                                                                           >  1:1041230/1‑150 (MQ=255)
                                                          gTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAg                                                                                           <  2:1501948/150‑1 (MQ=255)
                                                           tAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGc                                                                                          <  2:543090/150‑1 (MQ=255)
                                                             aTTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGt                                                                                        >  1:1332651/1‑150 (MQ=255)
                                                             aTTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGt                                                                                        <  1:230258/150‑1 (MQ=255)
                                                               ttAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAg                                                                                      <  1:1687412/150‑1 (MQ=255)
                                                                  aaGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAgcgc                                                                                   >  1:538713/1‑150 (MQ=255)
                                                                    ggCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCAt                                                                                 >  2:402712/1‑150 (MQ=255)
                                                                    ggCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCAt                                                                                 >  1:235110/1‑150 (MQ=255)
                                                                      ccGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCa                                                                               <  2:2299795/150‑1 (MQ=255)
                                                                       cGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                              >  2:1398155/1‑150 (MQ=255)
                                                                       cGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                              <  1:942980/150‑1 (MQ=255)
                                                                                aaTACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTgg                                                                     <  1:1755024/150‑1 (MQ=255)
                                                                                aaTACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTgg                                                                     <  1:323536/150‑1 (MQ=255)
                                                                                aaTACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTgg                                                                     <  1:38766/150‑1 (MQ=255)
                                                                                   aCCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGt                                                                  <  2:1043450/150‑1 (MQ=255)
                                                                                    ccATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGtt                                                                 >  2:1546001/1‑150 (MQ=255)
                                                                                     cATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGttt                                                                <  1:977624/150‑1 (MQ=255)
                                                                                       tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGt                                                              <  2:1332641/150‑1 (MQ=255)
                                                                                         cGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGtca                                                            >  1:1900654/1‑150 (MQ=255)
                                                                                          ggCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGtcat                                                           >  2:1163783/1‑150 (MQ=255)
                                                                                           gCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGtcatc                                                          >  2:800460/1‑150 (MQ=255)
                                                                                             cAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAg                                                        <  2:1603354/150‑1 (MQ=255)
                                                                                             cAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAg                                                        >  2:353250/1‑150 (MQ=255)
                                                                                              aaaGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGt                                                       <  1:2071133/150‑1 (MQ=255)
                                                                                              aaaGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGt                                                       >  2:1264701/1‑150 (MQ=255)
                                                                                                 gATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGtctc                                                    <  1:1141301/150‑1 (MQ=255)
                                                                                                 gATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGtctc                                                    <  2:231018/150‑1 (MQ=255)
                                                                                                   ttGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAg                                                  >  2:504835/1‑150 (MQ=255)
                                                                                                     gTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGc                                                >  2:194272/1‑150 (MQ=255)
                                                                                                     gTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGc                                                <  2:2278298/150‑1 (MQ=255)
                                                                                                     gTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGc                                                >  2:879009/1‑150 (MQ=255)
                                                                                                       cAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCg                                              >  2:2442001/1‑150 (MQ=255)
                                                                                                        aagtaagtGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                             >  1:2256987/1‑150 (MQ=255)
                                                                                                              gtGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtgggt                                       <  2:1785479/150‑1 (MQ=255)
                                                                                                                gcgcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTaa                                     >  2:1450469/1‑150 (MQ=255)
                                                                                                                  gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATa                                   >  1:1282980/1‑150 (MQ=255)
                                                                                                                        ttgATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATcc                             >  1:1866821/1‑150 (MQ=255)
                                                                                                                          gATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCgg                           >  1:575427/1‑150 (MQ=255)
                                                                                                                          gATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCgg                           <  1:2462039/150‑1 (MQ=255)
                                                                                                                           aggacGGATAAGGCGCTCGCGCCGCATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGc                          >  1:1463596/5‑150 (MQ=25)
                                                                                                                           aTGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGc                          <  1:1264908/150‑1 (MQ=255)
                                                                                                                           aTGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGc                          >  2:456298/1‑150 (MQ=255)
                                                                                                                             gCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCAATTGTAGGCCTGATAAGCGTAGCGGATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCa                        >  2:1287704/1‑150 (MQ=255)
                                                                                                                              acGGATAAGGCGCTCGCGCCGCATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCAt                       >  2:2249744/2‑150 (MQ=37)
                                                                                                                                 gataaggcgCTCGCGCCGCATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCAtttt                    >  2:2387883/6‑150 (MQ=37)
                                                                                                                                 gataaggcgCTCGCGCCGCATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCAtttt                    >  1:1971483/6‑150 (MQ=37)
                                                                                                                                   tGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCa                  >  2:2063608/1‑150 (MQ=255)
                                                                                                                                     cggcgCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAtt                >  1:2019515/1‑150 (MQ=255)
                                                                                                                                     cggcgCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAtt                >  2:2367673/1‑150 (MQ=255)
                                                                                                                                      ggcgCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAttt               >  2:1220432/1‑150 (MQ=255)
                                                                                                                                      ggcgCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAttt               <  1:2090409/150‑1 (MQ=255)
                                                                                                                                         gcgAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATATATTGCAATTCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAg            <  1:2390645/150‑1 (MQ=255)
                                                                                                                                          cgAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGc           <  1:2365654/150‑1 (MQ=255)
                                                                                                                                          cgAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGc           >  1:1638463/1‑150 (MQ=255)
                                                                                                                                             gcgcCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCAc        >  1:405502/1‑150 (MQ=255)
                                                                                                                                             gcgcCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCAc        >  1:1250916/1‑150 (MQ=255)
                                                                                                                                                ccTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTc     <  1:520704/150‑1 (MQ=255)
                                                                                                                                                 cTTATCCGTCCTACAAAACCTCGCAAATTCAATAGATTTCAATTCACTTGTTGGCCTGATAAGCTTATAGAATCATTCAATTTGTCTTTTTTCATCAGTCTCATGCCTTGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTcc    <  2:1251960/150‑1 (MQ=255)
                                                                                                                                                  ttATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCa   <  2:1638463/150‑1 (MQ=255)
                                                                                                                                                  ttATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCa   <  1:1159297/150‑1 (MQ=255)
                                                                                                                                                   tATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAg  >  1:2440372/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAG  >  NC_000913/3012332‑3012628

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: