Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,776,138 G→T 6.3% G74C (GGC→TGC)  yfjX → CP4‑57 prophage; putative antirestriction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,776,1380GT 6.3% 191.6 / 2.4 61G74C (GGC→TGC) yfjXCP4‑57 prophage; putative antirestriction protein
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base T (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

TCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACC  >  NC_000913/2775848‑2776433
                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                       
ccccAACCGCGCTTTATTCCCTCATACGCCAGCCTTCGCGGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCCGCCATTGCTGAGCGtg                                                                                                                                                                                                                                                                                               <  1:601754/300‑2 (MQ=255)
ccccAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGtt                                                                                                                                                                                                                                                                                               <  1:86603/300‑1 (MQ=255)
       cgcgTTTTATTACCTCATACTCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGATACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCGGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCt                                                                                                                                                                                                                                                                                        <  1:470070/301‑1 (MQ=255)
              tattaCCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATg                                                                                                                                                                                                                                                                                 >  2:392003/1‑301 (MQ=255)
                    aaaaaacgccagcggcgatggcTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGACCCCACTTTCAGTTGAGCGTAACGCAGGTCCCCGATGACCATCATATCCGCTTCTGGCCGCCGCACTTTGCCACCATCCCGCCGTGGCTAACGCTGGAACCGCGCCTTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGCACGTTCACGGCCACTGATGCCGACg                                                                                                                                                                                                                                                                           >  1:404350/21‑301 (MQ=255)
                      actgtagttcagacgtgtgctcttccgatctATTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGt                                                                                                                                                                                                                                                                         <  2:336140/271‑1 (MQ=255)
                                                 tttgtgtCACGGGGTCATACCCATGACAACACAGACAAAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCAGCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGc                                                                                                                                                                                                                                              <  1:381223/295‑1 (MQ=255)
                                                  ttttttCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCa                                                                                                                                                                                                                                             <  2:451517/301‑1 (MQ=255)
                                                  tttATTCACGGAGACATACCCATGACAACACATACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCa                                                                                                                                                                                                                                             <  2:218022/301‑1 (MQ=255)
                                                    tatggtcgggtacgtACCCATGACTACACAGACACTGGACGAGCTCGTATCCGCTAAGTGGTCTGAAGTGGGGTTGAGCGTAACGCAGGTACCTGATGAACAGCATATCCTCGTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGGGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATg                                                                                                                                                                                                                                           <  1:77902/287‑1 (MQ=255)
                                                     aTTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCACGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGAAGCTGACGGTggatcggaagcggggagggtgggggaggagg                                                                                                                                                                                                                                          >  1:335735/1‑271 (MQ=255)
                                                      ttCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGaa                                                                                                                                                                                                                                         >  2:301536/1‑301 (MQ=255)
                                                        cACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAACGTTTATGGCTCCTGATGCTGAAGGTGACGATAAAAGGCAAGTGCCCAAAGGCATGaaat                                                                                                                                                                                                                                       >  1:266826/1‑299 (MQ=255)
                                                              gaCATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCCCAACGGCAAGAATGGCAATg                                                                                                                                                                                                                                 >  1:191201/1‑301 (MQ=255)
                                                                cATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAACTGGCATTTGATCAACGGCATGAATGGCAATgga                                                                                                                                                                                                                               >  2:493856/1‑300 (MQ=255)
                                                                    cccATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTCTCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGACCCTGATGCTGACGGTGACGATAAATGGCATTTGCCCAACGGCATGAATGGCAATGGAGCgg                                                                                                                                                                                                                           >  1:416176/1‑301 (MQ=255)
                                                                      cATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCCGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGaa                                                                                                                                                                                                                         >  2:606539/1‑301 (MQ=255)
                                                                                  gacacagTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGaa                                                                                                                                                                                                             <  2:191444/301‑1 (MQ=255)
                                                                                    cacagTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCACTGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCACTGGTGCGGAAATGAGAGCGGAAGc                                                                                                                                                                                                           >  1:533463/1‑301 (MQ=255)
                                                                                    cacagTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTACGGAAATGAGCGCGGAAGc                                                                                                                                                                                                           >  2:557481/1‑301 (MQ=255)
                                                                                       agTACGACATCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGAATATCCGCTTCTGGCCGCAGCACTTTGGCACAATCCCGCCGTGGATAACGAAGGAACCGCGCATTTTCGCATGGATGGACCGATTCTGTGATGAGTACTGCGGGGGTATCGGGTACTTTTACACGCGCAGACATGGCGGAGCGTTTATGGCTGCTGATGCTGCAGGTGACGATAAATGACATTAGCTCAACGGAATGAATGGCACAGGTGCAGAACGGAGAGCAGAAGcagg                                                                                                                                                                                                        >  1:43055/1‑298 (MQ=255)
                                                                                        ttatgACCTCGTATCCGCTAACGGATTCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGTATTGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCa                                                                                                                                                                                                       <  2:43108/297‑1 (MQ=255)
                                                                                        gTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCa                                                                                                                                                                                                       >  2:122713/1‑301 (MQ=255)
                                                                                         tACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAg                                                                                                                                                                                                      <  2:416739/301‑1 (MQ=255)
                                                                                                 cGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGACCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTAACGCt                                                                                                                                                                                              >  1:164326/1‑301 (MQ=255)
                                                                                                     cccGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGACCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCt                                                                                                                                                                                          <  2:164478/301‑1 (MQ=255)
                                                                                                               taTCCGAATTTGAGTTGTGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGa                                                                                                                                                                                <  2:267094/300‑1 (MQ=255)
                                                                                                                aTCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGaa                                                                                                                                                                               >  2:441686/1‑301 (MQ=255)
                                                                                                                   cGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGTCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAataa                                                                                                                                                                            >  2:84185/1‑300 (MQ=255)
                                                                                                                           ggttgagCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCa                                                                                                                                                                    <  2:195976/300‑1 (MQ=255)
                                                                                                                                         gaGGTACCCGATGAACATCATATCCGCTTTTGGTGTCAGCAGTTTGTCTCCATCCCGCAGTGGATATGGCTGGAACAGCGCATTTTCGCTTGGATTGACCGCTTCTGTGATGAGTACGGCGGGGGTATCTGGTCCTTTGACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAAGGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATGACGCCTGCCTCACCg                                                                                                                                                      <  1:391489/300‑1 (MQ=255)
                                                                                                                                           ggTACCCGATGAACAGCATATCCGCGTCTGGCCGTAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGaa                                                                                                                                                    <  1:301132/301‑1 (MQ=255)
                                                                                                                                                 cGATTAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGAt                                                                                                                                              <  2:534136/301‑1 (MQ=255)
                                                                                                                                                      aaCAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCaa                                                                                                                                         >  2:347798/1‑300 (MQ=255)
                                                                                                                                                         aGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTTGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAc                                                                                                                                      <  2:421498/301‑1 (MQ=255)
                                                                                                                                                           cATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTNNCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCAGCAGGCATCGCTGTCTGCCTGCTTGAATAACACCAACACCCCCGCCCCCCCACATCCGATGcactaaccg                                                                                                                                    >  1:32166/1‑293 (MQ=255)
                                                                                                                                                                              agcaCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGTGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCg                                                                                                                 <  1:84061/301‑1 (MQ=255)
                                                                                                                                                                                    ttGGCACCATCCCGCAGTGGATAACGCTGGAACCGTGCTTTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTa                                                                                                           <  1:441122/301‑1 (MQ=255)
                                                                                                                                                                                       agccgtgtgatcttccgatctcGCTGGAACCGCGCATTTTGGCTGGGAGGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAAGGGCATGTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGc                                                                                                        <  2:238288/280‑1 (MQ=255)
                                                                                                                                                                                            agctcGCTGTGGATAACGCTGGAACCGCGCATTTTCGCCGGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCGTTTTACACGCGCAGCAATGGATGAGCGTTTATGGCTCCGGATGCTGAAGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCGGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTgc                                                                                                   <  1:493210/297‑1 (MQ=255)
                                                                                                                                                                                                cgtggggtgtatggCTGGAACCGTGCATTTTCGCTTTGATGGACCGTTTCTGTGAGGAGTACTGGGGGGGTTTCTGGTCCTTTTACACGCTAAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTgcagca                                                                                               <  1:347335/288‑1 (MQ=255)
                                                                                                                                                                                                      gggtAACGTTGGAGCTGGGCTTTTTCGCTTGGATGGACCGGTTCTTGGAGGAGTACTTGGGGGTGATCTGGTCAGTTTACACGCTCAGCAGTGGCGGAGCTTGTGTGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAGCGGCAGGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGAGTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGa                                                                                         <  1:545854/298‑1 (MQ=255)
                                                                                                                                                                                                          aaCGCTTGAATCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGAGGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTAGTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATcc                                                                                     <  2:214215/301‑1 (MQ=255)
                                                                                                                                                                                                            cGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCACTGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGCATGCGCTGCCATGACGTAGCCATATTACCGCCTGCGGGATTAcacagatcggaaacgccgcgcggt                                                                                   >  1:238021/1‑278 (MQ=255)
                                                                                                                                                                                                            cGCTGGAACCGCGCATTTTCGCCTGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAACACTATTAACGCCTGCGGGATTACGCCCTGCAACACCCTCAATCCAgc                                                                                  >  1:229365/1‑301 (MQ=255)
                                                                                                                                                                                                                                 gcctgatggaCCGCTTCTTTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGa                                                              <  2:229608/297‑1 (MQ=255)
                                                                                                                                                                                                                                  ctggctggACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGCCGATAAATGGCATGTGCTCAACGGCATGAAGGGCAATGGTACGGAAATCAGCGCCGAGGCCGCCGGTATCGCTGCTTGCCTGATTGCATATAGCCATCACACCTGCCTCCCCGAATGCGAGTCCATGACGGAAGACTAATCACCCCCGCGGGAATACGCTCTACCGCCCCCTAAAACCCAAGCCATTACGCCCaatactcat                                                             >  1:38686/1‑296 (MQ=255)
                                                                                                                                                                                                                                   ggtatggaCCGTTTTTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTGAGTAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGAtt                                                            <  1:122589/298‑1 (MQ=255)
                                                                                                                                                                                                                                         gaCCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCCGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGaa                                                      >  2:462278/1‑300 (MQ=255)
                                                                                                                                                                                                                                                tttgtgttgggttcggggtgggTATCTGGTCTTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCatg                                               <  1:605778/282‑1 (MQ=255)
                                                                                                                                                                                                                                                   acgcttttccggtctggggTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCatgatg                                            <  1:341540/286‑1 (MQ=255)
                                                                                                                                                                                                                                                       tgatgaCTGCGTGGGTTGTTGTTTCTTTGTCACGCTCAGCAATGGCGGTGCGTTTATGGCTCCTGATGCTGACGGGGACGATGAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCAGCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAAc                                        <  2:337281/296‑1 (MQ=255)
                                                                                                                                                                                                                                                         ggTACTGCGGGGGTATCTGGTCTTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAg                                      <  2:32217/300‑1 (MQ=255)
                                                                                                                                                                                                                                                          gTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAga                                     >  2:215000/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                  ggggTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGaaatcggaagagcaa                             >  2:342003/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                      tATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAAAAGTCACTCAacacac                         >  1:203535/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                           ggTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCCGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTCCGCACTGCAGCACCCTGAATCCAGCCACATTCAGCGCATTATTGATTAAGGATACTCATGATGGAACaataccccaacccccaaacc                    >  1:605942/1‑282 (MQ=255)
                                                                                                                                                                                                                                                                               gccataaacGCTCCGCCATTGCTGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCCTGACGGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCCCCCTGAATCCAGCGCCCTTATGCTCCATACTGATTAAATAAACTCACGACGGACAAATCACCCAcccccccacccccac                >  1:4007/8‑290 (MQ=255)
                                                                                                                                                                                                                                                                                     aacGCTCCGCCATTGCTGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGATACTTc          >  2:179299/2‑301 (MQ=255)
                                                                                                                                                                                                                                                                                         cTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTAACGCCTGCGGGATTCCGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACCCATGATGGAACAATCACTAATCCAACAGACACCGGTACTTCCCCt      >  1:457107/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                             gCAATTGCGTAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGAcc  <  2:404913/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                       
TCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACC  >  NC_000913/2775848‑2776433

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: