Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,636,008 G→T 7.8% intergenic (+115/+200) yhiM → / ← yhiN inner membrane protein with a role in acid resistance/putative oxidoreductase YhiN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,636,0080GT 7.8% 335.8 / 11.7 102intergenic (+115/+200)yhiM/yhiNinner membrane protein with a role in acid resistance/putative oxidoreductase YhiN
Reads supporting (aligned to +/- strand):  ref base G (44/50);  new base T (6/2);  total (50/52)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

Only 100 of 105 total aligned reads displayed.

TTTCAATATTAGAAGCGGGTTCTGCTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAA  >  NC_000913/3635859‑3636157
                                                                                                                                                     |                                                                                                                                                     
tttCAATATTAGAAGCGGGTTCTGCTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAg                                                                                                                                                       >  2:1384777/1‑150 (MQ=255)
    aaTATTAGAAGCGGGTTCTGCTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTAc                                                                                                                                                   >  1:180664/1‑150 (MQ=255)
     atatTAGAAGCTGGTTCTGCTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACt                                                                                                                                                  <  1:1152876/150‑1 (MQ=255)
      tatTAGAAGCGGGTTCTGCTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTg                                                                                                                                                 >  1:1020839/1‑150 (MQ=255)
      tatTAGAAGCGGGTTCTGCTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTg                                                                                                                                                 <  1:1011977/150‑1 (MQ=255)
       atTAGAAGCGGGTTCTGCTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGAGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGAATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGt                                                                                                                                                >  2:1111318/1‑150 (MQ=255)
          aGAAGCGGGTTCTGCTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGTAACTTTGCTGtttt                                                                                                                                             >  1:966404/1‑150 (MQ=255)
          aGAAGCGGGTTCTGCTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGtttt                                                                                                                                             >  1:333415/1‑150 (MQ=255)
           gAAGCGGGTTCTGCTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGttttt                                                                                                                                            <  2:1488707/150‑1 (MQ=255)
            aaGCGGGTTCTGCTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGtttttt                                                                                                                                           <  1:1322706/150‑1 (MQ=255)
              gCGGGTTCTGCTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGTAACTTTGCTGttttttat                                                                                                                                         <  1:1161868/150‑3 (MQ=255)
                gggTTCTGCTATAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTATGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCaa                                                                                                                                       <  1:333203/150‑1 (MQ=255)
                gggTTCTGCTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCaa                                                                                                                                       <  2:1560289/150‑1 (MQ=255)
                  gTTCTGCTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAAcc                                                                                                                                     <  2:474415/150‑1 (MQ=255)
                   ttCTGCTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGTAACTTTGCTGttttttatactta                                                                                                                                    <  2:1412983/150‑8 (MQ=255)
                        cTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCa                                                                                                                               <  2:70311/150‑1 (MQ=255)
                         tAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGCCCATTTGCGTAACGCTGGTTTTTCTTAGGCCTCATGAAATACCGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACGGTGGTTTTCAACCTGTTCat                                                                                                                              >  2:1267958/1‑150 (MQ=255)
                          aaaaaaTAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCata                                                                                                                             >  2:315300/1‑150 (MQ=255)
                          aaaaaaTAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCata                                                                                                                             >  2:1034760/1‑150 (MQ=255)
                          aaaaaaTAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCata                                                                                                                             >  2:967880/1‑150 (MQ=255)
                           aaaaaTAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCatat                                                                                                                            <  1:983224/150‑1 (MQ=255)
                              aataatTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTc                                                                                                                         <  1:1034760/150‑1 (MQ=255)
                               ataatTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCa                                                                                                                        <  1:86520/150‑1 (MQ=255)
                               ataatTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCa                                                                                                                        <  2:829431/150‑1 (MQ=255)
                               ataatTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCa                                                                                                                        <  2:1178484/150‑1 (MQ=255)
                                  atTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATaa                                                                                                                     <  1:103222/150‑1 (MQ=255)
                                     gCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGa                                                                                                                  >  2:1302206/1‑150 (MQ=255)
                                     gCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGa                                                                                                                  >  2:1365541/1‑150 (MQ=255)
                                     gCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGa                                                                                                                  >  2:219863/1‑150 (MQ=255)
                                      cAACGTACCGGATAAAACCAGCGTTGCCCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACAGTTTTTTTCAACCTGTTCATATTTCATAAAGAt                                                                                                                 >  2:1157035/1‑150 (MQ=255)
                                      cAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGAt                                                                                                                 >  1:1612912/1‑150 (MQ=255)
                                      cAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGAt                                                                                                                 >  1:1058962/1‑150 (MQ=255)
                                          gTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTgg                                                                                                             <  2:333415/150‑1 (MQ=255)
                                            aCCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGcc                                                                                                           >  2:490878/1‑150 (MQ=255)
                                             ccGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCa                                                                                                          <  1:950500/150‑1 (MQ=255)
                                                gATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACtg                                                                                                       <  1:666290/150‑1 (MQ=255)
                                                gATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACtg                                                                                                       <  1:302958/150‑1 (MQ=255)
                                                     aaCCAGCGTTGACGATTTGCGTAACGCTGGTTTTTCTTAGGAATCATCAAATAACGCAGCATTAATGCATAGTGGTTAAGTATTAAAAAAAGCAAAGTTAGTGTTTTTTTCAACCTGTTCATATTTCATAATGATCTGGCCAAtgagtga                                                                                                  >  2:817524/1‑150 (MQ=255)
                                                          gCGTTGACCATTTGCGTAACGCTGGTTTTTCTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTc                                                                                            >  1:1248205/1‑150 (MQ=255)
                                                            gTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGGATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCa                                                                                           >  1:759886/1‑150 (MQ=255)
                                                            gTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCa                                                                                           >  2:237977/1‑150 (MQ=255)
                                                              tGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAAt                                                                                         <  1:1509334/150‑1 (MQ=255)
                                                               gACCCTTTGCGTAACGCTGGTTTTTGTTATGCATCATTAAATAAAGCATCATTAATTCATAGTGGTTAACTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATc                                                                                        <  1:170924/150‑1 (MQ=255)
                                                                aCCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCa                                                                                       <  1:318334/150‑1 (MQ=255)
                                                                  cATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCAcc                                                                                     >  2:877761/1‑150 (MQ=255)
                                                                        cGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTtcca                                                                               <  1:315300/150‑1 (MQ=255)
                                                                        cGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGGTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATCTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTtcca                                                                               <  2:850100/150‑1 (MQ=255)
                                                                         gTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTtccat                                                                              >  1:1377245/1‑150 (MQ=255)
                                                                         gTAACGCTGGTTTTTCTTAGGCATCATGAAAAAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTtccat                                                                              <  1:1464355/150‑1 (MQ=255)
                                                                                  gTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTtat                                                                     <  2:1020839/150‑1 (MQ=255)
                                                                                   tttttCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTtata                                                                    >  2:1361527/1‑150 (MQ=255)
                                                                                      ttCTTATGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTa                                                                 <  1:1274009/150‑1 (MQ=255)
                                                                                      ttCTTAGGCATCATGAAATAACTCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTa                                                                 <  1:1194358/150‑1 (MQ=255)
                                                                                      ttCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTa                                                                 >  1:1539395/1‑150 (MQ=255)
                                                                                       tCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTaa                                                                <  1:1082046/150‑1 (MQ=255)
                                                                                        cTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAg                                                               >  2:855374/1‑150 (MQ=255)
                                                                                         ttAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGc                                                              >  1:940399/1‑150 (MQ=255)
                                                                                         ttAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGc                                                              >  1:992875/1‑150 (MQ=255)
                                                                                            ggCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATg                                                           >  1:143393/1‑150 (MQ=255)
                                                                                             gCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATgg                                                          >  2:299863/1‑150 (MQ=255)
                                                                                             gCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATgg                                                          >  2:1561992/1‑150 (MQ=255)
                                                                                               atcatGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAg                                                        <  1:1365541/150‑1 (MQ=255)
                                                                                               atcatGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAg                                                        <  1:1111318/150‑1 (MQ=255)
                                                                                                 catGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAggg                                                      >  2:1508939/1‑150 (MQ=255)
                                                                                                  atGAAATAACGCAGCATTAAAGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCGACCTTATATTAAGCATGGAGGGt                                                     >  2:543444/1‑150 (MQ=255)
                                                                                                    gAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGttt                                                   >  1:126890/1‑150 (MQ=255)
                                                                                                     aaaTAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTc                                                  <  1:855374/150‑1 (MQ=255)
                                                                                                       aTAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAg                                                <  1:992699/150‑1 (MQ=255)
                                                                                                       aTAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAg                                                <  1:265075/150‑1 (MQ=255)
                                                                                                       aTAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAg                                                <  1:745172/150‑1 (MQ=255)
                                                                                                       aTAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAg                                                <  1:219863/150‑1 (MQ=255)
                                                                                                        tAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGt                                               >  1:1602110/1‑150 (MQ=255)
                                                                                                        tAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGt                                               >  1:1619664/1‑150 (MQ=255)
                                                                                                                caTTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTtc                                                                                 >  1:478116/1‑108 (MQ=255)
                                                                                                                caTTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTtc                                                                                 <  2:478116/108‑1 (MQ=255)
                                                                                                                    aaTGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCAt                                   >  2:680156/1‑150 (MQ=255)
                                                                                                                     aTGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCAtt                                  >  2:125048/1‑150 (MQ=255)
                                                                                                                     aTGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCAtt                                  <  2:1612912/150‑1 (MQ=255)
                                                                                                                       gCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAg                                <  2:1177881/150‑1 (MQ=255)
                                                                                                                       gCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAg                                <  1:494090/150‑1 (MQ=255)
                                                                                                                         aTAGTGGTTTAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGaa                              <  2:905643/150‑1 (MQ=255)
                                                                                                                           aGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTGTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGaaaa                            <  2:331114/150‑1 (MQ=255)
                                                                                                                             tGGTTAAGTATAAAAAACAGCAAAGTTAATGTTTTTTTCAACCTGTACATATTTCAAAAAGATCTGGCGACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAaata                          <  1:817524/150‑1 (MQ=255)
                                                                                                                              ggTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAaataa                         >  1:55249/1‑150 (MQ=255)
                                                                                                                                ttAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCCGTTCACGGGCTCATTAGAAAATAATc                       <  1:237977/150‑1 (MQ=255)
                                                                                                                                 tAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATca                      <  2:143393/150‑1 (MQ=255)
                                                                                                                                 tAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATca                      >  1:468836/1‑150 (MQ=255)
                                                                                                                                   gttgaaaaaaaCAGTAACTTTGCTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATcaca                    >  1:1244181/6‑150 (MQ=255)
                                                                                                                                   aGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATaaca                    >  1:623681/1‑150 (MQ=255)
                                                                                                                                     tgaaaaaaaCAGTAACTTTGCTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAg                  >  2:1031621/4‑150 (MQ=255)
                                                                                                                                      gaaaaaaaCAGTAACTTTGCTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGa                 >  1:64141/3‑150 (MQ=255)
                                                                                                                                      ataAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGa                 >  1:111048/1‑150 (MQ=255)
                                                                                                                                       aaaaaaaCAGTAACTTTGCTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCAATAAACATTATATTAAGCATGGAGGGTTTCAGTTCACGGGATCATTAGAAAATAATCACAAGAt                >  1:1556716/2‑150 (MQ=255)
                                                                                                                                        aaaaaaCAGTAACTTTGCTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATa               >  1:1092884/1‑150 (MQ=255)
                                                                                                                                        aaaaaaCAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATa               <  1:648477/150‑1 (MQ=255)
                                                                                                                                             aCAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACCAGATAACCAc          <  1:456716/150‑1 (MQ=255)
                                                                                                                                               aGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTa        <  2:759886/150‑1 (MQ=255)
                                                                                                                                                gCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTAt       >  1:186807/1‑150 (MQ=255)
                                                                                                                                                 cAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTAta      >  1:1507487/1‑149 (MQ=255)
                                                                                                                                                     gTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCaaaa  <  1:1359599/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTCAATATTAGAAGCGGGTTCTGCTAAAAAATAATTGCAACGTACCGGATAAAACCAGCGTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAA  >  NC_000913/3635859‑3636157

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