Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,827,816 C→A 5.3% L148L (CTG→CTT astB ← succinylarginine dihydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,827,8160CA 5.3% 400.4 / 10.4 114L148L (CTG→CTTastBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base C (52/56);  new base A (4/2);  total (56/58)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-02

Only 100 of 116 total aligned reads displayed.

TTGGGGATTCACCTGATTCAGCCTTGCCACCGCCTCGCTGGCTTCGCGAGTCTGTCGCGCCGGATAACGGGAAGGCCGGGTATCATTGCCTTCTTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGGAATAAACGGACGCTCGTGCGGCGGGATCA  >  NC_000913/1827570‑1828078
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                      
ttGGGGATTCACTGGTCTCGGCTGTGACACCGCCCCGCTGGCTTCGCTGGTATTTCGCGCCGGATAACGGGAATGCCGGGTATCATTGCCTTCTTCGCGCCCATAGACAAAAAGTTGCATCCCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCg                                                                                                                                                                                                                                            <  4:247507/274‑1 (MQ=255)
ttGGGGATTCACCTGATTCAGCCTTGCCACCGCCTCGCTGGCTTCGCGAGTCTGTCGCGCCGGATAACGGGAAGGCCGGGTATCATTGCCTTCTTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGCCGTCAAAAATCGCTTTTAACACCGATTCAGTGACGGGCGCTTACGGCg                                                                                                                                                                                                                                            >  1:118887/1‑275 (MQ=255)
  ggggATTCACCTGATTCAGCCTTGCCACCGCCTCGCTGGCTTCGCGAGTCTGTCGCGCCGGATAACGGGAAGGCCGGGTATCATTGCCTTCTTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGaa                                                                                                                                                                                                                                          <  2:448188/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)
                                                                                        gcttattCGCACCCGTGGACAAAAAGTGGCCTACCCGTTGCACCGGAAGGACCACCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTATTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCAt                                                                                                                                                     <  4:214007/273‑1 (MQ=255)
                                                                                               ctggCCCGTAGACAAAAAGTTGCATAACCTGTTCACCGTAATGACCGCCGAGACGATTGTGGTTGGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGTCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTAc                                                                                                                                              <  2:222816/271‑1 (MQ=255)
                                                                                                             aaaaGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAg                                                                                                                                                                                                                                                              >  7:474021/1‑148 (MQ=255)
                                                                                                                aGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGa                                                                                                                                                                                                                                                           >  8:725085/1‑148 (MQ=255)
                                                                                                                  ttGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAANAATCGCTNTTANCAGCGANTCAGTGACGGGCGCTTCCAGCGAACGTGGAAATTTATTGTTCAgg                                                                                                                                                                                                                     >  3:24479/1‑184 (MQ=255)
                                                                                                                    gCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACggg                                                                                                                                                                                                                                                       <  7:725085/148‑1 (MQ=255)
                                                                                                                      aTACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTgg                                                                                                                                                                                                                 <  1:24479/184‑1 (MQ=255)
                                                                                                                       tACCCGGTTCACCTGAATGACCGCCGAGAAGGTTGTGGTTTGCCGCCCCCGCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGAc                                                                                                                      <  4:445635/274‑1 (MQ=255)
                                                                                                                         cccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACg                                                                                                                                                                                                                                                         >  5:926523/1‑141 (MQ=255)
                                                                                                                         cccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACg                                                                                                                                                                                                                                                         <  6:926523/141‑1 (MQ=255)
                                                                                                                          ccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGa                                                                                                                                                                                                                                                                   <  5:1000690/130‑1 (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)
                                                                                                                               tCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTCAATTTCTCTTCGTCGTTAAAAATCGCTTTCAACAGCGATTCAGTGACGGGCGCTTCCAGCGCACGGGTAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCCAATGGGGCGATGGTTGCAGCATTGGCAACCCACATTGCGGAGACGGAACGGAAGCGGGAAAg                                                                                                             >  4:277882/1‑275 (MQ=255)
                                                                                                                                taCCGTAAAGACCGGCTAGACGATTGGGGTTTGCCTCCCCCTCATCACCGAGCAACGGTACCTGTGGCAACGCCGAATGGAAGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGc                                                                                                            <  4:394780/274‑1 (MQ=255)
                                                                                                                                      aaTGACCGCCGAGACGAATATGGTTTGCCGCCCCCTCTTCACCGAGCCACCCTACCTGTGGCAACGCCGAATGGACGCTCAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCTATTCAGTGACGTGCGCTGCCAGCGACCCGTGTAATTCATTCTTCGGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTCACGGCGGATGGCGAGATCGTTGCAGCATTGGCCACCGACATTCGCGAACCGGACCTTACGATTGACAGCCAGTGc                                                                                                      >  4:164571/1‑275 (MQ=255)
                                                                                                                                          aCCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc                                                                                                                          >  2:304790/1‑251 (MQ=255)
                                                                                                                                            cgccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAAttt                                                                                                                                                                                                                               <  8:474021/148‑1 (MQ=255)
                                                                                                                                            cgccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGcc                                                                                                >  2:271044/1‑275 (MQ=255)
                                                                                                                                              ccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGAc                                                                                                                      <  4:304789/251‑1 (MQ=255)
                                                                                                                                                 agaCGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCa                                                                                                                                                                                                                                                               >  6:29033/1‑111 (MQ=255)
                                                                                                                                                 agaCGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCa                                                                                                                                                                                                                                                               <  5:29033/111‑1 (MQ=255)
                                                                                                                                                    cGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTg                                                                                         <  2:107923/274‑1 (MQ=255)
                                                                                                                                                      aTTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAgg                                                                                                                                                                                                                     >  7:801477/1‑148 (MQ=255)
                                                                                                                                                         ttgGTTTGCCGCCCCCCGATCACCGGGCAACGATACCTGGGGCAACGCCGAATGGACGCTAAATGTCGCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACt                                                                                    <  4:252775/273‑1 (MQ=255)
                                                                                                                                                             tttGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGa                                                                                                                                                                                                                                           <  5:1083336/119‑1 (MQ=255)
                                                                                                                                                             tttGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGa                                                                                                                                                                                                                                           >  6:1083336/1‑119 (MQ=255)
                                                                                                                                                                 ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAAt                                                                                                                                                                                                                                 >  7:724896/1‑125 (MQ=255)
                                                                                                                                                                 ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAAt                                                                                                                                                                                                                                 <  8:724896/125‑1 (MQ=255)
                                                                                                                                                                 ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAg                                                                            >  2:355253/1‑274 (MQ=255)
                                                                                                                                                                   gccCCCTCATCCCCGAGCAACGCTACCTGGGGCAAATCCGAATGGACGCCAAATGTCGTTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGTGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTAc                                                                         <  4:10947/275‑1 (MQ=255)
                                                                                                                                                                   gccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTga                                                                                                                                                                                                        <  8:801477/148‑1 (MQ=255)
                                                                                                                                                                      cccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCt                                                                       >  4:333154/1‑274 (MQ=255)
                                                                                                                                                                         tcatcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCtt                                                                                                                                                                                                                                                  >  6:1041364/1‑100 (MQ=255)
                                                                                                                                                                         tcatcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCtt                                                                                                                                                                                                                                                  <  5:1041364/100‑1 (MQ=255)
                                                                                                                                                                           atcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTcc                                                                                                                                                                                                                                                <  7:1036418/100‑1 (MQ=255)
                                                                                                                                                                           atcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTcc                                                                                                                                                                                                                                                >  8:1036418/1‑100 (MQ=255)
                                                                                                                                                                           atcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACt                                                                                                                                                                                                >  6:613988/1‑148 (MQ=255)
                                                                                                                                                                            tcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTc                                                                                                                                                                                                                        <  5:383131/123‑1 (MQ=255)
                                                                                                                                                                                  aGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAACTTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGCACGGTGAACTTTATTGTTCAGGTTGGCAACGGTGCGATCCACTTTGCAATCCAGCGTATCGGCCGATGGCGCGATCGTCGCCGCCTTGGCTACCCCCATTGGCGCAGAGGAACTGCCGCTGGAAAGCCAGGTAGGTGTCTGGGGTCAACCTTTTTCCAGTAACTGCCCATCGCCAc                                                           >  4:21696/1‑274 (MQ=255)
                                                                                                                                                                                    cAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAt                                                                                                                                                                                                                             <  7:702031/110‑1 (MQ=255)
                                                                                                                                                                                    cAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATcc                                                                                                                                                                                       <  5:613988/148‑1 (MQ=255)
                                                                                                                                                                                       cGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCCTCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACtta                                                     >  1:96053/1‑273 (MQ=255)
                                                                                                                                                                                       cGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGa                                                     >  1:3005/1‑275 (MQ=255)
                                                                                                                                                                                           aCCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTAt                                                                                                                                                                                >  8:290607/1‑148 (MQ=255)
                                                                                                                                                                                                  gCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGa                                                                                                                                                                         <  7:290607/148‑1 (MQ=255)
                                                                                                                                                                                                     aCGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACa                                                                                                                                                                                                                                                                       <  8:715513/51‑1 (MQ=255)
                                                                                                                                                                                                     aCGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACa                                                                                                                                                                                                                                                                       >  7:715513/1‑51 (MQ=255)
                                                                                                                                                                                                         cGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACgcagca                                   <  2:248219/275‑1 (MQ=255)
                                                                                                                                                                                                          gAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCAc                                  <  2:73674/275‑1 (MQ=255)
                                                                                                                                                                                                          gAATGGACGCTAAATTTCTCTTCGTCGTAAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCGGCGCACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCACTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCCCATAGCTACCACTGACTCCCAACTGACGcaccac                                  >  2:90084/1‑275 (MQ=255)
                                                                                                                                                                                                            aTGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCg                                <  2:470909/275‑1 (MQ=255)
                                                                                                                                                                                                                 cGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGGAAt                            >  2:317387/1‑274 (MQ=255)
                                                                                                                                                                                                                   cTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTga                                                                                                                                                                                                        <  8:710612/100‑1 (MQ=255)
                                                                                                                                                                                                                           ctctTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGt                                                                                                                                                             >  8:436401/1‑135 (MQ=255)
                                                                                                                                                                                                                           ctctTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGt                                                                                                                                                             <  7:436401/135‑1 (MQ=255)
                                                                                                                                                                                                                             ctTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGa                                                                                                                                                                                                                                           <  8:330383/55‑1 (MQ=255)
                                                                                                                                                                                                                               tcgtcgTTAAAAATCGCTTTTAACAGCGATTCAGTGTCGGGCGCTTCCAGCGAACGGTGAAATTTAGTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGGAATAAACGGACGCTCGt              <  1:67115/274‑1 (MQ=255)
                                                                                                                                                                                                                               tcgtcgTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACg                                                                                                                                                                                                                                        >  7:235077/1‑56 (MQ=255)
                                                                                                                                                                                                                               tcgtcgTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACg                                                                                                                                                                                                                                        <  8:235077/56‑1 (MQ=255)
                                                                                                                                                                                                                               tcgtcgTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTAccc                                                                                                                                            >  8:149344/1‑148 (MQ=255)
                                                                                                                                                                                                                                  tcgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCtt                                                                                                                                                                                                                                                  <  7:538550/43‑1 (MQ=255)
                                                                                                                                                                                                                                  tcgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCtt                                                                                                                                                                                                                                                  >  8:538550/1‑43 (MQ=255)
                                                                                                                                                                                                                                  tcgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATg                                                                                                                                                                                                    <  5:512682/89‑1 (MQ=255)
                                                                                                                                                                                                                                  tcgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATg                                                                                                                                                                                                    >  6:512682/1‑89 (MQ=255)
                                                                                                                                                                                                                                   cgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACAt                                                                                                                                        >  5:229432/1‑148 (MQ=255)
                                                                                                                                                                                                                                   cgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGGAATAAACGGACGCTCGTgcggc         <  4:271043/275‑1 (MQ=255)
                                                                                                                                                                                                                                      tAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACtt                                                                                                                                                                                               <  7:238343/90‑1 (MQ=255)
                                                                                                                                                                                                                                      tAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACtt                                                                                                                                                                                               >  8:238343/1‑90 (MQ=255)
                                                                                                                                                                                                                                      tAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTgg                                                                                                                                     >  7:289284/1‑148 (MQ=255)
                                                                                                                                                                                                                                         gaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCg                                                                                                                                                                              <  7:1040591/103‑1 (MQ=255)
                                                                                                                                                                                                                                         gaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCg                                                                                                                                                                              >  8:1040591/2‑104 (MQ=255)
                                                                                                                                                                                                                                         gaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGa                                                                                                                                  >  7:473340/2‑148 (MQ=255)
                                                                                                                                                                                                                                         aaaTCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGa                                                                                                                                  <  8:289284/148‑1 (MQ=255)
                                                                                                                                                                                                                                          aagcGCTTGTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGAAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGGAATAAACGGACGCTCGTGCGGCGGGATCa  <  2:176817/272‑1 (MQ=255)
                                                                                                                                                                                                                                          aaTCGCGTTTAACAGCGGTTCAGGGACGGGCGCTTCCAGGGAACGGGGAAATTTATTGTTCAGGTTGGCAACGGTGGGATGCACTTTGCCAACCAGCGTATCGGCAGGTGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGGAATAAACGGACGCGCGTGCGGCGGGATCa  <  4:288686/275‑1 (MQ=255)
                                                                                                                                                                                                                                           aTCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAg                                                                                                                                <  7:149344/148‑1 (MQ=255)
                                                                                                                                                                                                                                             cGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACAtt                                                                                                                                       <  8:405539/139‑1 (MQ=255)
                                                                                                                                                                                                                                             cGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACAtt                                                                                                                                       >  7:405539/1‑139 (MQ=255)
                                                                                                                                                                                                                                             cGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCg                                                                                                                              >  8:422702/1‑148 (MQ=255)
                                                                                                                                                                                                                                                ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGaa                                                                                                                           >  8:735992/1‑148 (MQ=255)
                                                                                                                                                                                                                                                 tttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGc                                                                                                                                                                                                   <  6:938116/75‑1 (MQ=255)
                                                                                                                                                                                                                                                 tttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCAt                                                                                                                                                                                         <  6:205769/85‑1 (MQ=255)
                                                                                                                                                                                                                                                 tttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCAt                                                                                                                                                                                         >  5:205769/1‑85 (MQ=255)
                                                                                                                                                                                                                                                 tttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCgg                                                                                                                                                                             <  7:915149/97‑1 (MQ=255)
                                                                                                                                                                                                                                                      cAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGAt                                                                                                                                                                <  8:206328/105‑1 (MQ=255)
                                                                                                                                                                                                                                                      cAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGAt                                                                                                                                                                >  7:206328/1‑105 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                      
TTGGGGATTCACCTGATTCAGCCTTGCCACCGCCTCGCTGGCTTCGCGAGTCTGTCGCGCCGGATAACGGGAAGGCCGGGTATCATTGCCTTCTTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGGAATAAACGGACGCTCGTGCGGCGGGATCA  >  NC_000913/1827570‑1828078

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