Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,532,279 G→A 11.1% intergenic (+55/‑130) cysZ → / → cysK sulfate transporter, sulfite inhibited/cysteine synthase A, O‑acetylserine sulfhydrolase A subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,532,2790GA11.1% 252.1 / 18.3 90intergenic (+55/‑130)cysZ/cysKsulfate transporter, sulfite inhibited/cysteine synthase A, O‑acetylserine sulfhydrolase A subunit
Reads supporting (aligned to +/- strand):  ref base G (37/42);  new base A (5/5);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

TATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCAAAATCACCAATATGCAGTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCC  >  NC_000913/2532009‑2532549
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                              
agacaggATTACCCCGCCTATAACCAAANNNTGCCGGTTAAAGAGGTGCGCACCGCCCTGCGCACACGAAAAATCACCAGTATGCAGTTTGGTGTTTTAACCGGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGccc                                                                                                                                                                                                                                                                             <  3:139537/268‑1 (MQ=255)
  ttGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCAACATCACCAATATGCAGTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGGTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCtt                                                                                                                                                                                                                                                                          >  2:108095/1‑275 (MQ=255)
          aCCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCAAAATCACCAATATGCATTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATATGGGTCGATTGCTATCGCGATAAAAACGAGATGGTGCGGTAAACATCTACCGGTCATTTGTTAAACCGTTTGTTTGAAACATGGGTGGCTTGTTTCACCACTTAttaccta                                                                                                                                                                                                                                                                  >  4:268348/1‑270 (MQ=255)
          aCCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCAAAATCACCAATATGCAGTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATc                                                                                                                                                                                                                                                                  >  2:299541/1‑275 (MQ=255)
          aCCCCTTCGATAACCACAAAGCGCCGTTTACAGAGATGCGCACCGCCCTTCGCACCCGCAAAATCACCAATATGCAGTTTGGTGCTTTAACACGCCAGTATACGCGGATCCCGCTGCTTAATCTATTCNTCATGCCCGTTCCCGTTTGTGGCGCGACGGATCTTTTGGTCGATCGCTGTCGCGATAAACACGAGATGCNGCTGTAACACTCTACCGGTGATTTTGTAAACGGTTTGTTTGGCACNANGTGTGCTTCTACCTCCCATTATTACATc                                                                                                                                                                                                                                                                  >  3:46827/1‑275 (MQ=255)
                tCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCAAAATCACCAATATGCAGTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCAt                                                                                                                                                                                                                                                            <  2:128889/275‑1 (MQ=255)
                            acaggCCGTTTAAAGAAATGCAACGCGACCCTCGTACACGCAAAATTACCAGTATGCAGTTTGGTGCTTCAACCAGCCGGTTTACGATGATCCCGCGGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACTGCGATGTGGGTCTATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTcc                                                                                                                                                                                                                                                <  4:482959/271‑1 (MQ=255)
                            aaGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCAAAATCACCAATATGCAGTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTcc                                                                                                                                                                                                                                                <  1:69316/275‑1 (MQ=255)
                             aGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCAAAATCACCAATAGGCAGTTTGGGGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTccc                                                                                                                                                                                                                                               <  2:59649/275‑1 (MQ=255)
                               taccatttAAAGAGATGCGCACCGCCCTGCGCACACGCAAAATCACCAATATGCAGTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCtt                                                                                                                                                                                                                                             <  2:496577/270‑1 (MQ=255)
                               tGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCAAAATCACCAATATGCAGTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCtt                                                                                                                                                                                                                                             <  1:43751/275‑1 (MQ=255)
                                 ccGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCAAAATCACTAATATGCAGTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCtt                                                                                                                                                                                                                                             >  2:27109/1‑273 (MQ=255)
                                              tGCGCACCGCCCTGCGCACACGCAAAATCACCAATATGCAGTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTc                                                                                                                                                                                                                                            >  2:65074/1‑261 (MQ=255)
                                                       cccTGCGCACACGCAAAATCACCAATATGCAGTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgccgcc                                                                                                                                                                                                                                                                              >  2:278550/1‑218 (MQ=255)
                                                            cgcACACGCAAAATCACCAATATGCAGTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTGATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACt                                                                                                                                                                                                                 <  2:105687/274‑1 (MQ=255)
                                                                     aaaaTCATCAATATGCAGTTTGTTGGTCTAACCAGACGGTTTAGGGTGATCCCGATGCTTAATCTGTTCAACATGCCCGTTTCCGTTTGTGGCGCGACGGCGAGTTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATAtt                                                                                                                                                                                                       <  4:102820/275‑1 (MQ=255)
                                                                               atatGCAGTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGNGGCGGTAACAATCTACCGGGTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTNTTCCATCTTGCATGTCATTATTTTCCTTCTTTATATAGATATGCTNANAGCTTAATTCCTCATAtttttttatcgg                                                                                                                                                                                             >  3:45685/1‑265 (MQ=255)
                                                                                       tttGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTCCATGTCattat                                                                                                                                                                                                                                                    >  1:144705/1‑212 (MQ=255)
                                                                                        ttGGTTCTTTAACCAGCCTGTTTACGATGATCCCGCGGCTTAATCGGTTCCTCATGCCCGTTGCCGTTGGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTAAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTAc                                                                                                                                                                                    <  2:341212/275‑1 (MQ=255)
                                                                                              cTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCtgttg                                                                                                                                                                              >  2:197631/1‑275 (MQ=255)
                                                                                                    ccAGCCTGTTTACGATGATCCCGCGGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATcc                                                                                                                                                                         <  2:219730/274‑1 (MQ=255)
                                                                                                        ccTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATTCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGGTGTGGCGGTGACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGTGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTTTATATAGCTATGCTACATCCTTCCTTCCTCATATTCTCTGAGCGGGTATTCTACCTGTTTTATCCCCAtt                                                                                                                                                                    >  1:121369/1‑275 (MQ=255)
                                                                                                        ccTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCattat                                                                                                                                                                                                                                                    >  2:528579/1‑195 (MQ=255)
                                                                                                        ccTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTc                                                                                                                                                                                                               >  2:190046/1‑232 (MQ=255)
                                                                                                            tttACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTcc                                                                                                                                                                                                                                                <  4:528574/195‑1 (MQ=255)
                                                                                                            tttACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCa                                                                                                                                                                                                           <  4:190044/232‑1 (MQ=255)
                                                                                                             ttACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATa                                                                                                                                                               <  1:45685/275‑1 (MQ=255)
                                                                                                             ttACGATGAGCCCGCTGCTTAATCTGTTCATCATGTCCGGTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATa                                                                                                                                                               <  1:46827/275‑1 (MQ=255)
                                                                                                                       cccGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGTTGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGa                                                                                                                                                     <  2:167945/275‑1 (MQ=255)
                                                                                                                          gctgctTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgccg                                                                                                                                                                                                                                                                                >  6:56120/1‑149 (MQ=255)
                                                                                                                           ctgctTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATACACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTTAGCGGGTATCCTCCCTGTTGTATACCAATTTCCTACAGTTAAGTCCAgg                                                                                                                                                 >  4:207191/1‑275 (MQ=255)
                                                                                                                               ttAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCt                                                                                                                                                                                                                                                                           <  6:159273/149‑1 (MQ=255)
                                                                                                                                 aaTCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCtgt                                                                                                                                                                                                                                                                         <  6:223000/149‑3 (MQ=255)
                                                                                                                                 aaTCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCtgt                                                                                                                                                                                                                                                                         <  6:414085/149‑3 (MQ=255)
                                                                                                                                          atcatGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCtt                                                                                                                                                                                                                                                                >  6:161237/1‑149 (MQ=255)
                                                                                                                                          atcatGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCtt                                                                                                                                                                                                                                                                >  6:47143/1‑149 (MQ=255)
                                                                                                                                                 ccGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTc                                                                                                                                                                                                                                                         <  6:192687/149‑1 (MQ=255)
                                                                                                                                                    ttGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGNCCCTTATTCCATCTTGCATGTCANTATTNc                                                                                                                                                                                                                                                 >  4:507232/1‑154 (MQ=255)
                                                                                                                                                    ttGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCatt                                                                                                                                                                                                                                                      <  5:56120/149‑1 (MQ=255)
                                                                                                                                                    ttGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGtatat                                                                                                                                                                                                                                     >  4:2519/1‑166 (MQ=255)
                                                                                                                                                        cGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTAttt                                                                                                                                                                                                                                                  <  5:234951/149‑1 (MQ=255)
                                                                                                                                                        cGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCtt                                                                                                                                                                                                                                             <  2:507237/154‑1 (MQ=255)
                                                                                                                                                        cGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGat                                                                                                                                                                                                                                 <  2:2519/166‑1 (MQ=255)
                                                                                                                                                            tgtgGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTTTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCNTGCATGTCNTTAttta                                                                                                                                                                                                                                                 >  3:30387/1‑145 (MQ=255)
                                                                                                                                                                gcgcGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCtt                                                                                                                                                                                                                                             <  1:30387/146‑1 (MQ=255)
                                                                                                                                                                gcgcGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTg                                                                                                                                                                                                                                          <  7:367312/149‑1 (MQ=255)
                                                                                                                                                                 cgcgACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCt                                                                                                                                                                                                                                                                           >  8:265214/1‑115 (MQ=255)
                                                                                                                                                                 cgcgACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCt                                                                                                                                                                                                                                                                           <  7:265214/115‑1 (MQ=255)
                                                                                                                                                                     aCGGCGATGTTGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTGACAATCTACCGGTTATTTTGTAAACGGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGtatat                                                                                                                                                                                                                                     >  6:432521/1‑149 (MQ=255)
                                                                                                                                                                                 gTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATg                                                                                                                                                                                                                              >  7:472008/1‑144 (MQ=255)
                                                                                                                                                                                 gTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATg                                                                                                                                                                                                                              <  8:472008/144‑1 (MQ=255)
                                                                                                                                                                                    gATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTTAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCtt                                                                                                                                                                                                                                             <  2:272787/126‑1 (MQ=255)
                                                                                                                                                                                     aTTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCt                                                                                                                                                                                                                     >  8:173055/1‑149 (MQ=255)
                                                                                                                                                                                            tCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGtata                                                                                                                                                                                                                                      >  2:435799/1‑125 (MQ=255)
                                                                                                                                                                                                gATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGa                                                                                                                                                                                                                                  <  4:435794/125‑1 (MQ=255)
                                                                                                                                                                                                gATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCat                                                                                                                                                                                                          >  7:228304/1‑149 (MQ=255)
                                                                                                                                                                                                gATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCat                                                                                                                                                                                                          >  7:57681/1‑149 (MQ=255)
                                                                                                                                                                                                       acGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGata                                                                                                                                                                                                                                >  7:325507/1‑120 (MQ=255)
                                                                                                                                                                                                       acGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGata                                                                                                                                                                                                                                <  8:325507/120‑1 (MQ=255)
                                                                                                                                                                                                       acGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTctct                                                                                                                                                                                                   >  5:332996/1‑149 (MQ=255)
                                                                                                                                                                                                         gcgATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGa                                                                                                                                                                                                 <  7:238292/149‑1 (MQ=255)
                                                                                                                                                                                                             tgtgGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCgg                                                                                                                                                                                             <  5:47143/149‑1 (MQ=255)
                                                                                                                                                                                                             tgtgGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCgg                                                                                                                                                                                             <  5:161237/149‑1 (MQ=255)
                                                                                                                                                                                                                             tACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCtgttgt                                                                                                                                                                             >  5:143592/1‑149 (MQ=255)
                                                                                                                                                                                                                                tggTTATTTGGTAACACCTTTGTGTGAAACAGGGGTGGTTGTTGCCGCCCCTTATTCCATTTTGCATGTCATTATTTCCCTGCTATATCTAGATATGCTAAATCCTGACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGa                                             <  4:299537/273‑1 (MQ=255)
                                                                                                                                                                                                                                      ttttGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGccc                                                                                                                                                                                                                                                                             >  5:331172/1‑44 (MQ=255)
                                                                                                                                                                                                                                      ttttGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGccc                                                                                                                                                                                                                                                                             <  6:331172/44‑1 (MQ=255)
                                                                                                                                                                                                                                      ttttGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAAtt                                                                                                                                                                    >  7:266921/1‑149 (MQ=255)
                                                                                                                                                                                                                                         tGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCa                                                                                                                                                                 <  6:381150/149‑1 (MQ=255)
                                                                                                                                                                                                                                             aaCCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACa                                                                                                                                                             >  8:89138/1‑149 (MQ=255)
                                                                                                                                                                                                                                              aCCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAg                                                                                                                                                            >  7:17856/1‑149 (MQ=255)
                                                                                                                                                                                                                                                cGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGtt                                                                                                                                                          >  7:237243/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        tgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCa                                                                                                                                                                 <  6:240882/134‑1 (MQ=255)
                                                                                                                                                                                                                                                        tgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCa                                                                                                                                                                 >  5:240882/1‑134 (MQ=255)
                                                                                                                                                                                                                                                        tgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAg                                                                                                                                                  <  7:173055/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        tgAAACAAGGGTGGCTTATGCCGCCTCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGCTATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATGTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGATTCGCATCGGTAACTGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGtt                    <  2:391742/275‑1 (MQ=255)
                                                                                                                                                                                                                                                         gAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAgg                                                                                                                                                 >  7:101651/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         gAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAgg                                                                                                                                                 >  6:373896/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         gAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAgg                                                                                                                                                 >  5:220142/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           aaCAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAg                 <  2:207194/275‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAg                                                                                                                                                       <  7:43471/140‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCAGGGGTGGCATATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCAGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAg                                                                                                                                                       >  8:43471/1‑140 (MQ=255)
                                                                                                                                                                                                                                                             cAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCAt                                                                                                                                             <  8:237243/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             cAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCAt                                                                                                                                             <  8:266921/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             aaGGGGCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCAt                                                                                                                                             >  6:249304/2‑149 (MQ=255)
                                                                                                                                                                                                                                                              aGGGGCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCACTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCCGTTGTATCCCAATTTCATACAGTTAAGGACAGACCATg                                                                                                                                            >  5:162653/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               ggggTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGa                                                                                                                                           >  6:43166/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               ggggCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACa                                                                                                                                                             >  5:304845/1‑131 (MQ=255)
                                                                                                                                                                                                                                                               ggggCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACa                                                                                                                                                             <  6:304845/131‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aTGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGAtttt                                                                                                                                 <  7:375705/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGAttttt                                                                                                                                <  8:101651/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tGCCGCCCCTTATTCCATATTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACCTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACCGGCCATGAGCAAGATTTTTGAAGATCACTCGCTGACTATCGGTCACACGTCGCTGGTTCGCCTGAGTCGCATCGGTAACGAACGAACTCTGCCGAAGGTGGCATCTCGTACCCCCAGTTCCAGCGTTATGGCCCGTAACGGTGca  >  4:140022/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                              
TATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCAAAATCACCAATATGCAGTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCC  >  NC_000913/2532009‑2532549

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: