Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,088,091 G→C 7.6% intergenic (+34/+16) yiiG → / ← frvR DUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,088,0910GC 7.6% 415.8 / 18.1 129intergenic (+34/+16)yiiG/frvRDUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain
Reads supporting (aligned to +/- strand):  ref base G (59/60);  new base C (6/4);  total (65/64)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 135 total aligned reads displayed.

CTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGC  >  NC_000913/4087943‑4088238
                                                                                                                                                    |                                                                                                                                                   
cTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggc                                                                                                                                                    >  2:213655/1‑150 (MQ=255)
cTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggc                                                                                                                                                    <  2:1410909/150‑1 (MQ=255)
   ggATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCa                                                                                                                                                 >  1:716352/1‑150 (MQ=255)
   ggATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCa                                                                                                                                                 >  2:1091285/1‑150 (MQ=255)
    gATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCAt                                                                                                                                                <  1:1203430/150‑1 (MQ=255)
       aaaGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCg                                                                                                                                             <  2:1531033/150‑1 (MQ=255)
         agagCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                                                                           <  1:88329/150‑1 (MQ=255)
          gagCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGaca                                                                                                                                          <  1:445802/150‑1 (MQ=255)
           agCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGacac                                                                                                                                         >  2:25485/1‑150 (MQ=255)
           agCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGacac                                                                                                                                         <  1:1211961/150‑1 (MQ=255)
            gCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGAcaca                                                                                                                                        <  1:897492/150‑1 (MQ=255)
             cAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAg                                                                                                                                       <  2:1142510/150‑1 (MQ=255)
               aCTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGtt                                                                                                                                     >  1:1181454/1‑150 (MQ=255)
                 ttCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCa                                                                                                                                   <  1:282824/150‑1 (MQ=255)
                  tCAGGCCGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCAGTACGCGAGTACAACGCAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGCCACAGTTCAc                                                                                                                                  >  2:1229956/1‑150 (MQ=255)
                  tCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAc                                                                                                                                  <  2:956763/150‑1 (MQ=255)
                   cAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcc                                                                                                                                 >  1:750908/1‑150 (MQ=255)
                    aGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGCCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCAGACACAGTTCAccc                                                                                                                                >  1:1117391/1‑150 (MQ=255)
                    aGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGAGGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAccc                                                                                                                                <  1:1508975/150‑1 (MQ=255)
                           aaaCTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCTTTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCCACACAGATCACCCATtgtta                                                                                                                         >  1:718152/1‑148 (MQ=255)
                            aaCTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACCCAGGTAACCCCGTGTGaa                                                                                                                        >  1:1299127/1‑150 (MQ=255)
                            aaCTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                                                                        >  1:1633229/1‑150 (MQ=255)
                            aaCTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACACCTTGTGaa                                                                                                                        >  1:1227069/1‑150 (MQ=255)
                             aCTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAg                                                                                                                       <  2:1635963/150‑1 (MQ=255)
                              cTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGt                                                                                                                      <  1:246307/150‑1 (MQ=255)
                              cTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGt                                                                                                                      <  2:927866/150‑1 (MQ=255)
                              cTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGacactgca                                                                                                                                     >  2:1244199/1‑131 (MQ=255)
                                cAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTaa                                                                                                                    >  1:704052/1‑150 (MQ=255)
                                  aGCGGGATGGTCGAAGACTCTTTCCCGAGAGCATTACTCGAGTACACCTAAATGGTTGATGACTATAATAGCATGCGTTAAACGTTTTTGATGTTAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACt                                                                                                                  <  1:1592784/150‑1 (MQ=255)
                                     gggAGGGTCTAAGACTGTTTCCCTAGAGCATTACGCGAGTACAACTAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCa                                                                                                               <  2:1322604/149‑1 (MQ=255)
                                       gATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATa                                                                                                             >  2:401678/1‑150 (MQ=255)
                                       gATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAAAAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACATAGTTCACCCCTTGTGAAGTAACTCCATa                                                                                                             <  2:1151490/150‑1 (MQ=255)
                                             cGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTAACCCCTTGTGAAATAACTCCATAAATGcc                                                                                                       >  1:1491865/1‑150 (MQ=255)
                                               aaGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTc                                                                                                     >  1:897034/1‑150 (MQ=255)
                                                aGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCa                                                                                                    <  1:1377620/150‑1 (MQ=255)
                                                aGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCa                                                                                                    <  2:358373/150‑1 (MQ=255)
                                                aGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCa                                                                                                    <  1:1384260/150‑1 (MQ=255)
                                                aGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCa                                                                                                    <  2:1181454/150‑1 (MQ=255)
                                                   gtctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCatat                                                                                                 >  1:547221/2‑150 (MQ=255)
                                                    tctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGACTCATATc                                                                                                >  1:51190/1‑150 (MQ=255)
                                                     ctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATcc                                                                                               <  1:1025318/150‑1 (MQ=255)
                                                         cccGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                                                                           <  1:620172/150‑1 (MQ=255)
                                                         cccGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                                                                           >  2:896457/1‑150 (MQ=255)
                                                             agagCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAcc                                                                                       <  1:844918/150‑1 (MQ=255)
                                                             agagCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTAGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAcc                                                                                       >  2:830677/1‑150 (MQ=255)
                                                                gCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCgg                                                                                    <  2:1364710/150‑1 (MQ=255)
                                                                 cATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                                                                                   <  1:213655/150‑1 (MQ=255)
                                                                  aTTACGCTAGTACAACGAAATGGTTGTTGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCa                                                                                  <  2:718152/150‑1 (MQ=255)
                                                                    tACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACACTTCACACCTTTTGAAGTAAATCCATAAATGCCGCATATACTGTTAACCCGGCaat                                                                                >  1:629798/1‑150 (MQ=255)
                                                                       gcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACACCTTGTGAAGTAACTCCATAAATGCCTCATAACCTGTTAAACCGGCaataat                                                                             >  2:448102/1‑150 (MQ=255)
                                                                        cgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtc                                                                            <  1:25485/150‑1 (MQ=255)
                                                                          aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtctc                                                                          <  1:1236353/150‑1 (MQ=255)
                                                                              cAACGAAATGGTTGATGACTATACTAGCCTGCGTTAAACGTTTTTGATGATAAATGCCGTGTCGGATTACGCGCATCCGCCACAGTTAACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAAACGGCAATAATCTccgcc                                                                      >  2:81186/1‑150 (MQ=255)
                                                                              cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgcc                                                                      <  2:1052341/150‑1 (MQ=255)
                                                                                   aaaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATg                                                                 <  1:369930/150‑1 (MQ=255)
                                                                                       ggTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGATGTAACTCCATAAATGCCTCAAATCCTGTTAACCCGGCAATAATCTCCGCCCGATGCTGa                                                             >  2:1641371/1‑150 (MQ=255)
                                                                                       ggTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa                                                             >  2:25543/1‑150 (MQ=255)
                                                                                        gTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTAGTGAAGTAACTCCATAAATGCCTCATATCCTCTTAACCCGGCAATAATCTCCGCCGGATGCTGac                                                            <  2:954763/150‑1 (MQ=255)
                                                                                           gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGAcaca                                                         <  2:1131842/150‑1 (MQ=255)
                                                                                            atgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTAGGAAGCGGCGCCTCCGACAAAGTTAACCCCTTGTGAAGTAACTCCAAAAATGCCTCATATACTGTTAACCCGGCAATAATAACCCCCGGATGCTGACACaa                                                        >  2:175230/1‑150 (MQ=255)
                                                                                                cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg                                                    >  1:1235204/1‑150 (MQ=255)
                                                                                                  ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCt                                                  <  1:348077/150‑1 (MQ=255)
                                                                                                   taataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACCCAGTTCACCCCTTTTGAATTAACTCCATAAATTCCTCAAATACTGTTAACCCCGCAATAATCTCCGCCGGATGCTGACACAAAATTCtt                                                 >  2:1512324/1‑150 (MQ=255)
                                                                                                    aataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTg                                                <  1:187373/150‑1 (MQ=255)
                                                                                                    aataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGGGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTg                                                <  2:1181171/150‑1 (MQ=255)
                                                                                                     ataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATACGACACAGTTAACACCTTGTGAAGTAACTCCATAAATGCCTCATATGCTGTTAACCCGGCAATAATCTCCGCCGGTTGCTGACACAATATGATTGc                                               >  2:1089547/1‑150 (MQ=255)
                                                                                                      taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGcc                                              <  2:704052/150‑1 (MQ=255)
                                                                                                        gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCACATATCCGGTTAACCCGGCAATAATCTACGCCGGCTGCTGACACAATATGCATGCCac                                            >  2:185382/1‑149 (MQ=255)
                                                                                                         ccTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCata                                           >  2:325089/1‑148 (MQ=255)
                                                                                                           tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATa                                         <  1:404824/150‑1 (MQ=255)
                                                                                                              gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGATAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCt                                      >  2:337948/1‑150 (MQ=255)
                                                                                                                tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTaa                                    >  1:127628/1‑150 (MQ=255)
                                                                                                                 aaaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTaaa                                   <  2:1437791/150‑1 (MQ=255)
                                                                                                                 aaaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTACATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTaaa                                   >  1:1031379/1‑150 (MQ=255)
                                                                                                                  aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAat                                  <  2:1633229/150‑1 (MQ=255)
                                                                                                                  aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAat                                  >  1:228878/1‑150 (MQ=255)
                                                                                                                       ttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATAtttt                             <  2:454938/150‑1 (MQ=255)
                                                                                                                        tttGATTATAAATTCTGTGTCGGATGCGGCGAATCCGACACGGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATAttttt                            <  2:439115/150‑1 (MQ=255)
                                                                                                                        tttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATAttttt                            <  1:277569/150‑1 (MQ=255)
                                                                                                                          tgatgaTAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTACATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAg                          <  1:49075/150‑1 (MQ=255)
                                                                                                                          tgatgaTAAATGCAGTGTCGGATGCGGCGCATACGACACATTTCACCCCTTGTGAAGTAACTACATAAATGCCTCATATACTGTTAACCCGGCAATAATATCCGCAGGATGCTGAAACAATATGATTGCCAGATATCTAAATCTTTTTAg                          >  2:35849/1‑150 (MQ=255)
                                                                                                                             tgatAAATGCAGTGTCGGATGCGGCGCATCCGAGCCAGTTCACGCCTTGTGAAGTAACTCCATAAATGCCTCAAATCCTGTTAACCCGGCAATATTCTCCGCCGGATGCTGACACACTATGATGGCCAGATAGCTAAATATTTTTAGCtg                       >  1:292442/1‑149 (MQ=255)
                                                                                                                              gatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCa                      <  2:29374/150‑1 (MQ=255)
                                                                                                                                   aTGCAGTGTCGGATGCGGCGCATCCGCCACAGTTCAACCCTTGTGAACTAACTCCATTAATGCCTCATATCCTGTTAACCCGTCAATAATCACCGCCGGCTGCTGCCACAATAAGCTTGCCAGAAATCGTAATATTTTTAGATCATGCcg                 >  2:1054522/1‑150 (MQ=255)
                                                                                                                                   aTGCAGTGTCGGATGCGGCGCATCCGACCCAGTTCACCCCTTGTGAAGTAACGCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAACTCCGCCGGAAGCTTCCACAATATGCTGGCCAGATAGAGAAAGAGTTTTTGCACATGAcg                 >  2:1054528/1‑150 (MQ=255)
                                                                                                                                   aTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGAcg                 <  1:25543/150‑1 (MQ=255)
                                                                                                                                     gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGAcgcg               >  1:1621013/1‑150 (MQ=255)
                                                                                                                                     gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGAcgcg               <  2:897034/150‑1 (MQ=255)
                                                                                                                                      cAGTGTCGGATGCGGCGCATCCGACACAGTTCCCCCCTTGTGAAGTAACTCCCTAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACgcgc              >  1:364900/1‑150 (MQ=255)
                                                                                                                                        gtgtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCAt            >  1:1460824/1‑150 (MQ=255)
                                                                                                                                        gtgtCGGATGCGCCGCATCCGACACAGTTAACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCAt            >  1:165250/1‑150 (MQ=255)
                                                                                                                                          gtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcg          >  1:1608271/1‑150 (MQ=255)
                                                                                                                                          gtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcg          >  2:887132/1‑150 (MQ=255)
                                                                                                                                          gtCGGATGCGCCGCATCCGACACAGTTCCCCCCTTGTGAAGTAACTCCATTAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcg          >  1:1257493/1‑150 (MQ=255)
                                                                                                                                          gtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcg          >  1:1354444/1‑150 (MQ=255)
                                                                                                                                            cGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcggc        >  2:1193492/1‑150 (MQ=255)
                                                                                                                                             ggATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcggcg       >  1:896178/1‑150 (MQ=255)
                                                                                                                                             ggATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcggcg       >  1:413287/1‑150 (MQ=255)
                                                                                                                                             ggATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcggcg       >  2:1388030/1‑150 (MQ=255)
                                                                                                                                              gATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcggcgg      >  2:510120/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGC  >  NC_000913/4087943‑4088238

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: