Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,213,156 G→A 5.1% noncoding (117/120 nt) rrfE → 5S ribosomal RNA of rrnE operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,213,1560GA 5.1% 297.5 / 5.3 97noncoding (117/120 nt)rrfE5S ribosomal RNA of rrnE operon
Reads supporting (aligned to +/- strand):  ref base G (46/46);  new base A (3/2);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.24e-01

TGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCAT  >  NC_000913/4212887‑4213416
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                    
tgtgtAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATc                                                                                                                                                                                                                                                                  >  1:67058/1‑274 (MQ=31)
tgtgtAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCCCCCAATGCGAGAGTAGGGAACTGCCATGCATc                                                                                                                                                                                                                                                                  >  4:220945/1‑274 (MQ=14)
         gcAGCGATGCGGTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAGTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGaa                                                                                                                                                                                                                                                         <  2:483331/274‑1 (MQ=17)
                 ggtttGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACGGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCg                                                                                                                                                                                                                                                 <  2:536212/271‑1 (MQ=25)
                  cGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACGGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGAGCGGTGGTACCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCTGAGCCCCAGTGGTAGGTGGGGGTCTCCCATTTCGAGAGTGGGGAACGTCAAGGCCGCAAATTCGAAAAAcccaga                                                                                                                                                                                                                                                >  4:536211/1‑271 (MQ=11)
                   acggAGCTAACATATACTAATGAACGNNGAGGCTTAACCTTACAACGCCGAAGGTGTTNTGACGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATNTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCCTGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCg                                                                                                                                                                                                                                                 <  3:68209/269‑1 (MQ=255)
                              cGGGACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCg                                                                                                                                                                                                                                    <  2:220946/274‑1 (MQ=255)
                                 tACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAAc                                                                                                                                                                                                                                                     >  1:68209/1‑254 (MQ=34)
                                 tACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTAAGAAGTGAAAAGCCGTAGCGCCGAGGGTAGTGTGGGGTGGCCCCATGCGAGAGTAGGGGACTGCCAGGAACGAANNNAGAAAAACCCCGGGCCAAAAGGACgggg                                                                                                                                                                                                                                 >  3:141661/1‑274 (MQ=11)
                                      aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGtttttt                                                                                                                                                                                                                           >  2:506298/1‑275 (MQ=255)
                                           ccGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCata                                                                                                                                                                                                                      >  2:313273/1‑275 (MQ=255)
                                           ccGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCata                                                                                                                                                                                                                      >  2:456708/1‑275 (MQ=255)
                                              tGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGTGTTTTTTGCATATCa                                                                                                                                                                                                                   >  2:124521/1‑275 (MQ=255)
                                                   cgtAACCTTCCAACGTCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTAAGCCTTATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTGCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGGGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAattat                                                                                                                                                                                                              <  4:400416/273‑1 (MQ=255)
                                                        ccTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCa                                                                                                                                                                                                         <  2:175331/275‑1 (MQ=255)
                                                            aCAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCatgat                                                                                                                                                                                                     >  2:467707/1‑275 (MQ=255)
                                                                cGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGACGTGAAACGCCGTAGCGCCGATGGTAGGGTGGGGTCTCCCCAGGCGAGAGGAGGGAACGGCCGGGCATCAAATTAGACAAACCCCGGTCCATAAGGCCGGGGTTTGTTGCCTACCAATTATTTGCATAATTAAg                                                                                                                                                                                                 >  4:81496/1‑275 (MQ=17)
                                                                  ccGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATGCAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTGCCAGCTGACCCCATGCCGAACTCAGAAGTGAAACGCCTTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGCAGGGAACTGGCAGGCATCAAATTGGAAAAACCCGGGTCCATAAGGCCGGGGTTTTTGGCATATCAATTATTTTCatgtttccttg                                                                                                                                                                                               >  4:398297/1‑267 (MQ=17)
                                                                      aGGAGTGTTGGAGGATTGAGAGAAGATTTTCAGCCTGATGCAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATATGCCTGGCGGCGGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGGGGAGGGAACTGCCAGGCATCAAATTAAAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAattat                                                                                                                                                                                                              <  2:398298/256‑1 (MQ=17)
                                                                          gtTTTGGCGGATTGAGAGAAGATTTTCAGCCGGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCAt                                                                                                                                                                                       <  2:434671/275‑1 (MQ=255)
                                                                           ttttGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGCAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGTAATCCCAt                                                                                                                                                                                       >  2:266327/1‑274 (MQ=255)
                                                                            tttGGGGAATTGAGAGAAGATTTTCAGCCTGATAAAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGACTTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGt                                                                                                                                                                                     <  1:151147/275‑1 (MQ=255)
                                                                            tttGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGt                                                                                                                                                                                     <  4:124520/275‑1 (MQ=255)
                                                                            tttGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGAGTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAAGTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCAAGCGAGAGGAGGGAACTGCCAGGAAACAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGt                                                                                                                                                                                     <  2:511688/275‑1 (MQ=255)
                                                                                        agtgaagaTTTTCAGACAGATACAGATNNAAAAAGAGAGCGGAAGCGAGCTGATAAAACNGAATTTGCCTGGGGGAAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGNCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGta                                                                                                                                                                           <  3:67058/270‑1 (MQ=18)
                                                                                                         tGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTaa                                                                                                                                                         <  2:430201/274‑1 (MQ=255)
                                                                                                                 aTTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTAc                                                                                                                                                >  4:388789/1‑275 (MQ=255)
                                                                                                                 aTTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGAACCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                                                                                                                                                                                                                                               <  2:558725‑M2/164‑2 (MQ=255)
                                                                                                                      aTCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCCTCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATACCCACATTCGTCGTAATCAAATTACTGTTc                                                                                                                                           >  1:132844/1‑275 (MQ=255)
                                                                                                                           aaCGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGAt                                                                                                                                      >  2:347507/1‑275 (MQ=255)
                                                                                                                                     cGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaa                                                                                                                                >  2:1713/1‑271 (MQ=255)
                                                                                                                                     cGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCa                                                                                                                            <  4:1713/275‑1 (MQ=255)
                                                                                                                                       gTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATTTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTTTTTGCGATCTTGAAAACaaa                                                                                                                          >  4:425477/1‑275 (MQ=255)
                                                                                                                                       gTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGACGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTTAGTTCTGTATATCCAAATTCGGCGTAATCAAATTATTGTTCGACATTTTGAAATCaaa                                                                                                                          >  4:486945/1‑275 (MQ=255)
                                                                                                                                              aaaaCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa                                                                                                                    <  2:388790/274‑1 (MQ=255)
                                                                                                                                              aaaaCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAg                                                                                                                   <  4:506297/275‑1 (MQ=255)
                                                                                                                                              aaaaCAGAAATTGCCGGGCGTCAGGATCGCGGGGGGCCCACCTGACCCCAGTCCGAACTCAAAAGGGAAACGCCGTAGCGCCGATGGTAGTGTGTGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa                                                                                                                    <  4:313272/274‑1 (MQ=37)
                                                                                                                                                    ccaTTTGCCGGGCGGCGGTAGCTGTGTGGTCCCACGTGACCCCAGTCCGAAATCAGAAGTGAAACGCAGTNGCGCCGATGGTAGTGGGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATGAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCGCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa                                                                                                                    <  3:115048/266‑1 (MQ=37)
                                                                                                                                                           ccggggggcAGGTGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCTGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACg                                                                                                       <  4:266326/268‑1 (MQ=255)
                                                                                                                                                            cTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCt                                                                                                     <  2:425478/275‑1 (MQ=255)
                                                                                                                                                                    aGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACATCAAataaaacaaaaggcccagtcttccgactgagccttt                                                                                                                                                                                                                            >  7:47975‑M2/1‑112 (MQ=255)
                                                                                                                                                                           cgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaa                                                                                                                                                                                                                                                           >  7:1034728‑M2/1‑105 (MQ=255)
                                                                                                                                                                           cgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaa                                                                                                                                                                                                                                                           <  8:1034728‑M2/110‑6 (MQ=255)
                                                                                                                                                                              tggtCCCACCTGACCCCATGCCGAACTCAGAGGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGGGTAGGGAACTGCTAGGCATCAAATTCGAAAAACCCCGGTACATAAGGCCGGGGTTTTTNGCATATCAATTATTTGAATGATGAAGGGAAGCTCATGTCAGTTCGGTATATCCAAATTCGTCGTAATCAAATTACTGGTCNCGATCTTGCAGGCAAACGTGAAATTNAAGGAGACGCCGAATTCCATAACCAgtgg                                                                                   >  3:91945/1‑272 (MQ=25)
                                                                                                                                                                              tggtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                                                                                                                                                                                                                                               <  6:343197‑M2/103‑2 (MQ=255)
                                                                                                                                                                              tggtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                                                                                                                                                                                                                                               >  5:343197‑M2/1‑102 (MQ=255)
                                                                                                                                                                                gtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaa                                                                                                                                                                                                                                                           <  5:265316‑M2/105‑6 (MQ=255)
                                                                                                                                                                                gtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaa                                                                                                                                                                                                                                                           >  6:265316‑M2/1‑100 (MQ=255)
                                                                                                                                                                                  cccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACTTCAAataaaacaaaaggcccagtcttccgactgagcc                                                                                                                                                                                                                               >  8:1038098‑M2/1‑98 (MQ=255)
                                                                                                                                                                                  cccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACATCAActaaaacaaaaggcccagtcttccgactgagcc                                                                                                                                                                                                                               <  7:1038098‑M2/131‑34 (MQ=255)
                                                                                                                                                                                   ccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCgg                                                                                                                                                                                                                                   >  5:128625/1‑126 (MQ=255)
                                                                                                                                                                                   ccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCgg                                                                                                                                                                                                                                   <  6:128625/126‑1 (MQ=255)
                                                                                                                                                                                     aCCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTGGCGCCGATGGTTGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTaa                                                                            <  1:66931/275‑1 (MQ=255)
                                                                                                                                                                                         gACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaa                                                                                                                                                                                                                                                           <  6:119114‑M2/96‑6 (MQ=255)
                                                                                                                                                                                         gACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaa                                                                                                                                                                                                                                                           >  5:119114‑M2/1‑91 (MQ=255)
                                                                                                                                                                                          aCCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTAGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGcc                                                                       >  4:92375/1‑275 (MQ=255)
                                                                                                                                                                                             ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                                                                                                                                                                                                                                               >  6:70762‑M2/1‑87 (MQ=255)
                                                                                                                                                                                             ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                                                                                                                                                                                                                                               <  5:70762‑M2/88‑2 (MQ=255)
                                                                                                                                                                                              cattcccaatccAGAAGTGAAACGCCGTAGCGCCAATGGTAGTTTGGGGTCCCCCCAGGCGGGAGTAGGGAACGGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAAt                                                                    <  4:347506/263‑1 (MQ=255)
                                                                                                                                                                                              cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggcccagtcttccgactga                                                                                                                                                                                                                                  <  6:491071‑M2/116‑31 (MQ=255)
                                                                                                                                                                                              cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggcccagtcttccgactga                                                                                                                                                                                                                                  >  5:491071‑M2/1‑86 (MQ=255)
                                                                                                                                                                                              cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                                                                                                                                                                                                                                               <  5:1103330‑M2/87‑2 (MQ=255)
                                                                                                                                                                                              cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                                                                                                                                                                                                                                               >  6:1103330‑M2/1‑86 (MQ=255)
                                                                                                                                                                                              cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGt                                                                                                                                                                                                                                               <  6:680080/103‑1 (MQ=255)
                                                                                                                                                                                              cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGt                                                                                                                                                                                                                                               >  5:680080/1‑103 (MQ=255)
                                                                                                                                                                                               aTGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTATGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGCTCTTTAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCCAttt                                                                  >  1:111886/1‑275 (MQ=255)
                                                                                                                                                                                               aTGCAGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGCAGTGTGGGGTCTCCCCATGCGAGAGCAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTAATGGNATGAGGCAGGGAATCTCATGTGAGTCCGGTATATGAAAATTCGTCGTAAACAAATTAATGTTCGAGATCTTTAAAGCAAANGAGAAGGGTAAGGGGACGGTGCCGGNAGTAACAAGAGCCGGGTAATCGCCAAatt                                                                  >  3:86102/1‑275 (MQ=25)
                                                                                                                                                                                                     aaCTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATa                                                                                                                                                                                                                                          >  6:918964/1‑101 (MQ=255)
                                                                                                                                                                                                     aaCTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATa                                                                                                                                                                                                                                          <  5:918964/101‑1 (MQ=255)
                                                                                                                                                                                                     aaCTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACTCTGCCTTCCGTACCCAGCGCCTGTTACTCGCCAATTTCttt                                                              >  4:49288/1‑273 (MQ=255)
                                                                                                                                                                                                           gAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACATCAAataaaacaaaaggcccagtcttccgactg                                                                                                                                                                                                                                   <  8:881736‑M2/102‑30 (MQ=255)
                                                                                                                                                                                                           gAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACATCAAataaaacaaaaggcccagtcttccgactg                                                                                                                                                                                                                                   >  7:881736‑M2/1‑73 (MQ=255)
                                                                                                                                                                                                            aaGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTaa                                                                                          >  1:119018/1‑238 (MQ=255)
                                                                                                                                                                                                                 aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTTTTTCAGAAAAAGt                                                >  2:27898/1‑275 (MQ=255)
                                                                                                                                                                                                                  aaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                                                                                                                                                                                                                                               >  7:965349‑M2/1‑66 (MQ=255)
                                                                                                                                                                                                                  aaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                                                                                                                                                                                                                                               <  8:965349‑M2/67‑2 (MQ=255)
                                                                                                                                                                                                                    cGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggcccagtctttcgactgagcct                                                                                                                                                                                                                              >  5:358880‑M2/1‑64 (MQ=255)
                                                                                                                                                                                                                    cGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggcccagtctttcgactgagcct                                                                                                                                                                                                                              <  6:358880‑M2/98‑35 (MQ=255)
                                                                                                                                                                                                                    cGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaa                                                                                                                                                                                                                                                           >  5:905020‑M2/1‑64 (MQ=255)
                                                                                                                                                                                                                    cGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaa                                                                                                                                                                                                                                                           <  6:905020‑M2/69‑6 (MQ=255)
                                                                                                                                                                                                                      ccGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAATCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTg                                            >  1:42447/1‑274 (MQ=255)
                                                                                                                                                                                                                       cGTAGGGCCAATGGTAGTGTGGGGTTTCCCCAAGCTGGATTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAATGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGGGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGGTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCa                                          <  4:113986/275‑1 (MQ=255)
                                                                                                                                                                                                                                     gaggtggtggTTTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTg                             <  4:386917/266‑1 (MQ=255)
                                                                                                                                                                                                                                              ttctcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTccacc                   <  1:141661/274‑1 (MQ=255)
                                                                                                                                                                                                                                               tctcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACg                                                                                                       >  2:564875/1‑190 (MQ=255)
                                                                                                                                                                                                                                                ctcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGc                                                                                                                                                                                                                         <  8:80187/75‑1 (MQ=255)
                                                                                                                                                                                                                                                ctcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGc                                                                                                                                                                                                                         >  7:80187/1‑75 (MQ=255)
                                                                                                                                                                                                                                                  cctcATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggcccagtcttccgactgagcc                                                                                                                                                                                                                               <  8:926911‑M2/67‑34 (MQ=255)
                                                                                                                                                                                                                                                  cctcATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggcccagtcttccgactgagcc                                                                                                                                                                                                                               >  7:926911‑M2/1‑34 (MQ=255)
                                                                                                                                                                                                                                                   cccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGc                                                                                                   <  4:564874/190‑1 (MQ=255)
                                                                                                                                                                                                                                                   cccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTc                                              >  4:362648/1‑243 (MQ=255)
                                                                                                                                                                                                                                                     cacgCGAGTGTAGGGAACTGCCCAGGATCAAATTAGACAAACCCCGGTCCATAAGGCCCGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGtt             <  4:27897/271‑1 (MQ=255)
                                                                                                                                                                                                                                                       tGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCa                                          <  2:362649/243‑1 (MQ=255)
                                                                                                                                                                                                                                                       tGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAAAGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAACGTTCTGCAAGATCTTGTTCTGCAACTCCACCCAAGTTCaa          >  1:112877/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                ggggAAGTGGCAGGGTGCAAATGAGGAAAACCCCGGTCGATGAGGCGGGGGTGTTTGGTGGATCAATTATGTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCAt  <  4:508642/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ggAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAAtctc                                                                                                                                                                                         >  3:25822/1‑89 (MQ=255)
                                                                                                                                                                                                                                                                      cTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGt                                                                                                                                                                                     <  1:25822/89‑1 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                    
TGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTTGTTCTGCAACTCCACCCACGTTCAACCTGGCAT  >  NC_000913/4212887‑4213416

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