Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 192,713 T→C 6.5% intergenic (+133/‑159) pyrH → / → frr uridylate kinase/ribosome recycling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913192,7130TC 6.5% 624.9 / 22.9 186intergenic (+133/‑159)pyrH/frruridylate kinase/ribosome recycling factor
Reads supporting (aligned to +/- strand):  ref base T (96/77);  new base C (6/6);  total (103/83)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 187 total aligned reads displayed.

CTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTCTGACTGAGACAAGTTTTCAAG  >  NC_000913/192564‑192861
                                                                                                                                                     |                                                                                                                                                    
cTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGt                                                                                                                                                      <  2:313520/150‑1 (MQ=255)
  ttAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGa                                                                                                                                                    <  1:610165/150‑1 (MQ=255)
  ttAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGa                                                                                                                                                    <  1:2128882/150‑1 (MQ=255)
  ttAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGa                                                                                                                                                    <  1:2097692/150‑1 (MQ=255)
  ttAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGa                                                                                                                                                    <  2:1361315/150‑1 (MQ=255)
   tAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAg                                                                                                                                                   <  1:415722/150‑1 (MQ=255)
   tAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAg                                                                                                                                                   <  2:1836593/150‑1 (MQ=255)
    aaTCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAgg                                                                                                                                                  <  2:1346152/150‑1 (MQ=255)
    aaTCACGGAATAAGGCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAgg                                                                                                                                                  <  1:510388/150‑1 (MQ=255)
      tCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGAt                                                                                                                                                >  1:378217/1‑150 (MQ=255)
           gAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGa                                                                                                                                           >  1:1590183/1‑150 (MQ=255)
              taatTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCa                                                                                                                                        >  2:1149622/1‑150 (MQ=255)
                atTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACt                                                                                                                                      >  1:1468869/1‑150 (MQ=255)
                 ttCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTg                                                                                                                                     >  2:2498392/1‑150 (MQ=255)
                    ccGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGccc                                                                                                                                  <  1:1148691/150‑1 (MQ=255)
                    ccGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGCGAGGATGGCGAGCACTGccc                                                                                                                                  >  2:477944/1‑150 (MQ=255)
                             ataaataaGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCaaa                                                                                                                         >  1:76049/1‑150 (MQ=255)
                                aataaGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATg                                                                                                                      <  2:1262482/150‑1 (MQ=255)
                                     gggTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCaaa                                                                                                                 <  1:650743/150‑1 (MQ=255)
                                       gTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCaaata                                                                                                               >  2:1949381/1‑150 (MQ=255)
                                        tAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCaaataa                                                                                                              >  1:933739/1‑150 (MQ=255)
                                        tAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCaaataa                                                                                                              >  2:252323/1‑150 (MQ=255)
                                        tAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCaaataa                                                                                                              >  1:38352/1‑150 (MQ=255)
                                         aaGATTCCGCTTAAGTATCGGGGGGGCGTAAGTCTGGTTATATGGCGTTATTGTTGAAGGCATTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGTATGGCGAGCACTGCCCGGGTCAAAAATGGCaaataaa                                                                                                             <  2:1847166/150‑1 (MQ=255)
                                         aaGATTCCGCGTACGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACACAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCaaataaa                                                                                                             <  2:1200515/150‑1 (MQ=255)
                                             ttCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAg                                                                                                         >  1:31065/1‑150 (MQ=255)
                                             ttCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAg                                                                                                         >  1:1295380/1‑150 (MQ=255)
                                                         aTCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCa                                                                                             <  2:95953/150‑1 (MQ=255)
                                                          tCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAg                                                                                            >  1:2179129/1‑150 (MQ=255)
                                                          tCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAg                                                                                            >  2:1136169/1‑150 (MQ=255)
                                                          tCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAg                                                                                            >  1:2179240/1‑150 (MQ=255)
                                                          tCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAg                                                                                            >  1:2038657/1‑150 (MQ=255)
                                                                  gCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTAc                                                                                    <  1:46476/150‑1 (MQ=255)
                                                                  gCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTAc                                                                                    <  1:1646606/150‑1 (MQ=255)
                                                                   cGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTAcc                                                                                   <  2:378217/150‑1 (MQ=255)
                                                                    gTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTAccc                                                                                  >  1:2293216/1‑150 (MQ=255)
                                                                      aaGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGt                                                                                <  2:2127864/150‑1 (MQ=255)
                                                                       aGTCTGGTTATCAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGGACAAAAGTCCGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTa                                                                               <  2:2585004/150‑1 (MQ=255)
                                                                        gTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTaa                                                                              <  2:719402/150‑1 (MQ=255)
                                                                         tCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAat                                                                             >  2:2144810/1‑150 (MQ=255)
                                                                           tGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAatat                                                                           >  2:2360667/1‑150 (MQ=255)
                                                                           tGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAatat                                                                           >  2:2404285/1‑150 (MQ=255)
                                                                             gTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGt                                                                         >  1:2120208/1‑150 (MQ=255)
                                                                             gTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGt                                                                         >  1:588796/1‑150 (MQ=255)
                                                                              ttATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGtt                                                                        <  1:2263110/150‑1 (MQ=255)
                                                                              ttATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGtt                                                                        <  2:1059239/150‑1 (MQ=255)
                                                                              ttATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGtt                                                                        <  2:2510922/150‑1 (MQ=255)
                                                                                ataAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTa                                                                      >  2:2293673/1‑150 (MQ=255)
                                                                                ataAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTa                                                                      >  2:1952952/1‑150 (MQ=255)
                                                                                ataAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTa                                                                      >  1:6840/1‑150 (MQ=255)
                                                                                  aaGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAAt                                                                    >  1:2558417/1‑150 (MQ=255)
                                                                                  aaGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAAt                                                                    >  1:1103795/1‑150 (MQ=255)
                                                                                    ggCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAAAGGCAAATAAAATATCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCa                                                                  >  1:248650/1‑150 (MQ=255)
                                                                                     gCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAg                                                                 >  1:735488/1‑150 (MQ=255)
                                                                                     gCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAg                                                                 >  2:2469813/1‑150 (MQ=255)
                                                                                     gCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAg                                                                 >  2:1124160/1‑150 (MQ=255)
                                                                                     gCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAg                                                                 >  1:2186460/1‑150 (MQ=255)
                                                                                       gTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAggg                                                               >  1:808833/1‑150 (MQ=255)
                                                                                       gTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAggg                                                               <  2:1461500/150‑1 (MQ=255)
                                                                                                 caggcagTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAg                                                     <  2:862283/150‑1 (MQ=255)
                                                                                                 caggcagTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAg                                                     >  1:1709505/1‑150 (MQ=255)
                                                                                                   ggcagTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGTCGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGca                                                   >  2:2035639/1‑150 (MQ=255)
                                                                                                   ggcagTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGTCGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGca                                                   >  2:1554875/1‑150 (MQ=255)
                                                                                                    gcagTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGcac                                                  >  1:2399790/1‑150 (MQ=255)
                                                                                                     cagTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGcaca                                                 <  1:1616697/150‑1 (MQ=255)
                                                                                                        tttGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACtt                                              >  1:260579/1‑150 (MQ=255)
                                                                                                         ttGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGCGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTc                                             <  1:477944/150‑1 (MQ=255)
                                                                                                         ttGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGCGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTc                                             <  2:57603/150‑1 (MQ=255)
                                                                                                         ttGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGCGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTc                                             <  1:5614/150‑1 (MQ=255)
                                                                                                             tCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACt                                         >  2:991258/1‑150 (MQ=255)
                                                                                                             tCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACt                                         <  1:1026957/150‑1 (MQ=255)
                                                                                                             tCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACt                                         >  2:136631/1‑150 (MQ=255)
                                                                                                                cGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACtgtg                                      <  1:2497785/150‑1 (MQ=255)
                                                                                                                       cgcgcAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACt                               <  2:1683297/150‑1 (MQ=255)
                                                                                                                       cgcgcAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACt                               >  2:646984/1‑150 (MQ=255)
                                                                                                                           cagcaACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTgtct                           >  1:108024/1‑150 (MQ=255)
                                                                                                                            agcaACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTgtctg                          <  1:971844/150‑1 (MQ=255)
                                                                                                                            agcaACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTgtctg                          >  1:897582/1‑150 (MQ=255)
                                                                                                                            agcaACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTgtctg                          <  1:477672/150‑1 (MQ=255)
                                                                                                                            agcaACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTgtctg                          >  1:1359655/1‑150 (MQ=255)
                                                                                                                            agcaACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTgtctg                          >  1:1685312/1‑150 (MQ=255)
                                                                                                                               aaCCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTgtctggtc                       >  2:1719183/1‑150 (MQ=255)
                                                                                                                                cccGGAGCGTACACACGGACGTGAGGCTGGCGAGCCCGGCCCGGGGCCAAAATGGCAAATACACTAGCCTAATAATCCAGACGAGTACCCGTAATATGTTTCATCAGGGCTCTCCTTAGCACACTTCCACTGTGTGTGACTgtctggtct                      <  2:2545929/149‑1 (MQ=255)
                                                                                                                                 ccGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTgtctggtctg                     >  2:624723/1‑150 (MQ=255)
                                                                                                                                  cGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTctga                    <  1:2432593/150‑1 (MQ=255)
                                                                                                                                  cGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTctga                    <  1:2530642/150‑1 (MQ=255)
                                                                                                                                      gCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTctgactga                <  1:2569155/150‑1 (MQ=255)
                                                                                                                                       cGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTCTGACTgag               >  1:2303881/1‑150 (MQ=255)
                                                                                                                                          aCAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTCTGACTGAGACa            >  1:1597079/1‑150 (MQ=255)
                                                                                                                                           cAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTCTGACTGAGACaa           <  1:1120079/150‑1 (MQ=255)
                                                                                                                                           cAAAAGGACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTCTGACTGAGACaa           <  2:2089250/150‑1 (MQ=255)
                                                                                                                                            aaaaGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTCTGACTGAGACAAg          >  1:1880955/1‑150 (MQ=255)
                                                                                                                                             aaaGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTCTGACTGAGACAAGt         <  1:1553115/150‑1 (MQ=255)
                                                                                                                                             aaaGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTCTGACTGAGACAAGt         <  2:792796/150‑1 (MQ=255)
                                                                                                                                                gTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTCTGACTGAGACAAGtttt      >  1:1201910/1‑150 (MQ=255)
                                                                                                                                                gTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTCTGACTGAGACAAGtttt      >  1:1175975/1‑150 (MQ=255)
                                                                                                                                                   cGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTCTGACTGAGACAAGTTTTCaa   <  2:943783/150‑1 (MQ=255)
                                                                                                                                                   cGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTCTGACTGAGACAAGTTTTCaa   <  2:1880955/150‑1 (MQ=255)
                                                                                                                                                    gTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGTTTACCCGTAATATGTTTTATCAGGGATATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTCTGACTGAGACAAGTTTTCAAg  <  1:1352498/150‑1 (MQ=255)
                                                                                                                                                    gTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTCTGACTGAGACAAGTTTTCAAg  <  2:664393/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTTTAATCACGGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGAGCACTGCCCGGGGCCAAAATGGCAAATAAAATAGCCTAATAATCCAGACGATTACCCGTAATATGTTTAATCAGGGCTATACTTAGCACACTTCCACTGTGTGTGACTGTCTGGTCTGACTGAGACAAGTTTTCAAG  >  NC_000913/192564‑192861

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