Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 940,751 T→C 5.1% intergenic (+31/‑208) serS → / → dmsA seryl‑tRNA synthetase/dimethyl sulfoxide reductase, anaerobic, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913940,7510TC 5.1% 422.6 / 8.6 138intergenic (+31/‑208)serS/dmsAseryl‑tRNA synthetase/dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (67/63);  new base C (3/4);  total (71/67)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 138 total aligned reads displayed.

GGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGA  >  NC_000913/940606‑940897
                                                                                                                                                 |                                                                                                                                                  
ggCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCgg                                                                                                                                                >  2:1288096/1‑150 (MQ=255)
ggCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGATGAAGGACGTATAGAAGTACCAGATGTTCTGCGTCCGTATATGAACGGACTGGAATATTTTGTCAAATACCCAATTTTTCTGAATCTAAAAAGCGGCAGCgg                                                                                                                                                >  1:1125356/1‑150 (MQ=255)
 gCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCggg                                                                                                                                               <  1:1449121/150‑1 (MQ=255)
     ttGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCt                                                                                                                                           <  2:1807681/150‑1 (MQ=255)
          cGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttta                                                                                                                                      <  1:936311/150‑2 (MQ=255)
          cGTACGCCGGTTGCAGTAATGGAACACTATCAGCAGGCTGATGGTCGTATTGAAGAACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttta                                                                                                                                      <  2:1704071/150‑2 (MQ=255)
           gTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGAttttata                                                                                                                                     >  1:895753/1‑147 (MQ=255)
            tACGCTGGTTGCAGTAATGGAAAACTATCATCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTg                                                                                                                                    <  2:1392110/150‑1 (MQ=255)
            tACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTg                                                                                                                                    <  2:726578/150‑1 (MQ=255)
            gaCGCTGGTTGCAGTAATGGAAAAGTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGAATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTg                                                                                                                                    <  1:122195/149‑1 (MQ=255)
             aCGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttagat                                                                                                                                   <  1:741030/150‑5 (MQ=255)
              cGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtc                                                                                                                                  <  2:958197/150‑1 (MQ=255)
                cTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtctc                                                                                                                                <  2:1631857/150‑1 (MQ=255)
                cTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttagattca                                                                                                                                <  1:499486/150‑8 (MQ=255)
                  ggTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTccc                                                                                                                              >  1:322434/1‑150 (MQ=255)
                  ggTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTccc                                                                                                                              <  1:233822/150‑1 (MQ=255)
                     tGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCttt                                                                                                                           <  2:865241/150‑1 (MQ=255)
                      gCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTg                                                                                                                          >  2:953121/1‑150 (MQ=255)
                         gTATTGGAAAACTATCAGCAGGCTGATGGTCGTATTGAATTACCTGAAGTTCTGCGTCCGTATATTAACGTACTGGAATATATTGGCTAATCCCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATa                                                                                                                       <  1:1259491/150‑1 (MQ=255)
                         gTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTtata                                                                                                                       >  2:824506/1‑150 (MQ=255)
                           aaTGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATAcc                                                                                                                     <  1:396087/150‑1 (MQ=255)
                             tGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGa                                                                                                                   >  1:632140/1‑150 (MQ=255)
                             tGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGa                                                                                                                   >  1:973058/1‑150 (MQ=255)
                              ggAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGaa                                                                                                                  >  2:1195287/1‑150 (MQ=255)
                               gAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAAc                                                                                                                 >  2:1088677/1‑150 (MQ=255)
                               gAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAAc                                                                                                                 <  1:89667/150‑1 (MQ=255)
                                    cTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataa                                                                                                            <  2:1412935/150‑1 (MQ=255)
                                        cagcagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGCACGGACTGGAATATTTTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCGTTTTTTGTCTCCCTTTGATACCGAACAATAATTaa                                                                                                        >  2:330187/1‑149 (MQ=255)
                                         agcagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCAATTTTTATTAAAAATAAAAACCGCATCCGGCGGCTTTTTTTTTAACCATTTATTACCGAAAAAAATTTaaa                                                                                                       >  2:1780265/1‑148 (MQ=255)
                                          gcagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActc                                                                                                      <  1:1026560/150‑1 (MQ=255)
                                           cagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcc                                                                                                     <  2:520017/150‑1 (MQ=255)
                                            agGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACTGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcct                                                                                                    <  1:296514/150‑1 (MQ=255)
                                              gCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCa                                                                                                  <  2:1564028/150‑1 (MQ=255)
                                               cTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCAc                                                                                                 <  1:2374/150‑1 (MQ=255)
                                                   tgttCGTATTTAAGTACCATAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTac                                                                                             <  2:1709863/147‑1 (MQ=255)
                                                    ggTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTaca                                                                                            <  1:1704787/150‑1 (MQ=255)
                                                     gTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATTTATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATTATTACTCCTCACTTacac                                                                                           >  1:1564590/1‑150 (MQ=255)
                                                     gTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATTCCGAACAATAATTAATCCGCACTTaccc                                                                                           >  1:1367481/1‑148 (MQ=255)
                                                     gTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTacac                                                                                           >  2:1703110/1‑150 (MQ=255)
                                                        gTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTTATACCGAACAATAATTACTACTAACTTACACGTa                                                                                        >  2:618090/1‑150 (MQ=255)
                                                        gTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTa                                                                                        >  2:133309/1‑150 (MQ=255)
                                                        gTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTa                                                                                        >  1:764284/1‑150 (MQ=255)
                                                         tATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTaa                                                                                       >  2:288695/1‑150 (MQ=255)
                                                             gAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAtact                                                                                   >  2:1620999/1‑150 (MQ=255)
                                                                   ccAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTc                                                                             >  2:1121658/1‑150 (MQ=255)
                                                                   ccAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAAGTACTCCTCACTTACACGTAATACTACTTTc                                                                             >  1:759645/1‑150 (MQ=255)
                                                                    cAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCg                                                                            <  2:322434/150‑1 (MQ=255)
                                                                     aGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGa                                                                           >  2:895724/1‑150 (MQ=255)
                                                                      taaGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAg                                                                          <  1:1611886/149‑1 (MQ=255)
                                                                      gAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATTTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAg                                                                          <  1:1259466/150‑1 (MQ=255)
                                                                       aaGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGt                                                                         <  2:1342114/150‑1 (MQ=255)
                                                                       aaGTTCTGCGTACGTATATGAACGGACTGGAATATATTGGCTAATACCAAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGt                                                                         <  1:283989/150‑1 (MQ=255)
                                                                        aGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTGTTGTATCCCTTTTTTTCCGAACAATAAATAATAATCACTTACCCGTAAAACTAATTTCGAGtt                                                                        >  1:370024/1‑149 (MQ=255)
                                                                         gTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGa                                                                       >  2:1187028/1‑150 (MQ=255)
                                                                         gTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGa                                                                       >  1:1657903/1‑150 (MQ=255)
                                                                         gTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGa                                                                       >  1:1287887/1‑150 (MQ=255)
                                                                                  gagTATATTAACTGACTGGAATATATTGGCTAATACCATATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGATCAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAcc                                                              <  1:776546/148‑1 (MQ=255)
                                                                                   cGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCt                                                             >  2:1132667/1‑150 (MQ=255)
                                                                                    gTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTa                                                            >  2:1441925/1‑150 (MQ=255)
                                                                                    gTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTATGATACCGAACAATAATTACTACTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAACtt                                                            >  1:140062/1‑149 (MQ=255)
                                                                                     tataTGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGTTACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTAt                                                           >  1:1821806/1‑150 (MQ=255)
                                                                                        atGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTAtctc                                                        >  2:419506/1‑150 (MQ=255)
                                                                                         tGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATctct                                                       >  2:571178/1‑150 (MQ=255)
                                                                                           aaCGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCttt                                                     <  2:1084730/150‑1 (MQ=255)
                                                                                            aCGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTg                                                    <  2:1258959/150‑1 (MQ=255)
                                                                                                aCTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGAttt                                                <  2:244996/150‑1 (MQ=255)
                                                                                                  tGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTc                                              <  2:1413487/150‑1 (MQ=255)
                                                                                                   ggAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTGTCCCTTTTATACCGAACAATAATTAATCCTCACTTACACGTAATACTAATTTCGAGTGAAAATCTAACTATCTCTTTGATTTTCa                                             >  1:1540048/1‑150 (MQ=255)
                                                                                                   ggAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCAATCTCTTTGATTTTCa                                             >  1:632597/1‑150 (MQ=255)
                                                                                                   ggAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTACCTTAAACGTAATACTAATTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCa                                             >  1:1156707/1‑150 (MQ=255)
                                                                                                   ggAATATATTGGCTAATACCCAATTTATCTGACTCTAAAAAGCGACTGCGGGCACTTATTTTGTCTCCCTTAGATACAGAACAATAATAACTCCTCACTTACACGTAATACTACTATCGAGTGAAAAACGACCTATCTCTGTGATATTCa                                             >  1:281937/1‑150 (MQ=255)
                                                                                                     aaTATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATTCTACTTTCGAGTGAAAATATACCTATCTCTTTGATTTTaaaa                                           >  1:444565/1‑150 (MQ=255)
                                                                                                     aaTATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTATTTGATTTTaaaa                                           >  1:1824026/1‑150 (MQ=255)
                                                                                                     aaTATATTGGCTAAAACCCAAATTTACTGAATCTAAAAAGAGCCTGCGGGCTCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCaaa                                           <  1:1780110/150‑1 (MQ=255)
                                                                                                     aaTATAATGGGTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTACCTTTGATACCGAACAATAATAACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCaaa                                           <  2:1102256/150‑1 (MQ=255)
                                                                                                          atTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAattatt                                      <  1:1253084/150‑1 (MQ=255)
                                                                                                           ttGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTc                                     <  2:897262/150‑1 (MQ=255)
                                                                                                           ttGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTc                                     <  2:51665/150‑1 (MQ=255)
                                                                                                            tGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCg                                    >  2:1293063/1‑150 (MQ=255)
                                                                                                            tGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGACTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAAATATTCg                                    >  2:934830/1‑150 (MQ=255)
                                                                                                             ggCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTAACTCTTTGATTTTCAAATTATTCGa                                   >  2:1134407/1‑150 (MQ=255)
                                                                                                              gCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGAt                                  >  2:487552/1‑150 (MQ=255)
                                                                                                               cTAATAACAAATTTTTCTTAATCAAAAAAGAGCCTGAGTGCGCTTTTTTTTTCTCCCTTTTATACCGAACAATAATTACTCCTAACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATg                                 <  2:736180/150‑1 (MQ=255)
                                                                                                                tAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTAATTTTCAAATTATTCGATGt                                <  2:1457190/150‑1 (MQ=255)
                                                                                                                   tACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTAAAATTATTCGATGtata                             >  1:1188792/1‑150 (MQ=255)
                                                                                                                    aCCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCAATCTCTTTGATTTTCAAATTATTCGATGTATAc                            <  2:632597/150‑1 (MQ=255)
                                                                                                                     cccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACa                           >  1:505934/1‑150 (MQ=255)
                                                                                                                      ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGACTCACTTAGATACCGAAGAATAATTACTCCTCACTAAAACGGAAAACTACTTTCGAGTGAAAATCTAGCTCTCTCTTAGATTTACAAATTATTCGATGTATACaa                          >  1:256504/1‑150 (MQ=255)
                                                                                                                       cAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAg                         <  2:1454248/150‑1 (MQ=255)
                                                                                                                         aTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTACTCACTTACACGTAATAATACTTTCGAGTGAAAATATACCTATCTCTTTTATTTTAAAATTATTAGATGTATACAAGcc                       >  2:1570516/1‑150 (MQ=255)
                                                                                                                                    caaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGaa            >  2:1531254/4‑150 (MQ=255)
                                                                                                                                     cTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGcagtt           >  1:161742/1‑146 (MQ=255)
                                                                                                                                     cTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATATACCTATCTCTTTGATTTTCAATTTATTCGATGTATACAAGCCTATATAGCTAAc           >  1:1192267/1‑150 (MQ=255)
                                                                                                                                     aaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAAc           >  2:1284739/3‑150 (MQ=255)
                                                                                                                                        aaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGc        <  1:1121658/150‑1 (MQ=255)
                                                                                                                                          aaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCTAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCATATATAGCGAACTGCta      <  2:1151767/150‑1 (MQ=255)
                                                                                                                                           aGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGATCAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCtat     <  1:1691480/150‑1 (MQ=255)
                                                                                                                                           aGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCtat     >  2:650169/1‑150 (MQ=255)
                                                                                                                                             cgcCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAg   <  2:121745/150‑1 (MQ=255)
                                                                                                                                              gcCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTTATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGa  >  1:297865/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGA  >  NC_000913/940606‑940897

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