Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,253,285 T→A 6.3% intergenic (+18/+33) yhaI → / ← yhaJ DUF805 family inner membrane protein/LysR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,253,2850TA 6.3% 336.0 / 7.5 115intergenic (+18/+33)yhaI/yhaJDUF805 family inner membrane protein/LysR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (49/58);  new base A (5/3);  total (54/61)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 119 total aligned reads displayed.

GTTCAGGTGCCTGGGCCTTGCTCTTTTTTGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTCCGGGCTGAC  >  NC_000913/3253137‑3253434
                                                                                                                                                    |                                                                                                                                                     
gTTCAGGTGCCTGGGCCTTGCTCTTTTTTGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTATGGTTCAAATTAATCTTAGAATTGGGGCGata                                                                                                                                                      >  2:1210587/1‑150 (MQ=255)
gTTCAGGTGCCTGGGCCTTGCTCTTTTTTGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTAATTCTGGCAGCAATCGTTACGGAACCTATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGata                                                                                                                                                      >  1:730050/1‑150 (MQ=255)
      gTGCCTGGGCCTTGCTCTTTTTTGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGc                                                                                                                                                <  1:1423511/150‑1 (MQ=255)
        gCCTGGGCCTTGCTCTTTTTTGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGccc                                                                                                                                              <  1:1301405/150‑1 (MQ=255)
          cTGGGCCTTGCTCTTTTTTGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGAAATATCGcccca                                                                                                                                            <  1:1526459/150‑2 (MQ=255)
          cTGGGCCTTGCTCTTTTTTGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGAAATATCGcccca                                                                                                                                            <  1:9131/150‑2 (MQ=255)
            gggCCTTGCTCTTTTTTGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGAAATATCGccccaat                                                                                                                                          <  1:1395980/150‑4 (MQ=255)
             ggCCTTGCTCTTTTTTGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCtttt                                                                                                                                         >  2:256346/1‑150 (MQ=255)
              gCCTTGCTCTTTTTTGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCttttt                                                                                                                                        >  1:1620649/1‑150 (MQ=255)
              gCCTTGCTCTTTTTTGTTACGTTTATTCGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAAAGATCCGAAGTTTGGTTCAAATTAATCATAGAATTGGGGCGATATTTCGCCCCttttt                                                                                                                                        >  1:1091858/1‑150 (MQ=255)
               ccTTGCTCTTTTTTGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTtta                                                                                                                                       >  2:1687968/1‑150 (MQ=255)
                 ttGCTCTTTTTTGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTttatt                                                                                                                                     <  1:1298406/150‑1 (MQ=255)
                   gCTCTTTTTTGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTaa                                                                                                                                   <  1:839590/150‑1 (MQ=255)
                    ctctTTTTTGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAAc                                                                                                                                  <  1:393442/150‑1 (MQ=255)
                      ctTTTTTGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACaa                                                                                                                                <  1:1053048/150‑1 (MQ=255)
                        tttttGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTTGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACaata                                                                                                                              <  1:66742/150‑1 (MQ=255)
                           ttGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAataata                                                                                                                           <  2:600973/150‑1 (MQ=255)
                             gTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATaataaa                                                                                                                         <  1:1434346/150‑1 (MQ=255)
                             gTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATaataaa                                                                                                                         >  2:272004/1‑150 (MQ=255)
                               tCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATaataaata                                                                                                                       >  1:895341/1‑150 (MQ=255)
                                  gTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAAtatta                                                                                                                    >  2:391102/1‑150 (MQ=255)
                                  gTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCAGAAGTTTGGTTCAAATTAATCTTAGAATTGGGTCGAAATTTATCACCTTTTTATTAACAATAATAAAtttta                                                                                                                    >  1:71433/1‑150 (MQ=255)
                                    ttATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATattatt                                                                                                                  >  2:1330851/1‑150 (MQ=255)
                                       ttGGCTGGCTCGTGCTTCTGGTTTTTTTCTTCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTTGTTATAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTc                                                                                                               <  2:1061417/150‑1 (MQ=255)
                                       ttGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTc                                                                                                               <  2:1178369/150‑1 (MQ=255)
                                        tggctggcTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTcc                                                                                                              <  2:1389760/150‑1 (MQ=255)
                                         ggctggctCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCg                                                                                                             >  1:1261642/1‑150 (MQ=255)
                                         ggctggctCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCg                                                                                                             <  2:1170621/150‑1 (MQ=255)
                                          gctggctCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGt                                                                                                            <  2:588612/150‑1 (MQ=255)
                                                    gCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGttt                                                                                                  <  2:26380/150‑1 (MQ=255)
                                                    gCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGttt                                                                                                  <  1:941788/150‑1 (MQ=255)
                                                     cTTCAGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGtttt                                                                                                 >  1:677256/1‑149 (MQ=255)
                                                      ttCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTgg                                                                                                <  2:569695/150‑1 (MQ=255)
                                                          ggTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCAGTTAAAAAGTTTGGGAAt                                                                                            >  1:600672/1‑150 (MQ=255)
                                                           gTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAAtt                                                                                           >  1:562465/1‑150 (MQ=255)
                                                             ttttttCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTc                                                                                         <  1:511569/150‑1 (MQ=255)
                                                                tttCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCg                                                                                      <  2:541051/150‑1 (MQ=255)
                                                                tttCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCg                                                                                      <  1:175286/150‑1 (MQ=255)
                                                                 ttCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGc                                                                                     <  2:1467378/150‑1 (MQ=255)
                                                                  tCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCa                                                                                    >  1:490629/1‑150 (MQ=255)
                                                                   atGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAg                                                                                   <  2:160149/149‑1 (MQ=255)
                                                                     gCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGac                                                                                 >  2:1355514/1‑150 (MQ=255)
                                                                     gCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGac                                                                                 >  2:719875/1‑150 (MQ=255)
                                                                      cacaGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAgccc                                                                                >  2:820601/1‑147 (MQ=255)
                                                                       acaGAAGGTACTTCTGGCAGCAATCGATACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGacac                                                                               >  1:243664/1‑150 (MQ=255)
                                                                        caGAAGGTACTTCTTGCAGCAATCGTTACGGAAACTATCCGAAGTTTGGTTCAAATTTATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACAcc                                                                              <  2:1156099/150‑1 (MQ=255)
                                                                           aaGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGATGTTTGGTTCAAATTAATCTTAGAATTGGGGCGAAATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCAGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAgg                                                                           >  2:118166/1‑150 (MQ=255)
                                                                           aaGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAgg                                                                           <  1:1475635/150‑1 (MQ=255)
                                                                             ggTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGAt                                                                         >  1:788653/1‑150 (MQ=255)
                                                                              gTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGAtt                                                                        <  2:895341/150‑1 (MQ=255)
                                                                               tACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGAttt                                                                       >  2:1494096/1‑150 (MQ=255)
                                                                                aCTTCTGGCAGCAATCTTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATATTAGAATTGGGTCTATATTTAGCCCGTTTTTATTAACAATAATAAATATTATTTTCCGTTATAAAGTTTGGGAATTTCCCGCAGACACCAGGAtttt                                                                      <  1:1210587/150‑2 (MQ=255)
                                                                                  ttCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGc                                                                    <  2:1441354/150‑1 (MQ=255)
                                                                                  ttCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGc                                                                    <  2:677256/150‑1 (MQ=255)
                                                                                       gcagcaATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTATCAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTccc                                                               >  2:1562293/1‑150 (MQ=255)
                                                                                        ctgcaATAGTTACGGAAACGATCCGAATTTTGGTTCAAATTAATCTTTGAATTTGTGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCTTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTcccc                                                              <  2:1278825/148‑1 (MQ=255)
                                                                                        cagcaATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTcccc                                                              <  1:1494096/150‑1 (MQ=255)
                                                                                              tCGTTACGGAAAAGATACGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCTATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACt                                                        <  2:71433/150‑1 (MQ=255)
                                                                                               cGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTg                                                       >  1:1016661/1‑150 (MQ=255)
                                                                                               cGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATATCGCCCCTTTTTCTTAACAATACTAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTg                                                       <  1:1204345/150‑1 (MQ=255)
                                                                                                 ttACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGCGTTCGCTTACCCCATACTGTc                                                     >  2:321735/1‑150 (MQ=255)
                                                                                                 ttACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTc                                                     >  2:1412155/1‑150 (MQ=255)
                                                                                                  tACGGAAACGATCAGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAGTTTGGGAATTTCCCTAAGACACAAGGATTTCGCTTCCCCCATACTGTCac                                                   >  1:950969/1‑150 (MQ=255)
                                                                                                   aCGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCac                                                   <  2:1620649/150‑1 (MQ=255)
                                                                                                    cGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCacg                                                  <  1:583373/150‑1 (MQ=255)
                                                                                                    cGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCacg                                                  >  1:739284/1‑150 (MQ=255)
                                                                                                    cGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCACGCAGACACCAGGATTTCGCTTCCCCCATACTGTCacg                                                  >  1:36975/1‑150 (MQ=255)
                                                                                                       aaaCGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCacgacg                                               <  2:347022/150‑1 (MQ=255)
                                                                                                        aaCGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGc                                              >  1:1676967/1‑150 (MQ=255)
                                                                                                        aaCGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGc                                              <  1:391102/150‑1 (MQ=255)
                                                                                                          cGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGATTCCCCCATACTGTCACGACGCCa                                            >  1:1129760/1‑150 (MQ=255)
                                                                                                          cGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATATTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCa                                            <  1:256346/150‑1 (MQ=255)
                                                                                                             tCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGc                                         >  2:1365098/1‑150 (MQ=255)
                                                                                                              ccGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGcc                                        >  1:1409231/1‑150 (MQ=255)
                                                                                                               cGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCa                                       >  1:475777/1‑150 (MQ=255)
                                                                                                                 aaGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCata                                     >  2:469800/1‑150 (MQ=255)
                                                                                                                  aGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCataa                                    <  1:1094830/150‑1 (MQ=255)
                                                                                                                  aGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCACAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCataa                                    >  2:742101/1‑150 (MQ=255)
                                                                                                                    tttGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCataata                                  <  2:334537/150‑1 (MQ=255)
                                                                                                                       ggTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCataataata                               >  2:957262/1‑150 (MQ=255)
                                                                                                                          tCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCg                            <  2:490629/150‑1 (MQ=255)
                                                                                                                          tCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCg                            >  1:1355289/1‑150 (MQ=255)
                                                                                                                           cAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGa                           >  2:368989/1‑150 (MQ=255)
                                                                                                                            aaaTTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCTGGCCATAATAATATCGAt                          <  2:191022/150‑1 (MQ=255)
                                                                                                                            aaaTTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCAGGCCATAATAATATCGAtc                         <  2:953436/150‑1 (MQ=255)
                                                                                                                            aaaTTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTAGGTTAAAAAGTTTTGGAATTTCCCGCAGACACCAGGTTTTCGCTTCCCCCATACTTTCACGACGCCAGGCCATAATAATATCGAt                          <  2:1091858/150‑1 (MQ=255)
                                                                                                                               ttAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGAtctc                       >  2:1356752/1‑150 (MQ=255)
                                                                                                                                  aTCTTATAATTTGGGCGATATTTCGCCACTGTTTATTAACAATAATAAATATTATTTTCCGTTAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTg                   <  2:950969/150‑1 (MQ=255)
                                                                                                                                  aTCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCt                    <  1:233585/150‑1 (MQ=255)
                                                                                                                                  aTCTTAGAATTGGGGCGATATTTAGCCCATTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCt                    <  2:101468/150‑1 (MQ=255)
                                                                                                                                    cTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCACCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTgg                  >  1:430247/1‑150 (MQ=255)
                                                                                                                                           ttGGGGCGATATTTCGCCCCTTTTTATTATCAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTc           <  1:1562293/150‑1 (MQ=255)
                                                                                                                                           aaggggCGAAATATCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTc           >  2:586324/3‑150 (MQ=255)
                                                                                                                                            aggggCGAAATATCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTcc          >  2:609446/2‑150 (MQ=255)
                                                                                                                                             tggtcGATATTTCGCACCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTCCg         <  1:308198/146‑1 (MQ=255)
                                                                                                                                               ggCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTCCggg       >  1:1105856/1‑150 (MQ=255)
                                                                                                                                               ggCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTCCggg       >  2:1409860/1‑150 (MQ=255)
                                                                                                                                                gCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTCCGGGc      >  1:775445/1‑150 (MQ=255)
                                                                                                                                                gCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTCCGGGc      >  2:954556/1‑150 (MQ=255)
                                                                                                                                                    tatTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTCCGGGCTgac  <  1:469800/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTCAGGTGCCTGGGCCTTGCTCTTTTTTGTTCCGTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTCCGGGCTGAC  >  NC_000913/3253137‑3253434

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