Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 699,722 A→G 44.8% intergenic (+115/‑129) PP_t06 → / → PP_23SD tRNA‑Ala/23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451699,7220AG44.8% 68.6 / 155.1 116intergenic (+115/‑129)PP_t06/PP_23SDtRNA‑Ala/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (28/36);  new base G (20/32);  total (48/68)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

GCCCAGACCTACCATTACATGGTTCAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTGGGGTTATATGGTCAAGTGAAGAAGCGCATACGGTGGATGCCTTG  >  AE015451/699494‑699885
                                                                                                                                                                                                                                    |                                                                                                                                                                   
gCCCAGACCTACCATTACATGGTTCAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAGTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTg                                                                                                                                                                  <  2:741213/232‑1 (MQ=25)
     gACCTACCATTACATGGTTCAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAg                                                                                                                                                          <  1:458344/235‑1 (MQ=34)
      aCCTACCATTACATGGTTCAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGtt                                                                                                                                                        <  1:267611/236‑1 (MQ=33)
        cTACCATTACATGGTTCAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTc                                                                                                                                                      >  1:288996/1‑236 (MQ=21)
        cTACCATTACATGGTTCAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTc                                                                                                                                                      <  1:2137/236‑1 (MQ=33)
             aTTACATGGTTCAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGGACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAGAATTCGGAGATGTAATAGATATAGACTGAACGcc                                                                                                                                                            >  2:596168/1‑225 (MQ=255)
              ttACATGGTTCAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGATTGAACAACAGGTTCActgc                                                                                                                                                 >  1:624409/1‑235 (MQ=18)
               tACATGGTTCAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCActgc                                                                                                                                                 <  2:570189/234‑1 (MQ=35)
                aCATGGTGCAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGtt                                                                                                                                                        <  1:454209/226‑1 (MQ=21)
                   tGGTTCAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCtggtgg                                                                                                                                           >  1:490696/1‑236 (MQ=25)
                   tGGTTCAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtg                                                                                                                                            <  1:413587/235‑1 (MQ=33)
                     gTTCAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtg                                                                                                                                            <  2:614303/233‑1 (MQ=33)
                       tccgCCGTCGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAg                                                                                                                                                          <  1:176492/214‑1 (MQ=21)
                       tCAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCa                                                                                                                                       <  1:160370/236‑1 (MQ=33)
                       tCAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCa                                                                                                                                       <  1:592335/236‑1 (MQ=34)
                        cAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAg                                                                                                                                      <  1:795198/236‑1 (MQ=34)
                        cAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAg                                                                                                                                      <  1:14089/236‑1 (MQ=33)
                         aGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAgg                                                                                                                                     >  1:670912/1‑236 (MQ=21)
                         aGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATc                                                                                                                                        >  2:435067/1‑233 (MQ=34)
                         aGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAgg                                                                                                                                     >  1:344227/1‑236 (MQ=34)
                            cGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTa                                                                                                                                  >  1:731624/1‑236 (MQ=21)
                                  aTACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGTATCAGGCTAAGGTaa                                                                                                                            >  1:172914/1‑236 (MQ=18)
                                  aTACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTaa                                                                                                                            >  1:655853/1‑236 (MQ=33)
                                   tACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtgg                                                                                                                                              >  2:459422/1‑217 (MQ=35)
                                           cATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAg                                                                                                                                                          <  1:189677/197‑1 (MQ=21)
                                           cATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAg                                                                                                                                                          >  2:189675/1‑197 (MQ=21)
                                           cATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGTTAAAAtt                                                                                                                        >  1:556719/1‑231 (MQ=21)
                                           cATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGCTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa                                                                                                                   >  1:90947/1‑236 (MQ=18)
                                              agatCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGt                                                                                                                 >  1:206388/4‑235 (MQ=18)
                                              aGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAGCTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTACGGTAAACTTTGTGAg                                                                                                                  >  2:235695/1‑234 (MQ=16)
                                              aGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGtt                                                                                                                >  1:243955/1‑236 (MQ=33)
                                                cTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTc                                                                                                               >  1:418291/1‑235 (MQ=35)
                                                 tCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGtt                                                                                                                <  2:776919/233‑1 (MQ=21)
                                                      tGGGAGAGCGCCGGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAAtt                                                                                                                        >  1:282571/1‑220 (MQ=21)
                                                        ggAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCg                                                                                                         <  2:501508/233‑1 (MQ=25)
                                                           gagCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAaga                                                                                                    >  1:328261/1‑235 (MQ=34)
                                                            agCGCCTGCCTTGCACGCAGGATGTCAGCGGTTCGATCCCGCTTGGCTCCACCTCTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa                                                                                                      <  2:224883/232‑1 (MQ=18)
                                                             gcgcCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGTCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCa                                                                                                 >  1:160068/1‑236 (MQ=17)
                                                             gcgcCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCActgc                                                                                                                                                 <  2:49235/188‑1 (MQ=25)
                                                                cctgcctTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAAc                                                                                               >  2:108936/1‑235 (MQ=18)
                                                                  tgcctTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTa                                                                                                                                  >  2:209243/1‑198 (MQ=21)
                                                                  tgcctTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTa                                                                                                                             <  1:578019/203‑1 (MQ=25)
                                                                  tgcctTGCACACAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCActgct                                                                                                                                                >  2:326213/1‑184 (MQ=25)
                                                                      ttGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACCCCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGaa                                                                                        <  1:199616/236‑1 (MQ=18)
                                                                      ttGCACGCAAGATGTCAGCTGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCt                                                                                                              <  2:266801/214‑1 (MQ=19)
                                                                          ccGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa                                                                                                      <  1:699973/217‑1 (MQ=21)
                                                                           cGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTTCTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTg                                                                                                                                                                  <  2:234311/157‑1 (MQ=25)
                                                                           cGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTGGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAg                                                                                                                                                          >  1:703394/1‑165 (MQ=25)
                                                                            gCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGc                                                                                                  >  2:804838/1‑220 (MQ=25)
                                                                             cAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTaa                                                                                                                                 >  2:187830/1‑188 (MQ=35)
                                                                               ggaggTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCg                                                                               <  1:241465/236‑1 (MQ=18)
                                                                                 aggTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGCATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGTTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCTGCGa                                                                              >  1:45333/1‑235 (MQ=16)
                                                                                     cAGCTGTTCTGTCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGtt                                                                                                                                                        <  2:22154/157‑1 (MQ=255)
                                                                                      aGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACgc                                                                                         >  1:76120/1‑218 (MQ=21)
                                                                                       gCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAAt                                                                            <  2:371954/231‑1 (MQ=34)
                                                                                        cGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATgtcgtc                                                                      <  1:214442/236‑1 (MQ=18)
                                                                                        cGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATgtcgtc                                                                      <  1:118686/236‑1 (MQ=33)
                                                                                        cGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATgtcgtc                                                                      <  1:459424/236‑1 (MQ=34)
                                                                                        cGGTTCGATACCGCTTCGCTCCACCTCTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATCTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtg                                                                                                                                            <  2:429733/166‑1 (MQ=255)
                                                                                         ggTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATCTGTGATAGATATAGGCTGAACACCAGTTTCACTGCtggt                                                                                                                                             >  2:28549/1‑164 (MQ=25)
                                                                                          gTTCGATCCCGCTTTGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCtt                                                                    <  2:176554/236‑1 (MQ=18)
                                                                                            tCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa                                                                                                      <  2:775590/200‑1 (MQ=21)
                                                                                             cGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGt                                                                                             >  2:183556/1‑208 (MQ=21)
                                                                                              gATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTcac                                                                 >  1:526442/1‑235 (MQ=34)
                                                                                                   cGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCg                                                                               <  1:354618/216‑1 (MQ=21)
                                                                                                     cTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCGTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATgtcgtc                                                                      <  2:70535/223‑1 (MQ=21)
                                                                                                        ggCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAg                                                      >  1:80999/1‑236 (MQ=33)
                                                                                                           tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATgtcgtc                                                                      <  2:526440/217‑1 (MQ=34)
                                                                                                            ccaccaCTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATgtcgtc                                                                      <  1:804840/216‑1 (MQ=21)
                                                                                                             caccacTTTCGCCGTACGCAGTTTCTTGTCAGAACATAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTGTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtg                                                                                                                                            <  2:238658/145‑1 (MQ=255)
                                                                                                             caccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCCCGAAGGTTGATTTCTGGCTTTTGTCAGAGCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCActgc                                                                                                                                                 >  2:308092/1‑140 (MQ=255)
                                                                                                               ccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTACAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCg                                                                               <  2:420219/204‑1 (MQ=25)
                                                                                                                  cTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATgtcgt                                                                       <  2:35790/209‑1 (MQ=34)
                                                                                                                       gCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGACAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAg                                                      >  2:135530/1‑221 (MQ=21)
                                                                                                                        ccGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAAtt                                                                                      >  2:282464/1‑188 (MQ=35)
                                                                                                                           tACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATgtcgt                                                                       <  2:425527/200‑1 (MQ=34)
                                                                                                                           gaCGCAGTTACTTGTCAGAACTTAGGAATGAACATTCAGCTCGAATGTTGAGTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa                                                                                                      <  2:622411/168‑1 (MQ=255)
                                                                                                                            aCGCAGTTACTTGTCAGCACTTAGAAATGAACATTCAGCTCGAATGTTGATGTCTGGCTTTTTTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTGGGGtt                                         <  2:20226/229‑1 (MQ=21)
                                                                                                                             cGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTGGGGTTATATGGt                                  <  1:505508/235‑1 (MQ=18)
                                                                                                                              gCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTc                                                                                                          >  1:487687/1‑162 (MQ=35)
                                                                                                                                aGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTGGGGTTATATGGTCAAg                              >  1:419067/1‑236 (MQ=21)
                                                                                                                                 gTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGAtt                                                   <  1:622890/214‑1 (MQ=21)
                                                                                                                                 gTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTGGGGTTATATGGTCAAGt                             <  1:190240/236‑1 (MQ=18)
                                                                                                                                  ttACTTGTCAGAACTTAGAAATGAACATTCAGCTCGCATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGAtt                                                   <  2:348629/213‑1 (MQ=18)
                                                                                                                                   tACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATgtcgtc                                                                      <  2:704385/193‑1 (MQ=18)
                                                                                                                                   tACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTGGGGTTATATGGTCAAGTga                           <  1:666377/236‑1 (MQ=33)
                                                                                                                                    aCTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTGGGGTTATATGGTCAAGTgaa                          <  1:439084/236‑1 (MQ=32)
                                                                                                                                     cTTGTCAGAACTTAGAAATGAACATTCAGCGCGAATGTTGATTGCTGTCTTTTGTCAGATCGTTCTCTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCa                                                                                                 <  2:661363/164‑1 (MQ=255)
                                                                                                                                          cAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGAt                                                    >  1:657791/1‑204 (MQ=18)
                                                                                                                                             aaCTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTGGGGTTATATGGTCAAGTgaagaa                       >  1:54657/1‑230 (MQ=32)
                                                                                                                                                ttAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCGCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTGGGGTTATATGGTCAAGTGAAGAAGCGCATACg              <  1:87660/236‑1 (MQ=17)
                                                                                                                                                ttAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTGGGGGTATATGGTCAAGTGAAGAAGCGCATACg              >  1:147431/1‑236 (MQ=17)
                                                                                                                                                 tAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATgtcgtc                                                                      <  2:275184/179‑1 (MQ=33)
                                                                                                                                                      aTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGaa                                                                                        <  1:748111/156‑1 (MQ=18)
                                                                                                                                                           cATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTg                                              <  2:665483/193‑1 (MQ=33)
                                                                                                                                                            aTTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATATTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTGGGGTTatat                                     >  2:449900/1‑201 (MQ=18)
                                                                                                                                                            aTTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATAAAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTGGGGGTATATGGTCAAGTGAAGAAGCGCATACGGTGGATGCCTTg  >  1:382678/1‑236 (MQ=18)
                                                                                                                                                                aGCTCGAATGTTTATTTCTTGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTGGGGTTATATg                                    <  2:160946/198‑1 (MQ=14)
                                                                                                                                                                aGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTGGGGtt                                         <  2:511185/193‑1 (MQ=14)
                                                                                                                                                                 gCTCGAATGTTGATTTCTTGCTTTTGTTAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTGGGGTTATATg                                    <  1:700703/197‑1 (MQ=14)
                                                                                                                                                                                                                                    |                                                                                                                                                                   
GCCCAGACCTACCATTACATGGTTCAGCCGTAGAATACGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTATAACCAGATTGCTTGGGGTTATATGGTCAAGTGAAGAAGCGCATACGGTGGATGCCTTG  >  AE015451/699494‑699885

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