Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,980,118 Δ1 bp 7.0% coding (366/429 nt) uspC → universal stress protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,980,1180G. 7.0% 277.2 / 14.3 86coding (366/429 nt)uspCuniversal stress protein C
Reads supporting (aligned to +/- strand):  ref base G (39/40);  new base . (2/4);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

TTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCG  >  NC_000913/1979970‑1980266
                                                                                                                                                     |                                                                                                                                                    
ttCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCgg                                                                                                                                                     >  2:85045/1‑150 (MQ=255)
 tCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCCGAGCAGGATGCTCgc                                                                                                                                                     >  2:1593335/1‑149 (MQ=255)
 tCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGc                                                                                                                                                    >  1:1114010/1‑150 (MQ=255)
 tCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGc                                                                                                                                                    >  2:1503071/1‑150 (MQ=255)
   aGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCCGAGCAGGATGCTCgcga                                                                                                                                                   <  2:297650/150‑4 (MQ=255)
    ggATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCCGAGCAGGATGCTCgcgaa                                                                                                                                                  <  2:379822/150‑5 (MQ=255)
    ggATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGaa                                                                                                                                                 <  2:141400/150‑1 (MQ=255)
      aTGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCTGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGaaaa                                                                                                                                               <  2:720818/150‑1 (MQ=255)
      aTGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCCGAGCAGGATGCTCGCGaaaa                                                                                                                                               <  2:1020791/150‑1 (MQ=255)
      aTGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGaaaa                                                                                                                                               >  1:1392907/1‑150 (MQ=255)
       tGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAag                                                                                                                                              >  2:1229464/1‑150 (MQ=255)
         cGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCCGAGCAGGATGCTCGCGAAAagaa                                                                                                                                            <  1:274348/150‑2 (MQ=255)
         cGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCCGAGCAGGATGCTCGCGAAAagaa                                                                                                                                            <  1:1593335/150‑2 (MQ=255)
          ggggTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGt                                                                                                                                           <  2:1214926/150‑1 (MQ=255)
           gtgTATCCCGTTGACAAGACGTTTATTGACTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGtt                                                                                                                                          >  2:102388/3‑150 (MQ=255)
               aTCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTg                                                                                                                                      >  1:79025/1‑150 (MQ=255)
               aTCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTg                                                                                                                                      <  2:807015/150‑1 (MQ=255)
                 cccGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGcc                                                                                                                                    >  1:1203836/1‑150 (MQ=255)
                         cAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGa                                                                                                                            <  1:1504741/150‑1 (MQ=255)
                           aGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAgg                                                                                                                          <  1:766209/150‑1 (MQ=255)
                           aGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAgg                                                                                                                          <  1:771950/150‑1 (MQ=255)
                           aGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAgg                                                                                                                          <  2:271239/150‑1 (MQ=255)
                                tttATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTACGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGAt                                                                                                                     <  1:395711/150‑1 (MQ=255)
                                 ttATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGAtg                                                                                                                    <  1:714919/150‑1 (MQ=255)
                                   aTTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATTATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGAtgtg                                                                                                                  <  2:1531476/150‑1 (MQ=255)
                                   aTTGCCTACGGTGAGTTAAGCGAACATATTCTAGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGAtgtg                                                                                                                  <  1:942483/150‑1 (MQ=255)
                                    ttGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCAACCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGc                                                                                                                 <  1:521923/150‑1 (MQ=255)
                                      gCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTg                                                                                                               >  1:1439689/1‑150 (MQ=255)
                                       ccTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGt                                                                                                              <  2:248226/150‑1 (MQ=255)
                                        cTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGtt                                                                                                             <  2:1392907/150‑1 (MQ=255)
                                         tACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCTGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTa                                                                                                            <  2:1483360/150‑1 (MQ=255)
                                            ggTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGtt                                                                                                         <  1:691587/150‑1 (MQ=255)
                                                aGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCAc                                                                                                     <  1:85045/150‑1 (MQ=255)
                                                  ttAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACtt                                                                                                   >  1:270055/1‑150 (MQ=255)
                                                    aaGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTAc                                                                                                 >  1:1511854/1‑150 (MQ=255)
                                                       cGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACggg                                                                                              >  1:173245/1‑150 (MQ=255)
                                                       cGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACggg                                                                                              >  2:665713/1‑150 (MQ=255)
                                                        gAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGga                                                                                             <  2:1439689/150‑1 (MQ=255)
                                                        gAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGga                                                                                             >  2:1201231/1‑150 (MQ=255)
                                                            atatTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGAtt                                                                                         >  2:308818/1‑150 (MQ=255)
                                                             tatTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTa                                                                                        >  2:1629110/1‑150 (MQ=255)
                                                              atTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTaa                                                                                       >  2:295038/1‑150 (MQ=255)
                                                                tCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAATAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAAcc                                                                                     <  2:580162/150‑1 (MQ=255)
                                                                tCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAAcc                                                                                     >  2:989608/1‑150 (MQ=255)
                                                                tCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAAcc                                                                                     <  1:615927/150‑1 (MQ=255)
                                                                 cTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCg                                                                                    <  1:230065/150‑1 (MQ=255)
                                                                  tGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGc                                                                                   >  2:95670/1‑150 (MQ=255)
                                                                    gAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTTCTGTTAGTTCCACTTACGGGAGATTAACCGCTc                                                                                 >  2:294506/1‑150 (MQ=255)
                                                                      aGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCt                                                                               <  2:1064946/150‑1 (MQ=255)
                                                                           gTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCa                                                                          >  1:1312177/1‑150 (MQ=255)
                                                                           gTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCa                                                                          >  1:1306439/1‑150 (MQ=255)
                                                                            tCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCaa                                                                         <  1:1503071/150‑1 (MQ=255)
                                                                              acaAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGc                                                                       >  1:805580/1‑150 (MQ=255)
                                                                                aaGCATCATTTCGATTTGGTTATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGATCATCATTCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCtt                                                                     >  2:665377/1‑150 (MQ=255)
                                                                                  acatcatTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTg                                                                   >  2:1533027/2‑150 (MQ=255)
                                                                                   catcatTTCGATTTGGTGATTTGCGGTAATCACAATAACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTgg                                                                  >  1:10323/1‑150 (MQ=255)
                                                                                     tcatTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGaa                                                                <  1:1189655/150‑1 (MQ=255)
                                                                                      catTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGaaa                                                               >  1:151863/1‑150 (MQ=255)
                                                                                       atTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAg                                                              >  1:1050042/1‑150 (MQ=255)
                                                                                        tttCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTTCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTCCGCAAGCTTTGGAAAgg                                                             <  1:294506/150‑1 (MQ=255)
                                                                                        tttCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAgg                                                             <  2:1314233/150‑1 (MQ=255)
                                                                                         ttCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGt                                                            >  2:1272995/1‑150 (MQ=255)
                                                                                              tttGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCaa                                                       <  1:1028558/150‑1 (MQ=255)
                                                                                                 ggTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACtt                                                    >  2:543143/1‑150 (MQ=255)
                                                                                                    gATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTAt                                                 >  1:257545/1‑150 (MQ=255)
                                                                                                      tttGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATcg                                               <  1:1229464/150‑1 (MQ=255)
                                                                                                         gCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCgcgc                                            >  2:656499/1‑150 (MQ=255)
                                                                                                            gTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgctgc                                         <  2:1114010/150‑1 (MQ=255)
                                                                                                             tAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgctgct                                        >  2:1262170/1‑150 (MQ=255)
                                                                                                                 cacaatcacaGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCt                                    <  2:532453/150‑1 (MQ=255)
                                                                                                                 cacaatcacaGTTTCTTTTCGCGAGCATACTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCAAATTACTGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCt                                    >  2:1187761/1‑150 (MQ=255)
                                                                                                                      tcacaGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcg                               <  1:546311/150‑1 (MQ=255)
                                                                                                                      tcacaGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTAACGCGCTGCTGGCTTTCcg                               >  1:1550537/1‑150 (MQ=255)
                                                                                                                       cagaGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAATAGTTATTTCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcgc                              <  2:493476/147‑1 (MQ=255)
                                                                                                                        acaGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcgcg                             <  1:857619/150‑1 (MQ=255)
                                                                                                                           gTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCtt                          >  1:1498376/1‑150 (MQ=255)
                                                                                                                             ctCTTTTCGCCGAGCAGGATGCTCGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTcg                        >  1:622745/2‑150 (MQ=255)
                                                                                                                               cTTTTCGCCGAGCAGGATGCTCGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTcgcg                      >  1:1162768/1‑150 (MQ=255)
                                                                                                                               cTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTcgcg                      <  1:989608/150‑1 (MQ=255)
                                                                                                                                    cgcgAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGAAGCGCGGaaag                 <  1:1262170/150‑3 (MQ=255)
                                                                                                                                     gcgAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTa                >  1:91597/1‑150 (MQ=255)
                                                                                                                                           caggatGCTCGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCt          <  2:1207063/145‑1 (MQ=255)
                                                                                                                                           cATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGc           <  2:977968/150‑1 (MQ=255)
                                                                                                                                             tCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCtt         >  1:410239/1‑150 (MQ=255)
                                                                                                                                              ccTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCttt        <  2:1203836/150‑1 (MQ=255)
                                                                                                                                              ccTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCttt        <  1:450673/150‑1 (MQ=255)
                                                                                                                                                 gCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAgg     <  2:1424230/150‑1 (MQ=255)
                                                                                                                                                 gCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAgg     >  1:561292/1‑150 (MQ=255)
                                                                                                                                                  cTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGt    <  1:665377/150‑1 (MQ=255)
                                                                                                                                                    cGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCg  >  2:394777/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCG  >  NC_000913/1979970‑1980266

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: