Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,827,811 T→G 5.5% K150T (AAA→ACA)  astB ← succinylarginine dihydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,827,8110TG 5.5% 374.7 / 13.1 110K150T (AAA→ACA) astBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base T (51/53);  new base G (4/2);  total (55/55)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.42e-02

Only 100 of 123 total aligned reads displayed.

CGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGATTCAGCCTTGCCACCGCCTCGCTGGCTTCGCGAGTCTGTCGCGCCGGATAACGGGAAGGCCGGGTATCATTGCCTTCTTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGGAATAAACGGACGCTCGTGCGGCGGGATCA  >  NC_000913/1827539‑1828078
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                            
  tCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGATTCAGCCTTGCCACCGCCTCGCTGGCTTCGCGAGTCTGTCGCGCCGGATAACGGGACGGCCGGGTACCCTTGCCTTCTTCTCGCCCGTAGCCAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGCCGATCGTGGTTTATCGCCCCTTCCTCCCCGCGCCACGCTACCCTGGGCAACGCCCACTGCACGCCCAATTTCTCTTCGCTCTTACAAATTCCTTTTa                                                                                                                                                                                                                                                                           >  4:22172/1‑274 (MQ=255)
  tCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGATTCAGCCTTGCCACCGCCTCGCTGGCTTCGCGAGTCTGTCGCGCCGGATAACGGGAAGGTCGGGTATCATTGCCTTCTTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCCCTTCGTCGTTAACAATCGCTTTTa                                                                                                                                                                                                                                                                           >  2:10947/1‑274 (MQ=255)
                               ttGGGGATTCACTGGTCTCGGCTGTGACACCGCCCCGCTGGCTTCGCTGGTATTTCGCGCCGGATAACGGGAATGCCGGGTATCATTGCCTTCTTCGCGCCCATAGACAAAAAGTTGCATCCCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGC‑TTTAACAGCGATTCAGTGACGGGCGCTTCCAGCg                                                                                                                                                                                                                                             <  4:247507/274‑1 (MQ=255)
                               ttGGGGATTCACCTGATTCAGCCTTGCCACCGCCTCGCTGGCTTCGCGAGTCTGTCGCGCCGGATAACGGGAAGGCCGGGTATCATTGCCTTCTTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGCCGTCAAAAATCGCTTTTAACACCGATTCAGTGACGGGCGCTTACGGCg                                                                                                                                                                                                                                             >  1:118887/1‑275 (MQ=255)
                                 ggggATTCACCTGATTCAGCCTTGCCACCGCCTCGCTGGCTTCGCGAGTCTGTCGCGCCGGATAACGGGAAGGCCGGGTATCATTGCCTTCTTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGaa                                                                                                                                                                                                                                           <  2:448188/275‑1 (MQ=255)
                                     aTTCACCTGATTCAGCCTTGCCACCGCCTCGCTGGCTTCGCGAGTCTGTCGCGCCGGATAACGGGAAGGCCGGGTATCATTGCCTTCTTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGt                                                                                                                                                                                                                                       <  2:122085/275‑1 (MQ=255)
                                         aCCTGATTCAGCCTTGCCACCGCCTCGCTGGCTTCGCGAGTCTGTCGCGCCGGATAACGGGAAGGCCGGGTATCATTGCCTTCTTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTGTCAGCGAACGGTTaat                                                                                                                                                                                                                                   >  2:34653/1‑274 (MQ=255)
                                              aTTCAGCCTTGCCACCGCCTCGCTGGCTTCGCGAGTCTGTCGCGCCGGATAACGGGAAGGCCGGGTATCATTGCCTTCTTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCACTGACGGGCGCTTCCAGCGAACGGGGAAATTTAt                                                                                                                                                                                                                              >  1:54842/1‑275 (MQ=255)
                                              aTTAAGCCTTGCCACCGCCTCGCTGGCGTCGCGAGTCTGTCGCGCCGGATAACGGGAAGGCCGGGTATCATTGCCTTCTTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTa                                                                                                                                                                                                                               <  2:22172/274‑1 (MQ=255)
                                                   acattGCCACCGGCTCGCTGGCTTCGCGGGTTTGTCGCGCCGGATAACGGGAAGGCCGGGTATCATTGCCTTCTTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGtt                                                                                                                                                                                                                          <  2:249758/271‑1 (MQ=255)
                                                      ttgaccccgccgCGGTGGCTTCGGGAGTCTGTCGCGACGGATAACCGGAAGGCCTGGGATCATTGCCCTCTTCTCGCCCGTAGACCAAAAGTTGCATACCCGGTTTGCCGTAAGGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAgg                                                                                                                                                                                                                      <  4:186852/269‑1 (MQ=255)
                                                                         ttCGCGAGTCTGTCGCGCCGGATAACGGGAAGGCCGGGTATCATTGCCTTCTTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCa                                                                                                                                                                                                   <  4:283132/275‑1 (MQ=255)
                                                                          tgggggggcTGTCTCGCCGGATAACGGGGGGGCCGGGTATCATTGCCTGCTTCGCGCCCGTCGACAAAAAGTTTCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCa                                                                                                                                                                                                   <  4:34652/266‑1 (MQ=255)
                                                                                                                   aTTGCCTTCTTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTTACGGGCTCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAAGGGTGGGCTGCACTTTGCCATCCAGCGTATCGGACAG‑TTGCCCGATCGTTTcc                                                                                                                                                         >  4:73673/1‑274 (MQ=255)
                                                                                                                              ctggCCCGTAGACAAAAAGTTGCATAACCTGTTCACCGTAATGACCGCCGAGACGATTGTGGTTGGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGTCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTAc                                                                                                                                              <  2:222816/271‑1 (MQ=255)
                                                                                                                                cGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTa                                                                                                                                                                                                                                                                           >  6:500521/1‑148 (MQ=255)
                                                                                                                                 gCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCAATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCcac                                                                                                                                          >  2:536165/1‑275 (MQ=255)
                                                                                                                                            aaaaGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCtt                                                                                                                                                                                                                                                                              >  5:927344/1‑133 (MQ=255)
                                                                                                                                            aaaaGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCtt                                                                                                                                                                                                                                                                              <  6:927344/133‑1 (MQ=255)
                                                                                                                                            aaaaGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAg                                                                                                                                                                                                                                                               >  7:474021/1‑148 (MQ=255)
                                                                                                                                               aGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGa                                                                                                                                                                                                                                                            >  8:725085/1‑148 (MQ=255)
                                                                                                                                                 ttGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAANAATCGCTNTTANCAGCGANTCAGTGACGGGCGCTTCCAGCGAACGTGGAAATTTATTGTTCAgg                                                                                                                                                                                                                      >  3:24479/1‑184 (MQ=255)
                                                                                                                                                    cATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCtt                                                                                                                                                                                                                                                                              <  8:489487/125‑1 (MQ=255)
                                                                                                                                                     aTACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTgg                                                                                                                                                                                                                  <  1:24479/184‑1 (MQ=255)
                                                                                                                                                      tACCCGGTTCACCTGAATGACCGCCGAGAAGGTTGTGGTTTGCCGCCCCCGCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGAc                                                                                                                      <  4:445635/274‑1 (MQ=255)
                                                                                                                                                        cccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACg                                                                                                                                                                                                                                                          >  5:926523/1‑141 (MQ=255)
                                                                                                                                                         ccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGa                                                                                                                                                                                                                                                                    >  6:1000690/1‑130 (MQ=255)
                                                                                                                                                          cGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGgcg                                                                                                                                                                                                                                                      <  5:1084416/143‑1 (MQ=255)
                                                                                                                                                          cGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGgcg                                                                                                                                                                                                                                                      >  6:1084416/1‑143 (MQ=255)
                                                                                                                                                            gTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAg                                                                                                                                                                                                                                               <  5:500521/148‑1 (MQ=255)
                                                                                                                                                             ttCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGc                                                                                                                                                                                                                                              >  6:260574/1‑148 (MQ=255)
                                                                                                                                                              tCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTCAATTTCTCTTCGTCGTTAAAAATCGCTTTCAACAGCGATTCAGTGACGGGCGCTTCCAGCGCACGGGTAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CCAATGGGGCGATGGTTGCAGCATTGGCAACCCACATTGCGGAGACGGAACGGAAGCGGGAAAg                                                                                                             >  4:277882/1‑275 (MQ=255)
                                                                                                                                                               taCCGTAAAGACCGGCTAGACGATTGGGGTTTGCCTCCCCCTCATCACCGAGCAACGGTACCTGTGGCAACGCCGAATGGAAGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGc                                                                                                            <  4:394780/274‑1 (MQ=255)
                                                                                                                                                                     aaTGACCGCCGAGACGAATATGGTTTGCCGCCCCCTCTTCACCGAGCCACCCTACCTGTGGCAACGCCGAATGGACGCTCAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCTATTCAGTGACGTGCGCTGCCAGCGACCCGTGTAATTCATTCTTCGGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTCACGG‑CGGATGGCGAGATCGTTGCAGCATTGGCCACCGACATTCGCGAACCGGACCTTACGATTGACAGCCAGTGc                                                                                                      >  4:164571/1‑275 (MQ=255)
                                                                                                                                                                         aCCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc                                                                                                                          >  2:304790/1‑251 (MQ=255)
                                                                                                                                                                           cgccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGcc                                                                                                >  2:271044/1‑275 (MQ=255)
                                                                                                                                                                                agaCGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCa                                                                                                                                                                                                                                                                <  5:29033/111‑1 (MQ=255)
                                                                                                                                                                                   cGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTg                                                                                         <  2:107923/274‑1 (MQ=255)
                                                                                                                                                                                     aTTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAgg                                                                                                                                                                                                                      >  7:801477/1‑148 (MQ=255)
                                                                                                                                                                                        ttgGTTTGCCGCCCCCCGATCACCGGGCAACGATACCTGGGGCAACGCCGAATGGACGCTAAATGTCGCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACt                                                                                    <  4:252775/273‑1 (MQ=255)
                                                                                                                                                                                            tttGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGa                                                                                                                                                                                                                                            <  5:1083336/119‑1 (MQ=255)
                                                                                                                                                                                            tttGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGa                                                                                                                                                                                                                                            >  6:1083336/1‑119 (MQ=255)
                                                                                                                                                                                                ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAAt                                                                                                                                                                                                                                  <  8:724896/125‑1 (MQ=255)
                                                                                                                                                                                                ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAAt                                                                                                                                                                                                                                  >  7:724896/1‑125 (MQ=255)
                                                                                                                                                                                                ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAg                                                                            >  2:355253/1‑274 (MQ=255)
                                                                                                                                                                                                  gccCCCTCATCCCCGAGCAACGCTACCTGGGGCAAATCCGAATGGACGCCAAATGTCGTTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGTGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTAc                                                                         <  4:10947/275‑1 (MQ=255)
                                                                                                                                                                                                  gccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTga                                                                                                                                                                                                         <  8:801477/148‑1 (MQ=255)
                                                                                                                                                                                                    ccccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTgaga                                                                                                                                                                                                       <  5:260574/148‑1 (MQ=255)
                                                                                                                                                                                                     cccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCt                                                                       >  4:333154/1‑274 (MQ=255)
                                                                                                                                                                                                        tcatcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCtt                                                                                                                                                                                                                                                   <  5:1041364/100‑1 (MQ=255)
                                                                                                                                                                                                          atcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTcc                                                                                                                                                                                                                                                 >  8:1036418/1‑100 (MQ=255)
                                                                                                                                                                                                          atcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTcc                                                                                                                                                                                                                                                 <  7:1036418/100‑1 (MQ=255)
                                                                                                                                                                                                          atcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACt                                                                                                                                                                                                 >  6:613988/1‑148 (MQ=255)
                                                                                                                                                                                                           tcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTc                                                                                                                                                                                                                         <  5:383131/123‑1 (MQ=255)
                                                                                                                                                                                                           tcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTc                                                                                                                                                                                                                         >  6:383131/1‑123 (MQ=255)
                                                                                                                                                                                                           taaCCGAGCAACGTTAACTCTGGCAACGCCGAATGGACTCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCAt                                                                 <  4:20082/272‑1 (MQ=255)
                                                                                                                                                                                                                 aGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAACTTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGCACGGTGAACTTTATTGTTCAGGTTGGCAACGGTGCGATCCACTTTGCAATCCAGCGTATCGG‑CCGATGGCGCGATCGTCGCCGCCTTGGCTACCCCCATTGGCGCAGAGGAACTGCCGCTGGAAAGCCAGGTAGGTGTCTGGGGTCAACCTTTTTCCAGTAACTGCCCATCGCCAc                                                           >  4:21696/1‑274 (MQ=255)
                                                                                                                                                                                                                   cAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAt                                                                                                                                                                                                                              >  8:702031/1‑110 (MQ=255)
                                                                                                                                                                                                                   cAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAt                                                                                                                                                                                                                              <  7:702031/110‑1 (MQ=255)
                                                                                                                                                                                                                   cAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATcc                                                                                                                                                                                        <  5:613988/148‑1 (MQ=255)
                                                                                                                                                                                                                      cGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGa                                                     >  1:3005/1‑275 (MQ=255)
                                                                                                                                                                                                                      cGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAACGCCAGTGCGGTGCCTGGCGTTCAACTTTTTCCAGTACCTGCTCAACCCTCCCACTg                                                      >  4:283259/1‑274 (MQ=255)
                                                                                                                                                                                                                                 gCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGa                                                                                                                                                                         <  7:290607/148‑1 (MQ=255)
                                                                                                                                                                                                                                    aCGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACa                                                                                                                                                                                                                                                                        >  7:715513/1‑51 (MQ=255)
                                                                                                                                                                                                                                        cGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACgcagca                                   <  2:248219/275‑1 (MQ=255)
                                                                                                                                                                                                                                         gAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCAc                                  <  2:73674/275‑1 (MQ=255)
                                                                                                                                                                                                                                           aTGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCg                                <  2:470909/275‑1 (MQ=255)
                                                                                                                                                                                                                                                cGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGGAAt                            >  2:317387/1‑274 (MQ=255)
                                                                                                                                                                                                                                                  cTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTga                                                                                                                                                                                                         <  8:710612/100‑1 (MQ=255)
                                                                                                                                                                                                                                                  cTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTga                                                                                                                                                                                                         >  7:710612/1‑100 (MQ=255)
                                                                                                                                                                                                                                                  cTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGAt                                                                                                                                                                <  7:1112427/140‑1 (MQ=255)
                                                                                                                                                                                                                                                          ctctTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGt                                                                                                                                                             <  7:436401/135‑1 (MQ=255)
                                                                                                                                                                                                                                                          ctctTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGt                                                                                                                                                             >  8:436401/1‑135 (MQ=255)
                                                                                                                                                                                                                                                            ctTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGa                                                                                                                                                                                                                                            <  8:330383/55‑1 (MQ=255)
                                                                                                                                                                                                                                                              tcgtcgTTAAAAATCGCTTTTAACAGCGATTCAGTGTCGGGCGCTTCCAGCGAACGGTGAAATTTAGTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGGAATAAACGGACGCTCGt              <  1:67115/274‑1 (MQ=255)
                                                                                                                                                                                                                                                              tcgtcgTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACg                                                                                                                                                                                                                                         <  8:235077/56‑1 (MQ=255)
                                                                                                                                                                                                                                                              tcgtcgTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACg                                                                                                                                                                                                                                         >  7:235077/1‑56 (MQ=255)
                                                                                                                                                                                                                                                              tcgtcgTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTAccc                                                                                                                                            >  8:149344/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                 tcgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCtt                                                                                                                                                                                                                                                   <  7:538550/43‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tcgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCtt                                                                                                                                                                                                                                                   >  8:538550/1‑43 (MQ=255)
                                                                                                                                                                                                                                                                 tcgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATg                                                                                                                                                                                                     >  6:512682/1‑89 (MQ=255)
                                                                                                                                                                                                                                                                     tAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACtt                                                                                                                                                                                                >  8:238343/1‑90 (MQ=255)
                                                                                                                                                                                                                                                                     tAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACtt                                                                                                                                                                                                <  7:238343/90‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTgg                                                                                                                                     >  7:289284/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                        gaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCg                                                                                                                                                                               >  8:1040591/2‑104 (MQ=255)
                                                                                                                                                                                                                                                                        gaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCg                                                                                                                                                                               <  7:1040591/103‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGa                                                                                                                                  >  7:473340/2‑148 (MQ=255)
                                                                                                                                                                                                                                                                        aaaTCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGa                                                                                                                                  <  8:289284/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aagcGCTTGTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGAAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGGAATAAACGGACGCTCGTGCGGCGGGATCa  <  2:176817/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aaTCGCGTTTAACAGCGGTTCAGGGACGGGCGCTTCCAGGGAACGGGGAAATTTATTGTTCAGGTTGGCAACGGTGGGATGCACTTTGCCAACCAGCGTATCGG‑CAGGTGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGGAATAAACGGACGCGCGTGCGGCGGGATCa  <  4:288686/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                            cGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCg                                                                                                                              >  8:29422/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                            cGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACAtt                                                                                                                                       <  8:405539/139‑1 (MQ=255)
                                                                                                                                                                                                                                                                            cGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACAtt                                                                                                                                       >  7:405539/1‑139 (MQ=255)
                                                                                                                                                                                                                                                                            cGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCg                                                                                                                              >  8:422702/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                               ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGaa                                                                                                                           >  8:735992/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                tttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGc                                                                                                                                                                                                    <  6:938116/75‑1 (MQ=255)
                                                                                                                                                                                                                                                                                tttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGc                                                                                                                                                                                                    >  5:938116/1‑75 (MQ=255)
                                                                                                                                                                                                                                                                                tttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCAt                                                                                                                                                                                          >  5:205769/1‑85 (MQ=255)
                                                                                                                                                                                                                                                                                tttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCAt                                                                                                                                                                                          <  6:205769/85‑1 (MQ=255)
                                                                                                                                                                                                                                                                                tttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCgg                                                                                                                                                                              <  7:915149/97‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                            
CGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGATTCAGCCTTGCCACCGCCTCGCTGGCTTCGCGAGTCTGTCGCGCCGGATAACGGGAAGGCCGGGTATCATTGCCTTCTTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGG‑CAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGGAATAAACGGACGCTCGTGCGGCGGGATCA  >  NC_000913/1827539‑1828078

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