Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 909,181 A→G 13.4% A38A (GCT→GCC ltaE ← L‑allo‑threonine aldolase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913909,1810AG13.4% 174.6 / 2.0 82A38A (GCT→GCCltaEL‑allo‑threonine aldolase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (35/33);  new base G (2/11);  total (37/45)
Fisher's exact test for biased strand distribution p-value = 3.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAA  >  NC_000913/908943‑909415
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               
tCGCCATTTTATCCAGGGGTAGCGGGCCGTCGGCCGCCGCGTCTATGGGTTGCGGTTGAATCCTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCaga                                                                                                                                                                                                                                                 <  1:365336/239‑1 (MQ=255)
         tacccAGCGGTAGCGTGCCGTCGGCAGCCGCTTCTATGGGTTGCGGTTGAATCCTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACg                                                                                                                                                                                                                                       <  1:1177956/236‑1 (MQ=255)
              aGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCCGCGCAGGCCCGCAGAGCATTCACGGTAgg                                                                                                                                                                                                                                  >  2:1290367/1‑239 (MQ=255)
                cGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGAAggg                                                                                                                                                                                                                                 >  1:846310/1‑238 (MQ=255)
                cGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTGCTTTACCGGAAAGCTCTGCTGCGTAGTCCTCCATAGCCTTAACGGTAtgg                                                                                                                                                                                                                                 >  1:482623/1‑238 (MQ=255)
                    aGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGgtcgt                                                                                                                                                                                                                             <  1:310660/238‑1 (MQ=255)
                          ccGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGGTATGCGCGGCCTGACCGACAATATACTCTTCGCCCCGTTCGCAGGGACTGGGCATAGCGACCAGGGTGGCCTGAGGGCCGGTCGGCAGAAAAATGGCGGCTGCTTTACCGGGAAACTGCTCGGCGGAATACTGCGGAGCATTAAAAGTAGGGGCGTCCCCgga                                                                                                                                                                                                                      >  2:979980/1‑237 (MQ=255)
                               ggCAGCCGCGTCTAGGGGTTGCGGTTGAATCCTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAAcgt                                                                                                                                                                                                                 <  1:1271404/239‑1 (MQ=255)
                                             aagagaccGGTTGAATCCTGCCCAGGCCCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCggg                                                                                                                                                                                                   <  1:183095/232‑1 (MQ=255)
                                                          aaTACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGCTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAACGACTTGCTGCGTAGTCCTGCCGAGCCGTACCGGGCGGGCCGTCCCCGCAAACGCCGCCCCCAACCGGGGCGGcccccctc                                                                                                                                                                                       >  2:944969/1‑232 (MQ=255)
                                                                 cccAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGa                                                                                                                                                                                >  1:808679/1‑238 (MQ=255)
                                                                 cccAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGGCGTCTCCGTAAACGTCGTTCCCAACCGCGGCGGCCATCATCCGTTCggg                                                                                                                                                                               >  1:1294553/1‑237 (MQ=255)
                                                                 cccAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGGCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAg                                                                                                                                                                               >  1:29171/1‑239 (MQ=255)
                                                                 cccAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGAATTAACCGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCCTCACCGCTTCggc                                                                                                                                                                               >  1:895268/1‑237 (MQ=255)
                                                                 cccAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGCTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTTCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACCTCGTCCCCAACCGGGGCGGCCATCCTCGCTTCGa                                                                                                                                                                                >  1:124550/1‑238 (MQ=255)
                                                                      cACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCCTCTTTACCGGAACAGTTCGCTGCGGGGGCCAGCAGGGACTTAACGGGAGAGGTGTACCCGCAAACGTCGTCCCCAACCGGGGGGGCCGCCATCCCCTCggagaggg                                                                                                                                                                          >  2:491670/1‑232 (MQ=255)
                                                                           ccgcGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGggcg                                                                                                                                                                     >  1:1201433/1‑238 (MQ=255)
                                                                           ccgcGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGGCCTGCAAGGCATTAACGGTAGGGTCGTTTCCGTAAACGGCGTCCCCAACCCGGGCGGCCATCATCGCTTCGAGCATGGcccg                                                                                                                                                                     >  1:867370/1‑238 (MQ=255)
                                                                           ccgcGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGACTTATTTATCGGGAAGCTCTGCTGCGTAGGCCTTCAGAGAATTAACGGTAGGGTCGTCTCCGTAACCGGCGTCCCCAACCGGGGCGGGCATCAACCCTTGAGTCATGGcgcg                                                                                                                                                                     >  2:244342/1‑238 (MQ=255)
                                                                                ccACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCCTCATCGCTTCGGGCATGGCGCGGCTTgg                                                                                                                                                               >  1:351709/1‑239 (MQ=255)
                                                                                 cACCGGCTTCAAACAGATAGTTACGCGCGGCCTGACCGACAATATACCCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGGGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTgg                                                                                                                                                               <  1:742001/238‑1 (MQ=255)
                                                                                      gttgcAAACAGTTTGTTAGGGGCGGCCTGGCCGACAATGTCCCCGTCGCCGCGGTCGCAGTGACTGGGCAGGGCGACCAGGTTGGCCTGAGTGCCGGGAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGt                                                                                                                                                         <  2:29171/235‑1 (MQ=255)
                                                                                                       aTGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGAAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGCCGGCCATCATCGCTTCTGGCATGGCGCGGCTTCGCCGGGGAACCGTCTCACTGctgc                                                                                                                                        >  1:865497/1‑236 (MQ=255)
                                                                                                        tGCGCGGCATGGCCGAAAATAGAACTTTCGCCGGGTTCGGATGGCCTTGGCAGGGGGACCAGGTTGGCCTGAGTGGCGGTAGGAAGAAAAAGGGCGGCTCCTTTACCGGAAAGCTCTGATGCCTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTaa                                                                                                                                       <  2:839692/239‑1 (MQ=255)
                                                                                                           gcgGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCCTCATCGCTTCGAGCATGGCGCGGCTTGGGCGGGGAACGGGCTCACTGCGTaaa                                                                                                                                      >  1:595113/1‑237 (MQ=255)
                                                                                                           gcgGACTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAaat                                                                                                                                     <  1:1241425/238‑1 (MQ=255)
                                                                                                             ggCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTACCGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGGAACGGTATCACTGCGGAaatcaa                                                                                                                                  >  1:219839/1‑239 (MQ=255)
                                                                                                               ccTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGGAACGGGATCACTGCGTAaataaatc                                                                                                                                >  1:408909/1‑239 (MQ=255)
                                                                                                                cTGACCGAAAATTTACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAaatcaatc                                                                                                                                <  1:16556/238‑1 (MQ=255)
                                                                                                                    ccGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGGCCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATgg                                                                                                                            <  1:244342/238‑1 (MQ=255)
                                                                                                                     cgaacatatGCTCTTCGCCGCGTTCGCAGGGACTGAGCAGAGCGACCATGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCCCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGACATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCa                                                                                                                          <  1:984687/234‑1 (MQ=255)
                                                                                                                     cGACAATATACCCTTCGCCGGGTTGGCAGTGACTGAGCCGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCa                                                                                                                          <  1:201122/239‑1 (MQ=255)
                                                                                                                     cGACAATATACCCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCa                                                                                                                          <  1:420921/239‑1 (MQ=255)
                                                                                                                     cGAAAATATTCTCTTCGCCGCGTTCGCAGTGACTGGGCAGGGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCTATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCa                                                                                                                          <  1:491670/239‑1 (MQ=255)
                                                                                                                       aCAATATACTCCTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGGGCATGGCGCGGCTTTGTCGGGGAACGGTATCACTGCGTAAATAAATCATGGCATg                                                                                                                        >  1:335772/1‑239 (MQ=255)
                                                                                                                        cAATATCCTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCAt                                                                                                                         <  1:166993/237‑1 (MQ=255)
                                                                                                                        cAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGc                                                                                                                           <  1:778101/235‑1 (MQ=255)
                                                                                                                         aaTATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATAACTGCGTAAATCAATCATGGCATGTc                                                                                                                     <  1:888476/239‑1 (MQ=255)
                                                                                                                                  ttCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTCCGAGCATGGCGCGGCCTGGTCGGGGAACGGGTAAACTGCCGAAATCAATCAGGGCATGCTCCtttttg                                                                                                              >  2:888476/1‑237 (MQ=255)
                                                                                                                                  ttCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGCCGTCCCCAACCGGGGCGGCCATCATCGCTTCGCGCATGGCGCGGCCTGGTCGGGGTACGGGATAACTGCGGAAATCCACAATGGCGCTTCCTtattatt                                                                                                            >  2:166993/1‑238 (MQ=255)
                                                                                                                                  ttCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGATCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCCTCATCGCTTCGGGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATAAATCATGGCATGACCTTATTATg                                                                                                            >  1:825855/1‑239 (MQ=255)
                                                                                                                                   tCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGa                                                                                                           >  1:200960/1‑239 (MQ=255)
                                                                                                                                    cgccgcGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGGATCACTGCGTAAATCAATCATGGCAGGTCCTTATTATGAc                                                                                                          >  1:721596/1‑239 (MQ=255)
                                                                                                                                       cgggTTGGCAGTGACTGGTCAGGGCGACCAGGTGGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGGCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACggg                                                                                                       <  2:1201433/236‑1 (MQ=255)
                                                                                                                                                        aGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGGGCATGGCGCGGGTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCCTGTCCTTATTTTGACGGGAACTGCCACCCctttt                                                                                      >  1:605429/1‑238 (MQ=255)
                                                                                                                                                             agCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGAATTAACGGTAGGGTCGTCTCCGTAAACGGCGTCCCCAACCCGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCCTGGTCGGGGAACGGTGTCACACGCTAACTCCATCATGGCATGTCCTTATCTTGACGGGCGATGCCACCCCTTTTAc                                                                                    >  2:1039398/1‑235 (MQ=255)
                                                                                                                                                                   tatgTTGGCCTGAGGGGCGGTGGGTAGAAAAATTGTGGTGTCTTTACCGGAAAGCTCTGCTGCGTGGTCCTGAAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCTCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAgt                                                                         <  2:1128057/236‑1 (MQ=255)
                                                                                                                                                                              tggggCCGGTGGGCATAAAAGTGGCGGCTTCTTTGCCGTAAAGCCTTGTTGCGTGTGCCTGCAGGGCATTAACGGAAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCAGCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGcc                                                              <  2:167372/235‑1 (MQ=255)
                                                                                                                                                                               aGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGgc                                                                                                                                                                        >  1:1062335/1‑135 (MQ=255)
                                                                                                                                                                               aGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGgc                                                                                                                                                                        <  2:1062335/135‑1 (MQ=255)
                                                                                                                                                                               aGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCCATCATGGCATGACCTTCTTATTACGGGGAATTGCACCCTTTTTAACTTAACCAGTTTGTTTTTGc                                                               >  1:1128057/1‑237 (MQ=255)
                                                                                                                                                                                  gCCGGGGGGGAGAAAAGTGGGGCTGTTTTTCCCTGAACGCTTTGCTGCGTGGGCCCGGAGACGATTAACGGTTGGGTCGTCTCCGTAAACGTTGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGt                                                           <  2:904505/238‑1 (MQ=255)
                                                                                                                                                                                      gtgggCCGAAAAATGGGGGTTTCTTTTCCGGAAATCCCGGGTGGGTCGTCCCGACGCGCATTAACGGTAGGGTCGGCCCCGTAAACGGCCTCCCCAACCGGGGCGGCCATCCTCCGTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGAt                                                      <  2:1294553/236‑1 (MQ=255)
                                                                                                                                                                                       tAGGGAGAAAAATGGCGGCTTCCTTACCGGAAAGCTGCGCCGCGTAGTGCTGAAGAGGATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGAt                                                      <  1:946728/238‑1 (MQ=255)
                                                                                                                                                                                       tAGGCATGAAAATGGCGGGTTTCTTACCGGAACGCTCTGCTGCGTAGTCCTGCAGAGCTTTAACGGTATGGTCGTCTCCGTAAACGTCCTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGAt                                                      <  1:649434/238‑1 (MQ=255)
                                                                                                                                                                                        aGGCGGAAAAACGGCGGTTTTGTTACCGACACTGTCTGCTGCGTATGCCTGAAGGGCATTAACGGGAGGGTCGTCTCCGTAACCGTCGTCCCCAACAGGGGCGGCCATCATCGCTGCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCa                                                    <  2:1310624/239‑1 (MQ=255)
                                                                                                                                                                                         ggCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCCACCGGGTCGGCCATCATCGCTTGGAGCATGGCGCGGCTTGGTCGGGGACCGGGATCACACCGGTAATACACCACTGCATGTACTTATTTGGACGGGAGATTGCCCCCCTTTTACCTTCACCCAGTTTTTTTTcccagagccaga                                                      >  2:1063272/1‑226 (MQ=255)
                                                                                                                                                                                         ggCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGt                                                                                                                                                              <  2:123747/135‑1 (MQ=255)
                                                                                                                                                                                         ggCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGt                                                                                                                                                              >  1:123747/1‑135 (MQ=255)
                                                                                                                                                                                         ggCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCCCTGCGTAAATAAATCATgg                                                                                                                            >  1:801135/1‑169 (MQ=255)
                                                                                                                                                                                         ggCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCCTAAATCAATCATGGCCGGGCCTTATTTTTGCGGGGAAAGCACCCCCTTTTACCCTTACCCAGTTGTTTTTGCCaagtgggacac                                                   >  2:1054167/1‑229 (MQ=255)
                                                                                                                                                                                              aaaaaTGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGa                                                                                                           <  1:1014207/180‑1 (MQ=255)
                                                                                                                                                                                                     gcggcTTCTTTACCGGAAAGCTCTGCTGCGTAGGCCTGCAGGGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACcgc                                        <  1:1039398/238‑1 (MQ=255)
                                                                                                                                                                                                      cggcTTCTGTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCTCCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACcgc                                        <  1:834060/237‑1 (MQ=255)
                                                                                                                                                                                                         tttcttTCCCGGAAATCTTTTCTGCGTAGTCCTGCAGGGCATTAACGGTAGGGTCGGTGCCGTAAACGGCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTcc                                    <  2:351709/237‑1 (MQ=255)
                                                                                                                                                                                                         cttcttTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTcc                                    <  1:213481/238‑1 (MQ=255)
                                                                                                                                                                                                                       aaGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCCCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTgcgc                      <  1:26493/239‑1 (MQ=255)
                                                                                                                                                                                                                       aaGCTCTGCTGCGTAGGCCTGCAGGGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCATTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTgcgc                      <  1:1290367/239‑1 (MQ=255)
                                                                                                                                                                                                                        atgCTCTGCCGCGGTGTCCCGCAGGGCTTTACACGTTAGGACGTCTGCGTAAACGTGCTCCCCAACCGGGGCGGCCATCCTCGCTTCGAGCATGGCGGGGCTTGGTCGGGTAACGGTATCACGGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCATCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTgcgc                      <  2:482623/237‑1 (MQ=255)
                                                                                                                                                                                                                        aatctctGCTGCTGAGCCCTGCAGGGCATTAACGGTAGGTTCGTCCCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTgcgc                      <  1:483050/236‑1 (MQ=255)
                                                                                                                                                                                                                         acctctGCCTCGTCGTCCTGCAGAGCATTGACGTTAGGGTCGTTTCCGTAAGCGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGGGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTgcgc                      <  2:895268/236‑1 (MQ=255)
                                                                                                                                                                                                                         aGCCCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTgcgc                      <  1:380338/238‑1 (MQ=255)
                                                                                                                                                                                                                          actcttgttcGTTGTCCTGGAGGGGATTAACGGGAGGGTCGTTGCCGGAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGGGCGGTTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTgcgcgc                    <  2:808679/230‑1 (MQ=255)
                                                                                                                                                                                                                           ctctGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGCTTGTTTTCGCCCGTTCGATCACTTCATCACCTCGTCCCCCTATGATTGCGCgca                   >  1:849134/1‑239 (MQ=255)
                                                                                                                                                                                                                            tctGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTgcgcgc                    <  1:1054167/237‑1 (MQ=255)
                                                                                                                                                                                                                            tctGCTGCGTAGTCCCGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCgca                   <  1:648363/238‑1 (MQ=255)
                                                                                                                                                                                                                                  tcGTTGTCCTGCAAAGCATTAACGGTAGGGTCGTCTCCGGAAACGCCGTCCCCAACCGGGGCGGCCAGCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCAtata            <  1:708076/238‑1 (MQ=255)
                                                                                                                                                                                                                                  gcgttttCCTGCAGGGCTTTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCAtata            <  2:721596/233‑1 (MQ=255)
                                                                                                                                                                                                                                  gcGTTGGCCGGCAGAGCATTAACGGTAGGGTCGTCCCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCAtata            <  1:1145622/239‑1 (MQ=255)
                                                                                                                                                                                                                                     tAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGGGCGGCATATAAAg         >  1:301316/1‑239 (MQ=255)
                                                                                                                                                                                                                                        tCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCCCCGCGTCCGCTGATGATTGTGCCCAGCAATTACAg         >  1:866057/1‑236 (MQ=255)
                                                                                                                                                                                                                                          cTGCAGGGCATTAACGGTTGGGTCGTCTCCGTAAACGGCGGCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGa    <  1:979438/239‑1 (MQ=255)
                                                                                                                                                                                                                                             cAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGGCCGCAGCATATACAGGCTGaaa  >  1:811633/1‑238 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               
TCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGTCTCCGTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAA  >  NC_000913/908943‑909415

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