Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,891,085 G→A P272S (CCG→TCG)  cysJ ← sulfite reductase, flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,891,0850GA98.2% 178.8 / NA 56P272S (CCG→TCG) cysJsulfite reductase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/29);  minor base T (0/1);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

ATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCA  >  NC_000913/2890840‑2891323
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              
aTGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCGGATACCAGACGCCCAGCGCGTCACCCGACTGGt                                                                                                                                                                                                                                           <  1:121507/251‑1 (MQ=255)
         cgtcgatttagcggccaccaGCGGCAGCTGTGTTTCACGGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGTCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACTCATGc                                                                                                                                                                                                                                  <  2:290601/237‑1 (MQ=255)
             gCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGGGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTGCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTGCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGa                                                                                                                                                                                                                              <  2:201303/251‑1 (MQ=255)
                 tATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTcc                                                                                                                                                                                                                          >  1:78086/1‑250 (MQ=255)
                        caccaGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAgg                                                                                                                                                                                                                   >  2:133298/1‑250 (MQ=255)
                         accaGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTc                                                                                                                                                                                                                  >  2:109261/1‑251 (MQ=255)
                          ccaGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCa                                                                                                                                                                                                                 <  2:78086/251‑1 (MQ=255)
                           caGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCaa                                                                                                                                                                                                                <  1:348454/251‑1 (MQ=255)
                             cgggactcagtgtttgttggCGGGTAGGCGTGGCAGAAGTCTCTGCAATGTTGGCGGTGTGGACAGTCAGTTCGAAGGGCCACTGTATGGCTCCTGTCAGATCCACCGGTTTGCCCCCGACGGTGCCAGGTTCAGCGCCTTTCCGCCACAGCACTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCACAGCGAGTCACCCGACTGGTAAAGCAGGCCCTAGTCACCTAAGTCAAtc                                                                                                                                                                                                              <  2:188292/232‑2 (MQ=255)
                                      gtgtTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATgg                                                                                                                                                                                                     <  1:150939/251‑1 (MQ=255)
                                      gtgtTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCAACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATgg                                                                                                                                                                                                     >  1:313884/1‑251 (MQ=255)
                                          ttCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCCGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGGTCCACGCCTTTCAGCCCCAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGTTCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCCTGCCCGAGTCACCTAAGTCAATTTCGATATGGCGa                                                                                                                                                                                                  >  2:6845/1‑250 (MQ=255)
                                             ccTGCGGGTAAGCGTGGCGTAATTCCCAACAATGTTGGCGGTGTGGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGTGAACGCAGGCCCGAGGCACCTAAGTCAATTTCGAGATGCCGAACGt                                                                                                                                                                                              <  2:43944/250‑1 (MQ=255)
                                              cTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGATCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTc                                                                                                                                                                                             >  1:197892/1‑251 (MQ=255)
                                                gCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCtt                                                                                                                                                                                           <  1:259246/251‑1 (MQ=255)
                                                gCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCt                                                                                                                                                                                            <  2:125657/250‑1 (MQ=255)
                                                    gTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGGGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATCCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCtttttt                                                                                                                                                                                       >  2:212367/1‑250 (MQ=255)
                                                      aaGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTTACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGTCGAACGTCTTTTTCag                                                                                                                                                                                     <  1:20734/251‑1 (MQ=255)
                                                       aGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGGCTTTTTc                                                                                                                                                                                       >  1:361602/1‑248 (MQ=255)
                                                       aGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGCCAATTTCGATATGGCGAACGCCTTTTTc                                                                                                                                                                                       >  2:269950/1‑248 (MQ=255)
                                                                       caacaaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCGCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTGTTCAGAGTTACGCCCGGTAAtt                                                                                                                                                                    <  2:225628/251‑1 (MQ=255)
                                                                            aTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTg                                                                                                                                                               <  2:304215/251‑1 (MQ=255)
                                                                              gTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGt                                                                                                                                                             <  1:220044/251‑1 (MQ=255)
                                                                              gTTGGCGGTGTTGACGGTCAGCTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGCCAGGCTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCCGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTACGTCAATTTCGATATGGCGAACGTCTTTTTCCGAGTTACGCCCGGTAATTTTCTGGt                                                                                                                                                             >  2:161305/1‑251 (MQ=255)
                                                                              gTTGGCGGAGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACCGGTTCATCGCCTTTCCGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGt                                                                                                                                                             >  2:294641/1‑251 (MQ=255)
                                                                                       gtTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACTCCCAGCGCGTCACCCGTCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACagaga                                                                                                                                                    <  1:346801/251‑1 (MQ=255)
                                                                                             ggTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAg                                                                                                                                              <  1:2519/251‑1 (MQ=255)
                                                                                              gTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGatat                                                                                                                                                                                                       <  1:69740/193‑1 (MQ=255)
                                                                                              gTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGatat                                                                                                                                                                                                       >  2:69740/1‑193 (MQ=255)
                                                                                                 aGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGccac                                                                                                                                          >  1:171658/1‑251 (MQ=255)
                                                                                                      gAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAgcg                                                                                                                                     >  1:314555/1‑251 (MQ=255)
                                                                                                                 gTAGCGCTTCGGTCAGAGGCAACGTTTTGCCATCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCtt                                                                                                                            <  2:351648/249‑1 (MQ=255)
                                                                                                                   aGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCt                                                                                                                        <  1:42958/251‑1 (MQ=255)
                                                                                                                                      aCGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGAtt                                                                                                     <  1:354303/251‑1 (MQ=255)
                                                                                                                                           ttGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCGGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCGGGTGTGGATTTCAt                                                                                                 <  2:314555/250‑1 (MQ=255)
                                                                                                                                              cccTCGACGGTGACAGGTTCATCGCCTTTCCGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATTCCCGAGGCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTATCCACCAGCGGCGCGTCTTTGCTGTTCGGGCTGGTGTGGATTTCATTTAc                                                                                             <  2:221505/251‑1 (MQ=255)
                                                                                                                                              cccTCGACGGTGACAGGTTCATCGCCTTTCAGCCCCAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATCCCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTAc                                                                                             <  1:77325/251‑1 (MQ=255)
                                                                                                                                                cTCGACGGTGACAGGGTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGCACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCAGGGGTGGATTTCATTTAccc                                                                                           >  2:230014/1‑250 (MQ=255)
                                                                                                                                                         gACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGt                                                                                    >  2:272521/1‑249 (MQ=255)
                                                                                                                                                                 cATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGGGTGGATTTCATTTACCCCGCCAGTAGCGACGGAt                                                                          >  2:121972/1‑251 (MQ=255)
                                                                                                                                                                      ccTTTCAGCCACAGCAGTTCGCCAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGGGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCgc                                                                     >  2:259728/1‑250 (MQ=255)
                                                                                                                                                                         ttCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAgg                                                                  >  1:9526/1‑251 (MQ=255)
                                                                                                                                                                          tCAGCCACAGCAGTTAGACAAGTTCTTTCACCAGTGCAGGATCGTTCAGATCCCAGACGCCCAGCACCTCACCCGACTGGCAACGCATGCCCGAGTCACCTCAGTCAATTTCGATCCGGCGAACGTCTTTTTCAGAGGTACGCCCGGTACTTTTCTGGTTAACAGAGCGGATAGCAAACAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCTGTAGCGACGGATTGCGAAGGt                                                                 >  1:364406/1‑251 (MQ=255)
                                                                                                                                                                           cAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTg                                                                >  1:314180/1‑251 (MQ=255)
                                                                                                                                                                                acaGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAgcggc                                                                                                                                   >  2:287482/1‑179 (MQ=255)
                                                                                                                                                                                acaGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAgcggc                                                                                                                                   <  1:287482/179‑1 (MQ=255)
                                                                                                                                                                                       ttCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGc                                                               <  1:162970/240‑1 (MQ=255)
                                                                                                                                                                                       ttCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGc                                                               >  2:162970/1‑240 (MQ=255)
                                                                                                                                                                                                   tttCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGGGCCGCGACAGGCGCACGCGGTTTa                                         >  2:44340/1‑250 (MQ=255)
                                                                                                                                                                                                      cACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAAGTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAgcg                                     <  1:294641/251‑1 (MQ=255)
                                                                                                                                                                                                       aCCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCATCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAgcgc                                    >  1:65754/1‑251 (MQ=255)
                                                                                                                                                                                                         catggCCGGAACGTTCGGAGTCCAGACGCCCAGCGTGTATCCCTACTTGTAACGCATGACCGGGTTACCTAAGTCAATTTCGATATGGCGAACGTCTTTTGCAGAGTTACGCCCGTTAAGTTTCTGGTTAGCAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCGGTACGGGGTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCCCGCGATTTAAGCGCAt                                  <  2:278000/247‑1 (MQ=255)
                                                                                                                                                                                                                 gATCGTTCTGATACAAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGCGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTCGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCCGTAGCGACGGATTGCGAAGGTGCCCCGACATGCGCACGCGATTTAAGCGCCTCAACCAcg                          >  2:114906/1‑251 (MQ=255)
                                                                                                                                                                                                                               cAGACGCCCAGCGCGTCCCCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTc            <  2:9526/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTg        <  1:269950/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gcgcGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCTTCTTTGCGGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGAGTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGgcagca  <  1:37976/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              
ATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCA  >  NC_000913/2890840‑2891323

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