Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 832,364 A→G 6.1% intergenic (+128/+104) 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,3640AG 6.1% 925.7 / 33.5 281intergenic (+128/+104)rhlE/ybiAATP‑dependent RNA helicase/DUF1768 family protein
Reads supporting (aligned to +/- strand):  ref base A (93/169);  new base G (4/13);  total (97/184)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 287 total aligned reads displayed.

AAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAG  >  NC_000913/832216‑832513
                                                                                                                                                    |                                                                                                                                                     
 aaCCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAg                                                                                                                                                     <  1:2902173/150‑1 (MQ=255)
 aaCCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAg                                                                                                                                                     <  1:1903499/150‑1 (MQ=255)
    cTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCg                                                                                                                                                  <  1:3392680/150‑1 (MQ=255)
    cTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCg                                                                                                                                                  <  1:1027437/150‑1 (MQ=255)
      gCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGtt                                                                                                                                                <  2:1644838/150‑1 (MQ=255)
           tGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                           >  2:4032066/1‑150 (MQ=255)
           tGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                           <  1:38495/150‑1 (MQ=255)
            gcgcAGTAATCTTTTAGGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                          <  1:3647191/150‑1 (MQ=255)
             cgcAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                         <  1:2224745/150‑1 (MQ=255)
                aGTAATCTCTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAACACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                      <  1:3848746/150‑5 (MQ=25)
                 gTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCACAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTAAAcgccttatc                                                                                                                                     <  2:4041271/150‑6 (MQ=25)
                     tCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                 <  1:2849350/150‑1 (MQ=21)
                     tCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                 <  2:2426719/150‑1 (MQ=21)
                       ttttATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATATGCATGCTTTTGTAGGCCGGATAAGGCGTTTACTCCGCATCCGGCAt                                                                                                                               >  2:3944484/1‑150 (MQ=255)
                        tttATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAtt                                                                                                                              <  1:269693/150‑1 (MQ=255)
                           ctGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTTCCTGATGCGCTTCGCTTATCAGGCCTACGTGAACTCTTCTATATATTGAATTTGCATGCTTTTGTAGTCCGGATAAGGCGTTTACGCCTCATCCGGCATTTca                                                                                                                           <  2:3482477/149‑1 (MQ=21)
                           aTGCCGGGCTTTGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTca                                                                                                                           <  1:4051438/150‑1 (MQ=255)
                           aTGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGAAAAAGTGCCTGTTGCTCTAGGCTTATCTGGCCTACTTGAACTCTGCAATATTTTTAATTTGCTTGCTTTTGTAGGCGGGTTTAGGGGTTTACGCCGCATCCGGCATTTca                                                                                                                           <  1:1969804/150‑1 (MQ=21)
                              ccGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCacaa                                                                                                                        <  2:2051373/150‑1 (MQ=255)
                                   cTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGc                                                                                                                   <  1:4105954/150‑1 (MQ=255)
                                    tATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCa                                                                                                                  <  1:1377599/150‑1 (MQ=255)
                                       gCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACtt                                                                                                               <  2:2208540/150‑1 (MQ=255)
                                         ccGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGt                                                                                                             <  1:143194/150‑1 (MQ=255)
                                           ggCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGACTTTGCATGCTTTGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCa                                                                                                           <  2:1121155/150‑1 (MQ=255)
                                           ggCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCa                                                                                                           <  2:3275005/150‑1 (MQ=255)
                                             cATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGc                                                                                                         <  1:1290380/150‑1 (MQ=255)
                                              aTCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCa                                                                                                        >  2:3149999/1‑150 (MQ=255)
                                              aTCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCa                                                                                                        >  1:895831/1‑150 (MQ=255)
                                               tCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCaa                                                                                                       <  2:937468/150‑1 (MQ=255)
                                                   gCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTg                                                                                                   >  1:1570639/1‑150 (MQ=255)
                                                    cTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTga                                                                                                  <  1:1176733/150‑1 (MQ=255)
                                                    cTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTga                                                                                                  <  2:285535/150‑1 (MQ=255)
                                                    cTGAGGAACAAACGCATAACTGCGTGATGCGCTACGCTTATCAGGCCTTCGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTga                                                                                                  <  2:1444759/150‑1 (MQ=255)
                                                     tgatgaACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTgag                                                                                                 >  1:3691520/1‑150 (MQ=255)
                                                       ctgaACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGaa                                                                                               <  1:652216/149‑1 (MQ=255)
                                                         gaACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAaca                                                                                             <  1:2572046/150‑1 (MQ=255)
                                                          aacacccGCAACACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATCTTGAATGTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAacac                                                                                            <  2:2247528/144‑1 (MQ=21)
                                                          aaCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAacac                                                                                            >  1:529992/1‑150 (MQ=255)
                                                           aCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACg                                                                                           <  2:1623635/150‑1 (MQ=255)
                                                            cAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATCTTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACgg                                                                                          <  1:2775610/150‑1 (MQ=255)
                                                             aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACggg                                                                                         >  2:2978288/1‑150 (MQ=255)
                                                                 gCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGaaaa                                                                                     <  1:3085449/150‑1 (MQ=255)
                                                                     aaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAAtt                                                                                 <  2:129393/150‑1 (MQ=255)
                                                                      aCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATGTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAAttt                                                                                <  1:464317/150‑1 (MQ=255)
                                                                        tGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTAt                                                                              >  1:303227/1‑150 (MQ=255)
                                                                            tGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGt                                                                          <  2:1054320/150‑1 (MQ=255)
                                                                             gATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGtt                                                                         >  2:306058/1‑150 (MQ=255)
                                                                             gATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGtt                                                                         >  1:869284/1‑150 (MQ=255)
                                                                             gATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGtt                                                                         >  2:514278/1‑150 (MQ=255)
                                                                               tGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGtttt                                                                       >  2:3015756/1‑150 (MQ=255)
                                                                                  gcTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTAcc                                                                    <  2:1739880/150‑1 (MQ=255)
                                                                                  gcTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTAcc                                                                    <  1:961769/150‑1 (MQ=255)
                                                                                  gatACGATTATCAGGCCTACGTGAACTCTGCAAGATATTGAATTTGCATGCTTTTGTAGGCCGGATTTGGCGTTTACGCCGCTTCCGGCATTTCACAACAAGCTCTTGTCAGCTATTTGAGAACACGGGAAAAAATTTATCTGTTTTAcc                                                                    <  1:3944484/148‑1 (MQ=37)
                                                                                   cTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCg                                                                   >  1:760893/1‑150 (MQ=255)
                                                                                   cTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCg                                                                   >  1:247999/1‑150 (MQ=255)
                                                                                       gCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCa                                                               >  2:3566216/1‑150 (MQ=255)
                                                                                       gCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCa                                                               >  1:3111683/1‑150 (MQ=255)
                                                                                        cTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCAt                                                              >  1:1877257/1‑150 (MQ=255)
                                                                                        cTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTGTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTCTCTGTTTTACCGcccct                                                              >  2:3895492/1‑148 (MQ=255)
                                                                                           aTCAGGCCTACGTGAACTCTGCAATATCTTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAg                                                           <  2:772123/150‑1 (MQ=255)
                                                                                            tCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGc                                                          <  2:1541071/150‑1 (MQ=255)
                                                                                               ggCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCAc                                                       <  1:2267831/150‑1 (MQ=255)
                                                                                                     cGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGAc                                                 <  1:3149271/150‑1 (MQ=255)
                                                                                                      gTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACg                                                >  1:809675/1‑150 (MQ=255)
                                                                                                         aaCTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTa                                             <  1:1061268/150‑1 (MQ=255)
                                                                                                         aaCTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGATa                                             >  1:110739/1‑150 (MQ=255)
                                                                                                            tctGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACtt                                          >  2:6009/1‑150 (MQ=255)
                                                                                                            tctGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACtt                                          <  2:2120733/150‑1 (MQ=255)
                                                                                                              tGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTtct                                        >  1:1261919/1‑150 (MQ=255)
                                                                                                              tGCAATATATTGAATGTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTtct                                        <  2:3307425/150‑1 (MQ=255)
                                                                                                               gCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTtctc                                       <  1:2391037/150‑1 (MQ=255)
                                                                                                               gCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTtctc                                       >  2:2073244/1‑150 (MQ=255)
                                                                                                               gCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTtctc                                       <  1:2125724/150‑1 (MQ=255)
                                                                                                                cAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTctct                                      >  2:1369265/1‑150 (MQ=255)
                                                                                                                cAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTctct                                      <  2:3447363/150‑1 (MQ=255)
                                                                                                                cAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTctct                                      <  2:1377604/150‑1 (MQ=255)
                                                                                                                cAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTctct                                      <  2:3148710/150‑1 (MQ=255)
                                                                                                                cAATATACTGACTTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTctct                                      <  2:1962064/150‑1 (MQ=255)
                                                                                                                  atatatTGAATGTGCATGCGTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCtat                                    <  1:1705988/150‑1 (MQ=255)
                                                                                                                   tagatTGACTGTGCATGCTTTTGGAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCtata                                   <  2:2336373/147‑1 (MQ=255)
                                                                                                                     tatTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGc                                 <  1:110070/150‑1 (MQ=255)
                                                                                                                       ttGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCa                               >  1:1631865/1‑150 (MQ=255)
                                                                                                                       ttGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCa                               >  1:1862181/1‑150 (MQ=255)
                                                                                                                             ttGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGtt                         <  1:2808708/150‑1 (MQ=255)
                                                                                                                                 aTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTcgcg                     >  1:3751925/1‑150 (MQ=255)
                                                                                                                                  tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCgcgc                    >  2:3742023/1‑150 (MQ=255)
                                                                                                                                   gagtttGTAGTCCGGCTAAGTATTTTACTCCGAATACTTCATTTAACAACAATCACTTGTCAGCAATTTGAGAACACGTGAAAAAATTTATCTTTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCa                   <  2:2677689/147‑1 (MQ=255)
                                                                                                                                    cTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCaa                  <  1:1795724/150‑1 (MQ=255)
                                                                                                                                       ttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTc               <  1:2618332/150‑1 (MQ=255)
                                                                                                                                        tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTcc              <  1:4061822/150‑1 (MQ=255)
                                                                                                                                        tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTcc              <  1:3180618/150‑1 (MQ=255)
                                                                                                                                        tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTcc              <  1:2978288/150‑1 (MQ=255)
                                                                                                                                          tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCAt            >  1:3527784/1‑150 (MQ=255)
                                                                                                                                           atgCCGGATGCGGCGTAAACGCCTTATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCAtt           >  2:1224276/3‑150 (MQ=255)
                                                                                                                                           aGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCAtt           <  1:1515823/150‑1 (MQ=255)
                                                                                                                                             gCCGGATAAGGCGAATACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGCGAAAAAATTAATCTGTTATACCTCCCATAAGCCACCGAGACGTTACATCTCTATACCCAATTATTCCCGCAACACCATTaa         >  2:4113966/1‑150 (MQ=255)
                                                                                                                                               cGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTaaaa       >  2:893440/1‑150 (MQ=255)
                                                                                                                                                 gATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAgg     >  1:3327485/1‑150 (MQ=255)
                                                                                                                                                   tAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTa   >  2:3333779/1‑150 (MQ=255)
                                                                                                                                                    aaGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAg  <  1:3963270/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAG  >  NC_000913/832216‑832513

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: