Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,532,275 T→A 9.5% intergenic (+51/‑134) 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,2750TA 9.5% 302.7 / 19.5 95intergenic (+51/‑134)cysZ/cysKsulfate transporter, sulfite inhibited/cysteine synthase A, O‑acetylserine sulfhydrolase A subunit
Reads supporting (aligned to +/- strand):  ref base T (41/45);  new base A (4/5);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-02

TCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCAAAATCACCAATATGCAGTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCC  >  NC_000913/2532004‑2532549
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                  
tcaacttttGCGGTTACCCCTTCGGTAACCACAAAGTGCCGTTTAAAGATGTTCGCACCGCCCTGCGCACACGCAAAATTACCAATATGCAGTTTGGGGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGGGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgcc                                                                                                                                                                                                                                                                                 <  4:467788/268‑1 (MQ=255)
tCCAGGATTGAGATTACCCCTTCGATAACCCCAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCAAAATCACCAATATGCAGTTTGGTGCATTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgcc                                                                                                                                                                                                                                                                                 <  4:164566/275‑1 (MQ=255)
     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)
                      cGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCAAAATCACCAATATGCAGTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgcc                                                                                                                                                                                                                                                                                 >  1:139537/1‑253 (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)
                                                                                                                               gctgctTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgcc                                                                                                                                                                                                                                                                                 <  1:115321/148‑1 (MQ=255)
                                                                                                                                ctgctTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATACACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTTAGCGGGTATCCTCCCTGTTGTATACCAATTTCCTACAGTTAAGTCCAgg                                                                                                                                                 >  4:207191/1‑275 (MQ=255)
                                                                                                                                    ttAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCt                                                                                                                                                                                                                                                                           <  6:159273/149‑1 (MQ=255)
                                                                                                                                      aaTCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCtgt                                                                                                                                                                                                                                                                         <  6:414085/149‑3 (MQ=255)
                                                                                                                                      aaTCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCtgt                                                                                                                                                                                                                                                                         <  6:223000/149‑3 (MQ=255)
                                                                                                                                               atcatGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCtt                                                                                                                                                                                                                                                                >  6:47143/1‑149 (MQ=255)
                                                                                                                                               atcatGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCtt                                                                                                                                                                                                                                                                >  6:161237/1‑149 (MQ=255)
                                                                                                                                                     cccGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgcc                                                                                                                                                                                                                                                                                 <  4:424397/126‑1 (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:57681/1‑149 (MQ=255)
                                                                                                                                                                                                     gATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCat                                                                                                                                                                                                          >  7:228304/1‑149 (MQ=255)
                                                                                                                                                                                                            acGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGata                                                                                                                                                                                                                                <  8:325507/120‑1 (MQ=255)
                                                                                                                                                                                                            acGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGata                                                                                                                                                                                                                                >  7:325507/1‑120 (MQ=255)
                                                                                                                                                                                                            acGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTctct                                                                                                                                                                                                   >  5:332996/1‑149 (MQ=255)
                                                                                                                                                                                                              gcgATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGa                                                                                                                                                                                                 <  7:238292/149‑1 (MQ=255)
                                                                                                                                                                                                                  tgtgGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCgg                                                                                                                                                                                             <  5:161237/149‑1 (MQ=255)
                                                                                                                                                                                                                  tgtgGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCgg                                                                                                                                                                                             <  5:47143/149‑1 (MQ=255)
                                                                                                                                                                                                                                  tACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCtgttgt                                                                                                                                                                             >  5:143592/1‑149 (MQ=255)
                                                                                                                                                                                                                                     tggTTATTTGGTAACACCTTTGTGTGAAACAGGGGTGGTTGTTGCCGCCCCTTATTCCATTTTGCATGTCATTATTTCCCTGCTATATCTAGATATGCTAAATCCTGACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGa                                             <  4:299537/273‑1 (MQ=255)
                                                                                                                                                                                                                                           ttttGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGccc                                                                                                                                                                                                                                                                             <  6:331172/44‑1 (MQ=255)
                                                                                                                                                                                                                                           ttttGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGccc                                                                                                                                                                                                                                                                             >  5:331172/1‑44 (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                                                                                                                                                                 >  5:240882/1‑134 (MQ=255)
                                                                                                                                                                                                                                                             tgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCa                                                                                                                                                                 <  6:240882/134‑1 (MQ=255)
                                                                                                                                                                                                                                                             tgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAg                                                                                                                                                  <  7:173055/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             tgAAACAAGGGTGGCTTATGCCGCCTCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGCTATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATGTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGATTCGCATCGGTAACTGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGtt                    <  2:391742/275‑1 (MQ=255)
                                                                                                                                                                                                                                                              gAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAgg                                                                                                                                                 >  7:101651/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              gAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAgg                                                                                                                                                 >  5:220142/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              gAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAgg                                                                                                                                                 >  6:373896/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                                                                                                                                                             <  6:304845/131‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ggggCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACa                                                                                                                                                             >  5:304845/1‑131 (MQ=255)
                                                                                                                                                                                                                                                                              aTGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGAtttt                                                                                                                                 <  7:375705/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                               tGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGAttttt                                                                                                                                <  8:101651/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                               tGCCGCCCCTTATTCCATATTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACCTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACCGGCCATGAGCAAGATTTTTGAAGATCACTCGCTGACTATCGGTCACACGTCGCTGGTTCGCCTGAGTCGCATCGGTAACGAACGAACTCTGCCGAAGGTGGCATCTCGTACCCCCAGTTCCAGCGTTATGGCCCGTAACGGTGca  >  4:140022/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                  
TCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCAAAATCACCAATATGCAGTTTGGTGCTTTAACCAGCCTGTTTACGATGATCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCC  >  NC_000913/2532004‑2532549

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