Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 832,371 T→A 5.9% intergenic (+135/+97) rhlE → / ← ybiA ATP‑dependent RNA helicase/DUF1768 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913832,3710TA 5.9% 367.3 / 9.4 119intergenic (+135/+97)rhlE/ybiAATP‑dependent RNA helicase/DUF1768 family protein
Reads supporting (aligned to +/- strand):  ref base T (56/56);  new base A (2/5);  total (58/61)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

Only 100 of 120 total aligned reads displayed.

CGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGCCCAGTC  >  NC_000913/832224‑832520
                                                                                                                                                   |                                                                                                                                                     
cGCTGCGCAGTAATCTTTTCTGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTa                                                                                                                                                     <  1:1754867/150‑1 (MQ=255)
cGCTGCGCAGTAATCTTTTATGCCGGGCTATACCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTa                                                                                                                                                     <  1:1867623/150‑1 (MQ=255)
   tGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                                  >  1:202258/1‑150 (MQ=255)
   tGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                                  <  1:1513515/150‑1 (MQ=255)
       cAGTAATCTTTTATGCCGGGCTATGACCGGCATCAGGCTGATGAACAAACTCAAAACTGCCTGATGCGCTACGATTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACTCCGCa                                                                                                                                              >  2:242426/1‑150 (MQ=255)
          tAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                           <  2:1533482/150‑1 (MQ=25)
          aaaTCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                           <  1:1956506/149‑1 (MQ=25)
           aaTCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                          >  1:379265/1‑150 (MQ=255)
             tCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                        <  1:1059000/150‑1 (MQ=21)
             tCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                        <  1:1299311/150‑1 (MQ=21)
                  tATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTc                                                                                                                                   <  1:1654088/150‑1 (MQ=255)
                  tATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTc                                                                                                                                   <  1:1553131/150‑1 (MQ=255)
                   aTGCCGGGCTATGCCCGGCATCATGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTca                                                                                                                                  <  1:1290074/150‑1 (MQ=255)
                   aTGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTca                                                                                                                                  >  2:131627/1‑150 (MQ=255)
                    tGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTcac                                                                                                                                 <  1:347474/150‑1 (MQ=255)
                    tGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTcac                                                                                                                                 >  1:1472635/1‑150 (MQ=255)
                     gCCGGGCTATGCCCGGCATCAGGCTGATTAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTcaca                                                                                                                                <  1:139155/150‑1 (MQ=255)
                      ccGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCacaa                                                                                                                               >  2:720475/1‑150 (MQ=255)
                      ccGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCacaa                                                                                                                               >  2:93799/1‑150 (MQ=255)
                        gggCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCacaaca                                                                                                                             <  2:1121123/150‑1 (MQ=255)
                          gCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAg                                                                                                                           <  1:1440184/150‑1 (MQ=255)
                                cccGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTg                                                                                                                     >  2:1610077/1‑150 (MQ=255)
                                  cGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTc                                                                                                                   <  2:855649/150‑1 (MQ=255)
                                    gCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGAACTTGTCAg                                                                                                                 >  2:372516/1‑150 (MQ=255)
                                     cATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGc                                                                                                                <  1:1420229/150‑1 (MQ=255)
                                     cATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGc                                                                                                                <  2:1898618/150‑1 (MQ=255)
                                       acAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCAAATCAGGCCTACGTGACCTCTGCAATATATTGAATTTGCATGCTTATGTAGGCCGGATAAGGCGTTTACGCCGCATACGGCATTTCACAACAAGCACTTGTCAGCaa                                                                                                              >  1:1079128/2‑150 (MQ=21)
                                         aGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAAtt                                                                                                            <  1:416195/150‑1 (MQ=255)
                                          ggCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAAttt                                                                                                           <  2:1191047/150‑1 (MQ=255)
                                            cTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTga                                                                                                         <  2:131201/150‑1 (MQ=255)
                                            cTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTga                                                                                                         <  2:940433/150‑1 (MQ=255)
                                             tgatgaACAAACGCAAAACTGCCTGATGCGCTACGCTTCTCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTgag                                                                                                        <  1:586804/150‑1 (MQ=255)
                                              gatgaACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTgaga                                                                                                       >  2:1731531/1‑150 (MQ=255)
                                              gatgaACAAACGAAAAACTGCCTGATGCGCTAAGCTTATAATTCATACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTgaga                                                                                                       <  2:1131227/150‑1 (MQ=25)
                                               atgaACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGaa                                                                                                      <  2:708302/150‑1 (MQ=255)
                                                tgaACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGaaa                                                                                                     >  2:159619/1‑149 (MQ=255)
                                                   aCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACg                                                                                                  <  2:640170/150‑1 (MQ=255)
                                                   aCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGGCAGCAATTTGTGAGCACg                                                                                                  >  1:1869479/1‑150 (MQ=255)
                                                     aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACggg                                                                                                >  2:597169/1‑150 (MQ=255)
                                                      aaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGGCAGCAATTTGAGAACACGGGa                                                                                               >  2:450816/1‑150 (MQ=255)
                                                        cGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGGATCCGGCATTTCACAACAAGCACTTGTCAGCAATATGAGAACACGGGaaa                                                                                             >  1:196872/1‑150 (MQ=255)
                                                        cGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGaaa                                                                                             <  1:1871362/150‑1 (MQ=255)
                                                         gCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGaaaa                                                                                            >  2:539633/1‑150 (MQ=255)
                                                              aCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAAttt                                                                                       <  1:1806793/150‑1 (MQ=255)
                                                              aCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAAttt                                                                                       >  1:52712/1‑150 (MQ=255)
                                                               cTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTa                                                                                      <  1:1959973/150‑1 (MQ=255)
                                                                 gCCTGATGCGCTACGCTGATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATc                                                                                    <  2:103555/150‑1 (MQ=255)
                                                                  ccTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGTATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCt                                                                                   >  1:1063505/1‑150 (MQ=255)
                                                                   cTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTg                                                                                  <  1:717571/150‑1 (MQ=255)
                                                                     gATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGtt                                                                                <  2:676019/150‑1 (MQ=255)
                                                                     gATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGtt                                                                                >  1:1251358/1‑150 (MQ=255)
                                                                       tGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGtttt                                                                              <  2:352537/150‑1 (MQ=255)
                                                                        gcgcTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTa                                                                             <  1:1475828/150‑1 (MQ=255)
                                                                         cgcTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCGTCCGGCATTTCACAACAAGCACTTGACAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTAc                                                                            >  1:733994/1‑150 (MQ=255)
                                                                             aCGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGcc                                                                        <  1:168929/150‑1 (MQ=255)
                                                                              cGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGccc                                                                       >  2:1920854/1‑150 (MQ=255)
                                                                              cGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGccc                                                                       >  1:1669905/1‑150 (MQ=255)
                                                                                 ttATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACAGCCCATa                                                                    >  1:768257/1‑150 (MQ=255)
                                                                                   ttCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAg                                                                  <  1:20124/149‑1 (MQ=255)
                                                                                   aTCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAg                                                                  <  2:213285/150‑1 (MQ=255)
                                                                                   aTCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAg                                                                  <  1:739597/150‑1 (MQ=255)
                                                                                   aTCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAg                                                                  <  1:581937/150‑1 (MQ=255)
                                                                                     cAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGAATACGTCATTTCACAACAAGCACTAGTCAGCAATTAGAGAACACGGGAAAAAATTTATCTCTTTTACCGCCCATAAGcc                                                                >  1:87916/1‑150 (MQ=37)
                                                                                          cTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGATAAAATTTATCTGTTTTCCCGCCCATAAGCCACCga                                                           >  2:508350/1‑150 (MQ=255)
                                                                                           tACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCgag                                                          >  1:1004212/1‑150 (MQ=255)
                                                                                            aCGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCgaga                                                         <  1:810222/150‑1 (MQ=255)
                                                                                             cGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGAc                                                        <  1:1731531/150‑1 (MQ=255)
                                                                                                  aCTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTAc                                                   <  1:1620278/150‑1 (MQ=255)
                                                                                                   ctctGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACt                                                  >  1:1415058/1‑150 (MQ=255)
                                                                                                     ctGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACttg                                                <  1:93799/150‑2 (MQ=255)
                                                                                                     ctGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTtc                                                >  1:723761/1‑150 (MQ=255)
                                                                                                            atatTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAg                                         <  1:1969613/150‑1 (MQ=255)
                                                                                                               ttGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCa                                      >  1:395494/1‑150 (MQ=255)
                                                                                                                 gAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAAt                                    >  2:1189682/1‑150 (MQ=255)
                                                                                                                 gAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAAAAAGCACTTGTCAGCTATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAAt                                    >  2:1470743/1‑150 (MQ=255)
                                                                                                                   aTTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTg                                  >  1:786148/1‑150 (MQ=255)
                                                                                                                          tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCgcgc                           <  2:754949/150‑1 (MQ=255)
                                                                                                                            cTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCaa                         >  1:838490/1‑150 (MQ=255)
                                                                                                                             ttttGTAGGCCGGATAAGGCGTTTACGCCGCTTCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAAc                        <  1:1091165/150‑1 (MQ=255)
                                                                                                                             ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAAc                        <  1:1442662/150‑1 (MQ=255)
                                                                                                                                 gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCa                    >  2:1237604/1‑150 (MQ=255)
                                                                                                                                  tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCAt                   >  1:1706650/1‑150 (MQ=255)
                                                                                                                                  tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCAt                   >  1:324423/1‑150 (MQ=255)
                                                                                                                                  aatgCCGGATGCGGCGTAAACGCCTTATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCAt                   >  2:591920/4‑150 (MQ=255)
                                                                                                                                  aatgCCGGATGCGGCGTAAACGCCTTATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCAt                   >  2:1334237/4‑150 (MQ=255)
                                                                                                                                    ggCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTa                 >  2:768237/1‑150 (MQ=255)
                                                                                                                                    ggCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTa                 >  1:706530/1‑150 (MQ=255)
                                                                                                                                    ggCCGGATAAGGCGTTTACGCCGCATCCGGCAATTCACAACAAGCACTTGTCAGCAATTAGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTCCTTCTCTATAGCCAATTGTTCGCGCAACTCCATTa                 >  2:1312114/1‑150 (MQ=255)
                                                                                                                                     gCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTaa                >  1:256995/1‑150 (MQ=255)
                                                                                                                                       cGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTaaaa              <  1:131627/150‑1 (MQ=255)
                                                                                                                                         gATAAGGCGTTTACGCCTCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATGTGTTTTACGGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAgg            <  2:1079128/150‑1 (MQ=255)
                                                                                                                                           tAAGGCGTTTACGCCGCATCCGGCATTTCACAACATGCACTTGTCAGCAATTTGAGAACACGGGAAAAATGTTATCTGTTTTACCGCGCATAAGACACCGAGACGTTACTTCTCTTTAGCCAATTGTTCGCGCCTCTCCAATAAAAGGTa          >  1:325091/1‑150 (MQ=255)
                                                                                                                                             aGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGc        <  1:621009/150‑1 (MQ=255)
                                                                                                                                             aGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGc        >  2:1733415/1‑150 (MQ=255)
                                                                                                                                              ggCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGcc       <  2:572529/150‑1 (MQ=255)
                                                                                                                                              ggCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCAGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCGAATTGTTCGCGCATCTCCATTAAAAGGTAGcc       >  2:67002/1‑150 (MQ=255)
                                                                                                                                               gCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGccc      <  1:1285960/150‑1 (MQ=255)
                                                                                                                                               gCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGccc      <  2:1472635/150‑1 (MQ=255)
                                                                                                                                                cGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCTTAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGcaca     >  2:315689/1‑150 (MQ=255)
                                                                                                                                                   ttACGCCGCAGCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGACACCGAGACGTTACTTCTCTATAGCCAATTGTACGCGAAACTCAATAAAAAGGTAGCCCAGTc  >  1:1983160/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGCCCAGTC  >  NC_000913/832224‑832520

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: