Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 832,378 G→T 6.5% intergenic (+142/+90) 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,3780GT 6.5% 311.1 / 8.4 115intergenic (+142/+90)rhlE/ybiAATP‑dependent RNA helicase/DUF1768 family protein
Reads supporting (aligned to +/- strand):  ref base G (56/51);  new base T (2/6);  total (58/57)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01

Only 100 of 119 total aligned reads displayed.

CAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGCCCAGTCGATTCTT  >  NC_000913/832231‑832527
                                                                                                                                                   |                                                                                                                                                     
cAGTAATCTTTTATGCCGGGCTATGACCGGCATCAGGCTGATGAACAAACTCAAAACTGCCTGATGCGCTACGATTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACTCCGCa                                                                                                                                                     >  2:242426/1‑150 (MQ=255)
   tAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                                  <  1:1750969/150‑1 (MQ=25)
   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:1553131/150‑1 (MQ=255)
           tATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTc                                                                                                                                          <  1:1654088/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)
                    cTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGc                                                                                                                                 <  1:1324524/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                                                                                                                       <  2:1898618/150‑1 (MQ=255)
                              cATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGc                                                                                                                       <  1:1420229/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)
                                         tgaACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAac                                                                                                            >  1:1652526/1‑150 (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)
                                                  gCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGaaaa                                                                                                   >  2:539633/1‑150 (MQ=255)
                                                       aCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAAttt                                                                                              <  1:1806793/150‑1 (MQ=255)
                                                       aCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAAttt                                                                                              >  1:52712/1‑150 (MQ=255)
                                                        cTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTa                                                                                             >  1:564838/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)
                                                                tGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATATGtttt                                                                                     >  1:146094/1‑150 (MQ=255)
                                                                  cgcTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCGTCCGGCATTTCACAACAAGCACTTGACAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTAc                                                                                   >  1:733994/1‑150 (MQ=255)
                                                                       cGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGccc                                                                              >  1:1669905/1‑150 (MQ=255)
                                                                       cGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGccc                                                                              >  2:1920854/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)
                                                                              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)
                                                                                              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)
                                                                                                                gCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTcg                                     >  1:1888355/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:1334237/4‑150 (MQ=255)
                                                                                                                           aatgCCGGATGCGGCGTAAACGCCTTATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCAt                          >  2:591920/4‑150 (MQ=255)
                                                                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTa                        >  1:1993153/1‑150 (MQ=255)
                                                                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTa                        >  1:706530/1‑150 (MQ=255)
                                                                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTa                        >  2:768237/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             <  2:1472635/150‑1 (MQ=255)
                                                                                                                                         cGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCTTAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGcaca            >  2:315689/1‑150 (MQ=255)
                                                                                                                                                gccgcATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAGGTAGCCCAGTCGAtta    >  1:1593098/1‑149 (MQ=255)
                                                                                                                                                 ccgcATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGCCCAGTCGATTc    >  1:710539/1‑150 (MQ=255)
                                                                                                                                                 ccgcATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGCCCAGTCGATTc    <  2:395494/150‑1 (MQ=255)
                                                                                                                                                   gcATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTACACTATTGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGCCAAGTCGcttctt  >  2:1216296/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGCCCAGTCGATTCTT  >  NC_000913/832231‑832527

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