Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,012,472 A→C 5.1% intergenic (+79/+142) 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,4720AC 5.1% 594.8 / 20.0 196intergenic (+79/+142)yqeA/yqeBputative amino acid kinase YqeA/XdhC‑CoxI family protein YqeB
Reads supporting (aligned to +/- strand):  ref base A (95/91);  new base C (5/5);  total (100/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

Only 100 of 199 total aligned reads displayed.

CACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGC  >  NC_000913/3012323‑3012619
                                                                                                                                                     |                                                                                                                                                   
cACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTAAAGTAAGTGCGCTTTGTTGATGCCGGATGAGGCGCGa                                                                                                                                                     >  2:1285725/1‑150 (MQ=255)
  cTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGGAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAgc                                                                                                                                                   <  2:560279/150‑1 (MQ=255)
    gATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAgcgc                                                                                                                                                 >  1:1785479/1‑150 (MQ=255)
    gATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCAAAAGATTTTAAAGTAAGTGCGCTATGTTGATTAATGATGCGGCGCGAgcgc                                                                                                                                                 >  2:1264908/1‑150 (MQ=255)
     aTTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGGGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGcc                                                                                                                                                <  2:1307166/150‑1 (MQ=255)
     aTTACCGTAAAGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGcc                                                                                                                                                >  2:1917708/1‑150 (MQ=255)
      ttACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCt                                                                                                                                               >  1:218440/1‑150 (MQ=255)
      ttACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCt                                                                                                                                               <  1:2125614/150‑1 (MQ=255)
         ccGTATTGAGCAAAGCGAAAGAAAGGATTGAAGGCAAAACAGGAACGGTGATTTTTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGAGCTTTGTTGATTCCGGATGAGGCGCGACCGCCTTAt                                                                                                                                            >  2:1049899/1‑150 (MQ=255)
          cGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                                           <  1:43620/150‑1 (MQ=255)
             aTTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGt                                                                                                                                        <  1:919401/150‑1 (MQ=255)
              ttGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATAAGGCGCTCGCGCCGCATCCggc                                                                                                                                       <  1:2080621/150‑3 (MQ=255)
               tGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGTTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATAAGGCGCTCGCGCCGCATCCggca                                                                                                                                      <  1:196866/150‑4 (MQ=255)
               tGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTcc                                                                                                                                      >  1:1056045/1‑150 (MQ=255)
               tGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTcc                                                                                                                                      <  1:1800186/150‑1 (MQ=255)
                gAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCt                                                                                                                                     <  1:1823257/150‑1 (MQ=255)
                gAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATAAGGCGCTCGCGCCGCATCCggcat                                                                                                                                     <  1:1435077/150‑5 (MQ=255)
                 aGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTa                                                                                                                                    >  2:703293/1‑150 (MQ=255)
                  gCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTAc                                                                                                                                   <  1:1136282/150‑1 (MQ=255)
                  gCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTAc                                                                                                                                   <  1:1979160/150‑1 (MQ=255)
                  gCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTAc                                                                                                                                   >  1:1956132/1‑150 (MQ=255)
                   cAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGTGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACa                                                                                                                                  <  2:1253795/150‑1 (MQ=255)
                    aaaGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACaa                                                                                                                                 >  2:1755024/1‑150 (MQ=255)
                         gAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACAt                                                                                                                            <  1:779204/150‑1 (MQ=255)
                         gAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACAt                                                                                                                            >  2:1131805/1‑150 (MQ=255)
                         gAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACAt                                                                                                                            >  1:488663/1‑150 (MQ=255)
                           aagaagGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCg                                                                                                                          >  1:1410535/1‑150 (MQ=255)
                             gaagGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCa                                                                                                                        >  1:543090/1‑150 (MQ=255)
                             gaagGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCa                                                                                                                        >  2:1502023/1‑150 (MQ=255)
                                  gATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTc                                                                                                                   <  2:1155104/150‑1 (MQ=255)
                                   aTTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCa                                                                                                                  <  2:2036806/150‑1 (MQ=255)
                                      gAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGAGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATa                                                                                                               >  2:853883/1‑150 (MQ=255)
                                          gCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGAtt                                                                                                           >  1:583427/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)
                                                cAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATg                                                                                                     >  2:1839236/1‑150 (MQ=255)
                                                  ggAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCa                                                                                                   <  2:1406338/150‑1 (MQ=255)
                                                  ggAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCa                                                                                                   >  1:1043450/1‑150 (MQ=255)
                                                      cGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGGGCTTg                                                                                               >  1:1102854/1‑150 (MQ=255)
                                                      cGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTg                                                                                               <  2:541949/150‑1 (MQ=255)
                                                            tttGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGcc                                                                                         <  1:294615/150‑1 (MQ=255)
                                                              tGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTg                                                                                       <  2:2324324/150‑1 (MQ=255)
                                                              tGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTg                                                                                       >  2:520704/1‑150 (MQ=255)
                                                               gTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGa                                                                                      >  1:687535/1‑150 (MQ=255)
                                                                tCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGAt                                                                                     <  1:1160564/150‑1 (MQ=255)
                                                                 cAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATa                                                                                    >  1:1321404/1‑150 (MQ=255)
                                                                 cAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATa                                                                                    <  2:1056045/150‑1 (MQ=255)
                                                                   gTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAg                                                                                  >  1:1041230/1‑150 (MQ=255)
                                                                   gTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAg                                                                                  >  2:977624/1‑150 (MQ=255)
                                                                   gTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAg                                                                                  >  1:1352969/1‑150 (MQ=255)
                                                                     aaTTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCg                                                                                <  2:1429270/150‑1 (MQ=255)
                                                                      aTTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGt                                                                               >  1:1332651/1‑150 (MQ=255)
                                                                        ttAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAg                                                                             <  1:1687412/150‑1 (MQ=255)
                                                                             ggCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCAt                                                                        >  2:402712/1‑150 (MQ=255)
                                                                               ccGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCa                                                                      <  2:2299795/150‑1 (MQ=255)
                                                                                cGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                     >  2:1398155/1‑150 (MQ=255)
                                                                                  ggtgggtAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                   <  1:402712/150‑1 (MQ=255)
                                                                                     gggtAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAAt                                                                <  1:2006388/150‑1 (MQ=255)
                                                                                         aaTACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTgg                                                            <  1:1755024/150‑1 (MQ=255)
                                                                                         aaTACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTgg                                                            <  1:323536/150‑1 (MQ=255)
                                                                                             ccATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGtt                                                        >  2:1546001/1‑150 (MQ=255)
                                                                                               aTCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTg                                                      <  1:2297548/150‑1 (MQ=255)
                                                                                                tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGt                                                     <  2:1332641/150‑1 (MQ=255)
                                                                                                   ggCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGtcat                                                  >  2:1163783/1‑150 (MQ=255)
                                                                                                    gCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGtcatc                                                 >  2:800460/1‑150 (MQ=255)
                                                                                                      cAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAg                                               <  2:1457100/150‑1 (MQ=255)
                                                                                                      cAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAg                                               <  2:1603354/150‑1 (MQ=255)
                                                                                                      cAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAg                                               >  2:1645552/1‑150 (MQ=255)
                                                                                                      cAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAg                                               >  2:353250/1‑150 (MQ=255)
                                                                                                          gATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGtctc                                           <  1:1141301/150‑1 (MQ=255)
                                                                                                          gATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGtctc                                           <  2:231018/150‑1 (MQ=255)
                                                                                                             tGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTAATCAGTCTCAgg                                        <  2:1706792/150‑1 (MQ=255)
                                                                                                              gTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGc                                       >  2:172462/1‑150 (MQ=255)
                                                                                                              gTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGc                                       >  2:194272/1‑150 (MQ=255)
                                                                                                              gTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGc                                       >  2:879009/1‑150 (MQ=255)
                                                                                                              gTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGc                                       <  2:2278298/150‑1 (MQ=255)
                                                                                                              gTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGAGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGc                                       >  2:1870597/1‑150 (MQ=255)
                                                                                                                cAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCg                                     >  2:2442001/1‑150 (MQ=255)
                                                                                                                cAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCg                                     >  1:2318708/1‑150 (MQ=255)
                                                                                                                     aagtGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtgg                                <  1:2062303/150‑1 (MQ=255)
                                                                                                                      agtGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtggg                               <  1:963334/150‑1 (MQ=255)
                                                                                                                       gggcgcTTTGTCGACGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtgggt                              <  2:1011716/148‑1 (MQ=255)
                                                                                                                           gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATa                          >  1:1282980/1‑150 (MQ=255)
                                                                                                                                  tgATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCg                   <  1:576252/150‑1 (MQ=255)
                                                                                                                                   gATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCgg                  >  1:523347/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)
                                                                                                                                       acGGATAAGGCGCTCGCGCCGCATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCAt              >  2:2249744/2‑150 (MQ=37)
                                                                                                                                          gataaggcgCTCGCGCCGCATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCAtttt           >  2:2387883/6‑150 (MQ=37)
                                                                                                                                            tGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCa         >  2:2063608/1‑150 (MQ=255)
                                                                                                                                             gcggcgCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAt        <  1:1917708/150‑1 (MQ=255)
                                                                                                                                              cggcgCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAtt       >  2:2367673/1‑150 (MQ=255)
                                                                                                                                               ggcgCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAttt      >  2:1220432/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)
                                                                                                                                                   cgAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGc  <  1:1881676/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGC  >  NC_000913/3012323‑3012619

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: