Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,074,842 C→T 5.1% H58H (CAC→CAT) ‡ yihW → putative DNA‑binding transcriptional regulator YihW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,074,8420CT 5.1% 573.9 / 11.5 176H58H (CAC→CAT) ‡yihWputative DNA‑binding transcriptional regulator YihW
Reads supporting (aligned to +/- strand):  ref base C (90/77);  new base T (7/2);  total (97/79)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

Only 100 of 178 total aligned reads displayed.

AACCCGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGACTGT  >  NC_000913/4074693‑4074991
                                                                                                                                                     |                                                                                                                                                     
aaCCCGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGGGCCATCAc                                                                                                                                                       <  2:355988/150‑1 (MQ=255)
  cccGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACgg                                                                                                                                                     >  1:14050/1‑150 (MQ=255)
   ccGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGt                                                                                                                                                    <  2:105309/150‑1 (MQ=255)
              ccAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCg                                                                                                                                         >  2:655178/1‑150 (MQ=255)
               cAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCgg                                                                                                                                        <  1:750608/150‑1 (MQ=255)
                 aCTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGc                                                                                                                                      <  2:168645/150‑1 (MQ=255)
                 aCTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGc                                                                                                                                      <  1:1142687/150‑1 (MQ=255)
                      tcATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCa                                                                                                                                 <  1:1182338/150‑1 (MQ=255)
                       cATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAg                                                                                                                                >  1:401187/1‑150 (MQ=255)
                        aTGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGc                                                                                                                               >  1:180965/1‑150 (MQ=255)
                            tGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcgtcg                                                                                                                           >  2:494512/1‑150 (MQ=255)
                                     aGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCggg                                                                                                                                            >  2:1072244/1‑124 (MQ=255)
                                      gCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGc                                                                                                                 <  1:802506/150‑1 (MQ=255)
                                      gCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGc                                                                                                                 <  1:608860/150‑1 (MQ=255)
                                        gTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGt                                                                                                               >  1:591765/1‑150 (MQ=255)
                                          gggTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTc                                                                                                             >  1:865862/1‑150 (MQ=255)
                                           ggTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCg                                                                                                            <  1:891620/150‑1 (MQ=255)
                                            gTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCTCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGa                                                                                                           <  2:1196581/150‑1 (MQ=255)
                                             tataTGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGag                                                                                                          <  1:272480/150‑1 (MQ=255)
                                             tataTGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGag                                                                                                          >  1:321040/1‑150 (MQ=255)
                                             tataTGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGag                                                                                                          <  1:676187/150‑1 (MQ=255)
                                              atatGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTAGGCGCCATCACGGTGGCGCGGGTAGGGCTTCCAGCGTCGTTAATACGGCGTTCGagc                                                                                                         <  1:1276427/150‑1 (MQ=255)
                                               tatGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGagca                                                                                                        <  2:687966/150‑1 (MQ=255)
                                                 tGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGagcagc                                                                                                      >  2:441203/1‑150 (MQ=255)
                                                  gAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCg                                                                                                     >  2:358901/1‑150 (MQ=255)
                                                      atTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAg                                                                                                 >  2:85353/1‑150 (MQ=255)
                                                       ttGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAgg                                                                                                >  1:79461/1‑150 (MQ=255)
                                                              gCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCa                                                                                         <  1:275310/150‑1 (MQ=255)
                                                               cTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCaa                                                                                        >  1:604775/1‑150 (MQ=255)
                                                                tGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCaaa                                                                                       >  2:103536/1‑150 (MQ=255)
                                                                 ggCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAAc                                                                                      >  2:1213563/1‑150 (MQ=255)
                                                                   cAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCg                                                                                    <  2:1177020/150‑1 (MQ=255)
                                                                          gctgGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGaaaa                                                                             >  2:597619/1‑150 (MQ=255)
                                                                           ctgGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGaaaaa                                                                            >  2:1139620/1‑150 (MQ=255)
                                                                           ctgGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGaaaaa                                                                            >  1:1155598/1‑150 (MQ=255)
                                                                             ggATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGaaaaaaa                                                                          >  2:218799/1‑150 (MQ=255)
                                                                              gATGTCTCCCCGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGaaaaaaaa                                                                         >  2:823645/1‑150 (MQ=255)
                                                                                tGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGc                                                                       <  1:1029757/150‑1 (MQ=255)
                                                                                   ctcGACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGAt                                                                    <  2:1223913/150‑1 (MQ=255)
                                                                                   ctcCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGAt                                                                    >  2:610631/1‑150 (MQ=255)
                                                                                   ctcCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGAt                                                                    >  2:260651/1‑150 (MQ=255)
                                                                                   ctcCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGAt                                                                    >  2:565859/1‑150 (MQ=255)
                                                                                      cACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGc                                                                 >  2:396698/1‑150 (MQ=255)
                                                                                       aCGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGcc                                                                >  1:515731/1‑150 (MQ=255)
                                                                                        cGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCg                                                               >  2:112159/1‑150 (MQ=255)
                                                                                         gCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGa                                                              <  1:121947/150‑1 (MQ=255)
                                                                                           aGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAg                                                            <  1:655178/150‑1 (MQ=255)
                                                                                              cGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCCCGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGAATGCCGAAGCgg                                                         >  2:1109754/1‑150 (MQ=255)
                                                                                                 tCCGCCGGGATATTCGTAAATTAAGCTAGCTAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGACGAAGCGGTgg                                                      <  2:862977/150‑1 (MQ=255)
                                                                                                 tCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTgg                                                      <  2:706525/150‑1 (MQ=255)
                                                                                                 tCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTgg                                                      <  1:301088/150‑1 (MQ=255)
                                                                                                  ccgccgGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGTGTTCGATCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGAGATTGCCGAAGCGGTGGc                                                     <  2:401187/150‑1 (MQ=255)
                                                                                                  ccgccgGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGc                                                     <  2:1155598/150‑1 (MQ=255)
                                                                                                    gccgGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGTCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGCTGGCac                                                   >  2:686325/1‑149 (MQ=255)
                                                                                                    gccgGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAg                                                   >  1:1004033/1‑150 (MQ=255)
                                                                                                             ttCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTc                                          >  1:477771/1‑150 (MQ=255)
                                                                                                             ttCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTc                                          <  1:494512/150‑1 (MQ=255)
                                                                                                              tCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTcc                                         <  2:20456/150‑1 (MQ=255)
                                                                                                               cGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCt                                        >  1:456732/1‑150 (MQ=255)
                                                                                                               cGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCt                                        >  1:929201/1‑150 (MQ=255)
                                                                                                                gTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTg                                       <  2:1004431/150‑1 (MQ=255)
                                                                                                                gTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTg                                       <  1:358901/150‑1 (MQ=255)
                                                                                                                gTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTg                                       <  1:82539/150‑1 (MQ=255)
                                                                                                                gTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTg                                       >  2:453667/1‑150 (MQ=255)
                                                                                                                  aaaTTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGAt                                     >  1:523124/1‑150 (MQ=255)
                                                                                                                       aaGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGACGCGGTGGCAGACTATATTCCTGATGGTTc                                >  1:1087237/1‑150 (MQ=255)
                                                                                                                       aaGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTc                                <  1:311650/150‑1 (MQ=255)
                                                                                                                            aGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTcaacaa                           >  1:526650/1‑150 (MQ=255)
                                                                                                                             gCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAat                          <  1:1357716/150‑1 (MQ=255)
                                                                                                                             gCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAat                          >  2:346572/1‑150 (MQ=255)
                                                                                                                             gCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAat                          >  2:1392834/1‑150 (MQ=255)
                                                                                                                              cAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAata                         <  1:429739/150‑1 (MQ=255)
                                                                                                                              cAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAata                         >  2:355290/1‑150 (MQ=255)
                                                                                                                                 ggCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATAtt                       <  1:394973/149‑1 (MQ=255)
                                                                                                                                  gCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTa                     >  1:693778/1‑150 (MQ=255)
                                                                                                                                    cTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATc                   >  2:1147333/1‑150 (MQ=255)
                                                                                                                                     tGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCa                  >  1:510387/1‑150 (MQ=255)
                                                                                                                                      gATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAg                                                                        >  2:743461/1‑95 (MQ=255)
                                                                                                                                      gATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTcaac                             >  1:1223237/1‑138 (MQ=255)
                                                                                                                                      gATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTc                                <  2:933583/135‑1 (MQ=255)
                                                                                                                                      gATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTc                                >  1:933583/1‑135 (MQ=255)
                                                                                                                                       gacccgcgCCACCGTGATGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGAAGGTTCAACAATATTTATCAcc                >  1:375068/5‑150 (MQ=255)
                                                                                                                                        acccgcgCCACCGTGATGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCa               >  1:993716/4‑150 (MQ=255)
                                                                                                                                         cccgcgCCACCGTGATGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCAt              >  1:735944/3‑150 (MQ=255)
                                                                                                                                          aCGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCAtt             <  2:974386/150‑1 (MQ=255)
                                                                                                                                           cgcgCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTg            <  2:318710/150‑1 (MQ=255)
                                                                                                                                           cgcgCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTg            >  1:1000664/1‑150 (MQ=255)
                                                                                                                                            gcgcCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTgg           <  2:1000664/150‑1 (MQ=255)
                                                                                                                                             cgcCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGt          >  1:1085971/1‑150 (MQ=255)
                                                                                                                                              gcCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTca                               <  2:1257294/128‑1 (MQ=255)
                                                                                                                                              gcCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTa         <  2:515731/150‑1 (MQ=255)
                                                                                                                                               ccATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTAc        <  1:565859/150‑1 (MQ=255)
                                                                                                                                                cATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACg       >  2:1380409/1‑150 (MQ=255)
                                                                                                                                                 aTCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGa      >  1:551419/1‑150 (MQ=255)
                                                                                                                                                 aTCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGa      >  1:1313279/1‑150 (MQ=255)
                                                                                                                                                  tCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGAc     <  1:845024/150‑1 (MQ=255)
                                                                                                                                                  tCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGAc     <  1:1070827/150‑1 (MQ=255)
                                                                                                                                                  tCACGGTGGCGCGGGTCGGGCTTCCAGCGACGTTAATACGGCCTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCCATTGCCGAAGCGGTCGCAGACTAAATTCCTGATGCTTCAACAATATTTATCACCATTGGTACGaa     >  2:279869/1‑149 (MQ=255)
                                                                                                                                                   cACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGACt    >  1:932557/1‑150 (MQ=255)
                                                                                                                                                     cGGTGGCGCGGGTCGGGCTTCAAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGACTGt  >  2:89396/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACCCGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGACTGT  >  NC_000913/4074693‑4074991

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