Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,572,319 T→G E395A (GAA→GCA)  glgB ← 1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,572,3190TG97.5% 239.1 / ‑4.2 80E395A (GAA→GCA) glgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (47/30);  total (49/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGACCGCCCATATCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACGGAAGTCGTCGCG  >  NC_000913/3572079‑3572560
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
aGACCGCCCATATCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGt                                                                                                                                                                                                                                          >  1:286920/1‑250 (MQ=255)
  cccGCCCATATCCTGCGGACGAGAAACGCCGGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCGCTCTTTACGGCTGTGGGCGCGATAAATTAGTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTaaa                                                                                                                                                                                                                                       <  1:37612/250‑1 (MQ=255)
       ccATATCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGCCTCCTCAGCCATTGTCACCGCACCGGAAACATGCTCACCAAGCATACGGTTGGTATTACGCCAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAACCGTGCAATACAGTAAAGagcg                                                                                                                                                                                                                                  >  2:56061/1‑251 (MQ=255)
             ccTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTAcc                                                                                                                                                                                                                             <  2:209138/250‑1 (MQ=255)
              cTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGGAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACcgc                                                                                                                                                                                                                           >  1:194410/1‑250 (MQ=255)
               tGCGGGCGAGAAACGCCAGAGAAATCGGTAGACTACTCAGCAATTGTCACCGCACCGGAAACATGATCACAAAGAATACGGTTGGTATTACGCAAGAATTCAATCGATTCAAGCTTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCACCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGAATCGACGGCAAGCGCATCAATACCAAAACATGGCATCCAGTAAAGcgcgataacga                                                                                                                                                                                                                           >  2:66071/1‑243 (MQ=255)
                gCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACcga                                                                                                                                                                                                                           <  1:331165/249‑1 (MQ=255)
                        gaAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGAGCGTTACCGACGAGGAGGtt                                                                                                                                                                                                                 >  1:351234/1‑251 (MQ=255)
                        gaAACGCCAGGGAAATCGGTAGACTCATCAGCAATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAgct                                                                                                                                                                                                                 >  2:302641/1‑249 (MQ=255)
                           aCGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGTGCGTTACCGACGAGGAAGTTACt                                                                                                                                                                                                              <  1:350086/251‑1 (MQ=255)
                             gCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCCGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGGTACCGACGCGGAAGATACTg                                                                                                                                                                                                             >  2:38370/1‑250 (MQ=255)
                              ccAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCCGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGCCGAGGAAGTTAATGAc                                                                                                                                                                                                           >  2:148648/1‑251 (MQ=255)
                                     aaTCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACgg                                                                                                                                                                                                    >  1:444150/1‑251 (MQ=255)
                                      aTCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAACTTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGAGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGAGAAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGc                                                                                                                                                                                                   >  1:389042/1‑251 (MQ=255)
                                      aTCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGc                                                                                                                                                                                                   >  1:320848/1‑251 (MQ=255)
                                             gacacctcattcATTGTCACCGCACCGTAAACCGGCTCACCAAGAATACTGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGGTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCAGCGCATCGACGCGCAGCGCATCAATTCCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATa                                                                                                                                                                                            <  2:355721/240‑1 (MQ=255)
                                                      gCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGAt                                                                                                                                                                                   <  1:388334/251‑1 (MQ=255)
                                                       ccATTGTCACCTCACCGGAACCCTGCTCCCCAAGAATCCGTTTGGTATTCCGCAAGAATGCAATCTCTTCAAGATTTCCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATc                                                                                                                                                                                  <  2:378275/251‑1 (MQ=255)
                                                         aTTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAACACGAGCAATCCAGTAAAGCGCGTAACCGACGAGGAAGTTACTGACTTCACGGCGCCCATAGTTGTAGATc                                                                                                                                                                                  >  2:391078/1‑249 (MQ=255)
                                                           tGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCTCGTTACCGTCGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGc                                                                                                                                                                               <  1:23244/250‑1 (MQ=255)
                                                              caccgcaccgGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCgtgt                                                                                                                                                                           <  2:444150/251‑1 (MQ=255)
                                                                  gcaccgGCAACCTGCACACAAAGAAAAAGGTTGGTATTACGCAAGAACTCAATCGCTTCAAGATCCTCGAGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTAAGGCTGTAGTCGCGATAAATCATTGACGCCCCCGCATAGACGCCCGGCGCATCAATAACAAAACGCTCAATCCAGTAAAACGCGTTACCGAAGAGGACGTTACTGACTTCACGGCCACCATAGTCTTAGATCAGCGTGttacc                                                                                                                                                                       >  2:449183/1‑248 (MQ=255)
                                                                      cgGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGCCGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTcc                                                                                                                                                                   >  2:297292/1‑251 (MQ=255)
                                                                       ggAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTAAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATACACTCCCCCTCTTTACGGCTGTAGTCGCGATACATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCAAGTAAAGCGCGATACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCt                                                                                                                                                                  >  1:226016/1‑251 (MQ=255)
                                                                        gAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtg                                                                                                                                                                 <  1:89631/251‑1 (MQ=255)
                                                                        gAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGAGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGCCGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtg                                                                                                                                                                 >  1:308076/1‑251 (MQ=255)
                                                                         aaaCCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtg                                                                                                                                                                 >  1:72680/1‑250 (MQ=255)
                                                                                cTCACCCAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGc                                                                                                                                                         >  1:368212/1‑251 (MQ=255)
                                                                                  cACCAAGAATACGGTTGTTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTGACTGACTGCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGcc                                                                                                                                                        <  1:428019/250‑1 (MQ=255)
                                                                                      aaGAATACGGTTGGTATTACTCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCAc                                                                                                                                                   <  1:438956/251‑1 (MQ=255)
                                                                                          aTACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAg                                                                                                                                                                                                                       <  2:273490/179‑1 (MQ=255)
                                                                                          aTACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAg                                                                                                                                                                                                                       >  1:273490/1‑179 (MQ=255)
                                                                                                     aTTACGCAAGAATTCAATCGCGTCAATATGCTCGCACCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCAACGACGCGCAGCGCATAAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCTGATCGCTGTGtt                                                                                                                                    <  2:75056/251‑1 (MQ=255)
                                                                                                         cGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTAGGGATCCACTCCCCCTCTTTACGACTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGAAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCAt                                                                                                                                 >  1:417490/1‑250 (MQ=255)
                                                                                                         cGCAAGAATACAATCGCTTCAAGATTCTCGAGCCCGCCAAACTCGTTAGGGATCCACTACCCCTCTTTAAGGCTGTAGTCGAGATAAATCATTGACGCCCCCGAATCGACGCGCAGCGCATCCATACCAACACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATCGTTGTAGATCTGCGTGTTCAAGTCCTGATGATAGCCTTCACGCGGATCACTGTGTTCAt                                                                                                                                 >  2:405511/1‑250 (MQ=255)
                                                                                                           cAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACa                                                                                                                              <  1:275061/251‑1 (MQ=255)
                                                                                                            aagaatcaaaTCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACa                                                                                                                              >  1:126614/9‑250 (MQ=255)
                                                                                                                 ttCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTCCTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCg                                                                                                                        >  1:297456/1‑251 (MQ=255)
                                                                                                                    aaTCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGCCGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCAAGCGGATCGCTGTGGTCATACAAGATCGTGc                                                                                                                     >  2:273609/1‑251 (MQ=255)
                                                                                                                            cAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATCGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAAt                                                                                                             >  1:267959/1‑251 (MQ=255)
                                                                                                                                 ttCTCGCGCCCGCCAAATTCGATCGGGATCCATTCCACCCCTTTCCGGCAGTAGTCGAGATAAATCATTAACGCCACCGCATCGACGCGCAGCGCAACAATACCAAAACGTGAAACCCAGTAACGCGCGTTACCGACGAGGAAGTTACTGAATTCACGGCGACAATAGATGTAGATCAGCGTGTTCCAGTCCAGATGATAGCCTTCACGCGGATCGCTGTGTTCATACACGTTCGTGCCATCAAATTCgg                                                                                                         >  1:446503/1‑250 (MQ=255)
                                                                                                                                   ctcGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCa                                                                                                       >  1:112131/1‑250 (MQ=255)
                                                                                                                                     cgcgCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCGGTGTTCATACAGGTTCGTGCCATCAAATTCGGCAAgc                                                                                                    <  2:286920/251‑1 (MQ=255)
                                                                                                                                      gcgcCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAgcg                                                                                                   >  1:43987/1‑251 (MQ=255)
                                                                                                                                       cgcccgccAAATTCGTTCGGGATCCACTCCCCATCTTTACGGCTCTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACAAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATCCGGCAAgcgc                                                                                                  >  1:429883/1‑251 (MQ=255)
                                                                                                                                           cgccAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCaa                                                                                                >  1:414840/1‑249 (MQ=255)
                                                                                                                                              caacttcgttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCGAAGtca                                                                                           <  1:119021/247‑1 (MQ=255)
                                                                                                                                                    cgttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACACGAAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTc                                                                                     >  1:271841/1‑251 (MQ=255)
                                                                                                                                                      ttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCCTAGACGCGCAGCGCATCACTACCAAACCGTGCAATCCAGTAACGCGCGTTACCGACGAGGAAGTTACTGAATTCACGGCGAACATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCgg                                                                                   >  2:59313/1‑251 (MQ=255)
                                                                                                                                                        cgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCggg                                                                                  >  1:317384/1‑250 (MQ=255)
                                                                                                                                                                aCTCCCCTTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGAAGCGAAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGCTAGCCTTCACGCGGGTCGCTGCGTTCATACAAGTTCGTGCCATCAAATTCGGCCAGCGCAACGTAATCAGGCGGGACGTGGCCt                                                                         >  2:129296/1‑251 (MQ=255)
                                                                                                                                                                aCTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCACTACCAAAACGTGAAATCCAGTAAAGCTCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCt                                                                         >  1:231424/1‑251 (MQ=255)
                                                                                                                                                                   cccccTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTg                                                                        <  1:57119/249‑1 (MQ=255)
                                                                                                                                                                      ccTCTTTACGGCTGTAGTCTCGATAAATCATTGTCGCCTCCGCATCGACTCGCAGCTCATCAATACCAAAACGGGCAATCCAGTAAAGCGCTTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCATCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCAGCATATGCTGCAAGGGCAAAGTCATCAGTCGGGAAGTGGCCTGGCAcc                                                                   <  1:374501/251‑1 (MQ=255)
                                                                                                                                                                        tctTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCAccc                                                                  >  1:18482/1‑250 (MQ=255)
                                                                                                                                                                           ttACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAATGCGCGTTACCGACGAGGAAGTTACTGACTTCTCGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAAGTCGGCAAGCGCAAAGTCAGCAGTCGGGAAGTGGCCTGGCACCCAg                                                                <  1:453640/249‑1 (MQ=255)
                                                                                                                                                                                gCTGTAGTCGCGATAAATCATTGACGCCACCGCATAGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCATGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACACAGTCgaga                                                          >  1:153184/1‑250 (MQ=255)
                                                                                                                                                                                 ttGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAAt                                                        >  2:237371/2‑251 (MQ=255)
                                                                                                                                                                                   gTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGGGGCCTGGCACCCAGTCGAGAATc                                                       <  2:414840/250‑1 (MQ=255)
                                                                                                                                                                                         gcgATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTc                                                 >  1:299186/1‑250 (MQ=255)
                                                                                                                                                                                                   actGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGACATCAGTCGGGAAGTGGCCTGGCACCACGTCGAGAATCACGTTCAGACCAGCTGc                                      >  2:329617/3‑251 (MQ=255)
                                                                                                                                                                                                   aTTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGc                                      >  1:69401/1‑251 (MQ=255)
                                                                                                                                                                                                     tGACGCCACCTCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTGACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCTAGAATCACGTTCAGACCAGCTGCg                                     <  2:165013/250‑1 (MQ=255)
                                                                                                                                                                                                      gACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCgt                                    <  2:448137/250‑1 (MQ=255)
                                                                                                                                                                                                         gCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCCATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGACGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGc                                >  1:342244/1‑251 (MQ=255)
                                                                                                                                                                                                         cccACCTCATCGACGCGCAGCGCATCAAGTCCAAAACGTGCAATCCATTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTTTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGTAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGTTGCTTGTGc                                <  2:439376/250‑1 (MQ=255)
                                                                                                                                                                                                          ccACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCg                               <  1:66071/251‑1 (MQ=255)
                                                                                                                                                                                                           cACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCgg                              >  1:397451/1‑251 (MQ=255)
                                                                                                                                                                                                               gCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATc                          >  1:172983/1‑251 (MQ=255)
                                                                                                                                                                                                                  acGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCTCGTTCAGACCAGCTGCGTGTGCGGCATCAAt                       <  1:412418/250‑1 (MQ=255)
                                                                                                                                                                                                                   cGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACTCGGATCGCTTTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAAt                       <  1:243530/250‑1 (MQ=255)
                                                                                                                                                                                                                    gACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTCAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATg                      >  1:263080/1‑250 (MQ=255)
                                                                                                                                                                                                                    gACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATg                      >  1:354365/1‑250 (MQ=255)
                                                                                                                                                                                                                    gACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGGCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTTTGCGGCATCAATGa                     <  1:348062/251‑1 (MQ=255)
                                                                                                                                                                                                                        cgcagcgcaTCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATa                 <  1:131032/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aTCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCAGAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACGGAAgt         <  2:84123/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aaaTACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACGGACTTCACGGCGACCATAGTTGTAGATCATCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGGTCATACAAGTGCGTGCCATCAAAGTCGGCAAGCGCAAAGTCAGCAGTCGGGAAGTGGCCTGGCACCGAGTCGAGAATCACGTTCAGACCAGCTGCGTGGGCGGCATCAATTAAATATCTGAAgtc        <  2:132854/249‑1 (MQ=255)
                                                                                                                                                                                                                                    aaaaCCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACAGACTTCACGGCGACCATAGTTGTAGATCGGCGAGTTCCAGTCCTGACGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGACGTGGCCTGGCACCAAGTCGAGAATAACGTTCAGACCAGCTGCGGGTGCGGCATCAATGAAATAACGGAAgacgt      >  2:354695/4‑251 (MQ=255)
                                                                                                                                                                                                                                    aaTACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACGGAAgtcgt      >  1:306983/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aaaaCGTACAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGTCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACGGAAGTCGTcgcg  >  1:162822/1‑249 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
AGACCGCCCATATCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACGGAAGTCGTCGCG  >  NC_000913/3572079‑3572560

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