Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,323 G→A D76D (GAC→GAT araD ← L‑ribulose‑5‑phosphate 4‑epimerase AraD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091366,3230GA99.0% 306.2 / NA 100D76D (GAC→GATaraDL‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (45/54);  minor base T (0/1);  total (45/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTG  >  NC_000913/66080‑66565
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
aCCCTGTTTTCCGAAGGTTTCTACGATGACGTTACCGTTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCggag                                                                                                                                                                                                                                              <  1:397596/250‑1 (MQ=255)
         tcagagggTTTCTACGATGTCTTTTCCGTTTTCCCACACATATTCACCTTTGATTTCTGCGTCGGTCTTTTTGCGGGTGCGGGGAATGGTGCCGTAGAAATAGTCGGCGTAGGTGGTGCCGGTTGCTGGAATCGACTGCCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGAGAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCttt                                                                                                                                                                                                                                     <  2:666789/244‑1 (MQ=255)
           cgaaGGTTTCTACGATGACGTTACCGGTTTCCCCCTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTc                                                                                                                                                                                                                                  <  1:121359/249‑1 (MQ=255)
           cgaaGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTc                                                                                                                                                                                                                                  <  1:219508/249‑1 (MQ=255)
             aaGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATTGAGCAGTCGGTGGGGTGGCGGATCGGGGGAGGGCTTTTTCgc                                                                                                                                                                                                                                >  1:275749/1‑250 (MQ=255)
                gTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCAcc                                                                                                                                                                                                                             <  1:562051/251‑1 (MQ=255)
                   tCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCtt                                                                                                                                                                                                                           <  1:574158/250‑1 (MQ=255)
                     tACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGGCTGGCCCGCCTGCGCCCAGATCGTGGCGTGGCGCGGGTGTGTGTGCACCATGCCGCCAATGGACGGGAATGCCTGATCGAGAAGTCGGTGGGTTGGCGTATCGGAGGAGGGATTTTTCGCCCCTTCa                                                                                                                                                                                                                         >  2:224452/1‑250 (MQ=255)
                      aCGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAAGGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCa                                                                                                                                                                                                                         <  1:251328/249‑1 (MQ=255)
                         aTGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGGCGGTGGGGTGGCGTATCGGAGGGGGGCTTTTTCTCACCTTCAACCAc                                                                                                                                                                                                                    >  2:253672/1‑251 (MQ=255)
                                 aCCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCgg                                                                                                                                                                                                            >  1:173880/1‑251 (MQ=255)
                                    ggTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGGATCGGAGGGGGGCTTTTTCGCACCTTCAACCACTTCACCGGttt                                                                                                                                                                                                         >  1:514228/1‑251 (MQ=255)
                                       ttCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGCATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGAAGTCGGGGAGTTGGCGTATCGGAGGGGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCg                                                                                                                                                                                                       >  1:408807/1‑250 (MQ=255)
                                               cATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGCGGTGGCGTATCGGGGGAGGGCTTTTTCGCACCTTCAACCACTTCACAGGTTTCGATGCTAAcc                                                                                                                                                                                              >  1:651535/1‑251 (MQ=255)
                                                  atTCTCCGTTGATTTCGGCGGCGGTCAGTTGGCGGGTGCAGGGAATGGGGCCGTAGAAATAGTCGGCGTGGGTGGGGCCGGTTGCTGGAATCGACGGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCAc                                                                                                                                                                                            <  2:461571/250‑1 (MQ=255)
                                                       ccgtttattccggCGTCGGGCATTTTGCGGGTTCGGGTAAGGGCTCCGGAGAAATAGCCGGCGTGGGTGGTGCCGGTTGCTGGAACCTACTGGCCCGCCTGCGCCCGGATAGTGGCGGGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTTATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCag                                                                                                                                                                                      <  2:72530/239‑2 (MQ=255)
                                                         gTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCatat                                                                                                                                                                                    <  1:458397/251‑1 (MQ=255)
                                                           tGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGGGTGGGTGTGCACAATGCCGCCAATGGACGGGGATGCCTGATAGAGCAGGCGGGGAGTTGGCGGATCGGAGGGGGGCTTTTTTTCCCCTTCAACCACTTCACCGGTTTCGAGGCTAACCCCGACCATATCg                                                                                                                                                                                  >  2:247820/1‑251 (MQ=255)
                                                               ttCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGGTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCGAATGGGCGGGGATGCCTGATAGAGCAGTCGGAGAGTTGGGGTATCGGAGGAGGGCTTTTTCGCAACTTCAACCACTTCACCGGTTTCGATGCGAACCACGACCCTATCGTCa                                                                                                                                                                               >  2:560252/1‑250 (MQ=255)
                                                                 cTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCg                                                                                                                                                                            >  1:570357/1‑251 (MQ=255)
                                                                       cGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTAGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCCTATCGTCAGCGGGCATg                                                                                                                                                                      >  1:639313/1‑251 (MQ=255)
                                                                                ggggggTGCAGGGAATGGTGCCTTAGAAATAGTCGGCGGGGGGTGTGCTGGTTGCGTGAATCGACGGGCCCGCCTGCGCCCAGATAGGGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTaa                                                                                                                                                             <  2:465789/248‑1 (MQ=255)
                                                                                 gCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCCATGGACTGGAATGCCTGATAGAGCAGTCGGGGAGTTGGCGTATCGGGGGGGGGCTTTTTCCGCCCTTCAACCACTTCACCGGTTTCGATGCTAACCCCGACCATATCGTCCGCGGtactgcatctggaat                                                                                                                                                            >  2:334751/1‑237 (MQ=255)
                                                                                 gCGGGTGCAGGGAAGGGTGCCGAAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACGGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGTCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATTCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTaa                                                                                                                                                             <  1:600551/250‑1 (MQ=255)
                                                                                      gtgaGGGAATGGGGCCGGATAAATAGGCGGGGTGGGTGGGGCGGGTTGCTGGAATCGAGTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTATGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACg                                                                                                                                                       <  2:60906/248‑1 (MQ=255)
                                                                                      ggCAGGGAATGGGGCCGTAGAAATGGCCGGCGTGGGTTGGGCCGGGTGCGGGATCCGACTGGCCCGCCGGCGCCCAGATAGTGGCGTGGCGCGAGTGGGTGTGCTCAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCGTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACg                                                                                                                                                       <  2:403651/250‑1 (MQ=255)
                                                                                       gCAGGGAATGGAGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGc                                                                                                                                                      <  1:97523/251‑1 (MQ=255)
                                                                                        cAGGGAATGGGGCCGTAGAAATAGTCGGCGTGGGTGGGGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGcc                                                                                                                                                     <  1:176328/251‑1 (MQ=255)
                                                                                           ggAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCgg                                                                                                                                                   >  1:458499/1‑250 (MQ=255)
                                                                                            gAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGGATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCCCTTCACCGGGTTCGATGCTAACCACGACCATATCGCCAGCGGGCATGATGCTGTAATCGCCGCCGGa                                                                                                                                                  >  1:6456/1‑250 (MQ=255)
                                                                                                ggTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGGGTGTGCACAAAGCCGCCAATGGACGGGAATGCCTGATAGAGCAGGCGGAGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCCCCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGGCATGATGCTGTAATCGCCGCCGGGAGGtt                                                                                                                                             >  2:675119/1‑251 (MQ=255)
                                                                                                ggTGCCGTAGAAATAGTCGGCGTGGGGGGGGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGtt                                                                                                                                             <  1:362852/251‑1 (MQ=255)
                                                                                                   gCCGTGGAAATCGTCGGGGGGGGTGGGGCCGGTTGGTGGAAGCGGCTGGCCCGCCTGCGCCCATATTGTGGCGGGGCGCGAGTGTGGGTGAACAATGCCGCCAATGGGCGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGGATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTg                                                                                                                                           <  2:169451/250‑1 (MQ=255)
                                                                                                   gCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCCGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTg                                                                                                                                           >  1:502097/1‑250 (MQ=255)
                                                                                                   gCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGa                                                                                                                                          <  1:291300/251‑1 (MQ=255)
                                                                                                      gtaaaagaagcggtggggggtggTGCCGTTGTCTGGATTCGAGTGCACCGCCTGCGCCCAGATAGGGGGGGGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGGCATGATGCTGTAATCGACGCCGGAAGGTTTGATca                                                                                                                                       <  2:35960/234‑1 (MQ=255)
                                                                                                      gTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATCAAGCAGTCGGTGAGTTGGCGTATCGGAGGGGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCCCGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGACGGATTGAta                                                                                                                                        >  2:353550/1‑249 (MQ=255)
                                                                                                           aaTAGTCGGCGTGGGGGGGGCGGGTTGCTGGAATCGCATGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGGGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAg                                                                                                                                  <  2:235005/251‑1 (MQ=255)
                                                                                                            ataggcgggggggggggggcgggTTGCGGGAATCGATGGCCCCGCCTGCGCCCAGATGGTGGCGTGGCGCGGTGGTGTGTGCCCAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGa                                                                                                                                 <  2:703321/230‑1 (MQ=255)
                                                                                                              aGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCt                                                                                                                                                                 >  1:54133/1‑217 (MQ=255)
                                                                                                              aGTCGGCGTGGGTGGGGCCGGTTGCTGGAATCGCCTGGCCCGCCTGCGCCCAGATAGGGGCGGGGCGCGAGGGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCTCCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCt                                                                                                                                                                 <  2:54133/217‑1 (MQ=255)
                                                                                                               gttgggttgggtggTGCCGGTTTCGTGAACCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGAcgc                                                                                                                              <  2:606641/244‑1 (MQ=255)
                                                                                                               gTCGGCGGGGGGGGGGCGGGTTGCGGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAACGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGAcgc                                                                                                                              <  2:37232/251‑1 (MQ=255)
                                                                                                                 cGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGAcgccg                                                                                                                            <  1:675119/251‑1 (MQ=255)
                                                                                                                    cGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCCCCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTTATCACAAAGACCCCgcgc                                                                                                                         >  2:5013/1‑251 (MQ=255)
                                                                                                                            gtgCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATc                                                                                                                 >  1:12957/1‑251 (MQ=255)
                                                                                                                            gtgCCGGTTGCTGGAACCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATc                                                                                                                 <  2:115399/251‑1 (MQ=255)
                                                                                                                            gtgCAGGTTGCGGGAACCGACTGCCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGATTGTGTGTGCACAATGCTGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGTGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCGCGACCATATCGTCAGCGGTCATGATGCTGTAAGCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATc                                                                                                                 <  2:342363/251‑1 (MQ=255)
                                                                                                                            gggCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATc                                                                                                                 <  1:225507/249‑1 (MQ=255)
                                                                                                                                   ttGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGc                                                                                                                                                                             >  1:75956/1‑184 (MQ=255)
                                                                                                                                     gCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCt                                                                                                        >  1:225269/1‑251 (MQ=255)
                                                                                                                                        gggaatggacgggcccGCCTGGGCCCAGAGAGGGGCGTGGCGCTATTGGTTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATCGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTGCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAAGCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGa                                                                                                      <  2:432448/238‑1 (MQ=255)
                                                                                                                                          gAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACACTTCCGCCCATGGACGGGAATGCCTGAAAGAGCAGTCGGGGAGTTGGCGTATCGGAGGGGGGCTTTTTCGCACCTTCAACCACTTCCCCGGTTTCGATGCTAACCACGACCATATCGGCAGCGGGCCTGATGCTGTAATCGCCGCCGGAATGTGTGATCACAAAGACGCCGCGCTccccaaccaccgcgctgacg                                                                                                    >  2:382460/1‑232 (MQ=255)
                                                                                                                                          gAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACg                                                                                                    <  1:253672/251‑1 (MQ=255)
                                                                                                                                              cGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGGATCGGAGGGGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGCCGCCGGAAGGTTTGATCAAAAAGACGCCGCGCTCGCGATCAACGCAGCTGACGTTGc                                                                                                >  1:447001/1‑251 (MQ=255)
                                                                                                                                                 ctggcccGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGcccc                                                                                             >  1:486273/5‑251 (MQ=255)
                                                                                                                                                   gtcccGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGTCGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCAt                                                                                           <  1:687683/249‑1 (MQ=255)
                                                                                                                                                   gccccgcttgccccctggtaggggggtggcgGGAGTGGGTGTGCAAAATGCCGACAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCGTCAACCACTTCACCGGTTTCGATGCTAACCACGTCCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCAt                                                                                           <  2:570357/226‑1 (MQ=255)
                                                                                                                                                     cccGCCTGCGCCCAGAGAGGGGCTTGGAGCGGGTGGGGGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCAtgt                                                                                         <  2:189197/251‑1 (MQ=255)
                                                                                                                                                                         ggagtggcgCGAGTGTGTGTGCTCAATGCAGCAAATGGACTGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGGATCGGAGGAGGGCTTTTTCGCTCCTTCAACCACTTCTCCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGGCATGATGCTGTAATCGACGACGTAAGGTGTGAGCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGtt                                                                     <  1:683929/248‑1 (MQ=255)
                                                                                                                                                                          gggtgggggggtgtgtgtgggtacaAGGCCGTCAATGGGCGGGGATGCCTGATATAGCAGTCGGTTAGTTTGCGTATCGGAGGAGTGCTTTTTCGCTCCTTCAACCTCTTCACCGGTGTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTatgt                                                                      <  2:490986/228‑1 (MQ=255)
                                                                                                                                                                               ggcgcgtgggtttgtGCACATTGCCGCCAATGGACGGGAATGCCTGAGAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCGCAAAGACGCCGCGCTCGCGATCAACTGCGGTGACGTTGCCCCATGTGAGCGTGACCAGGttattgtttggc                                                                <  2:651535/237‑11 (MQ=255)
                                                                                                                                                                                  cgAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGCGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGccc                                                                                              >  1:545638/1‑217 (MQ=255)
                                                                                                                                                                                   cgggtgtgtttgctgaaTGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTGTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAGGCTGTAATCGTCGCCGTAAGGTTTGATCACAAAGACGCCGCGCTCGCGTTCAACGGCGCTGACGTCGccc                                                                                              <  2:545638/202‑1 (MQ=255)
                                                                                                                                                                                      gtgtgtgtGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGTCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGGTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCa                                                         <  2:514228/244‑1 (MQ=255)
                                                                                                                                                                                       gtggtggtgcaaaTGCCGCAAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGATCGTGACCAGGTTATGTTTTGGCAGCGCCAg                                                        <  2:14487/240‑1 (MQ=255)
                                                                                                                                                                                        gtgtgtGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAACGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAgg                                                       >  2:264679/6‑251 (MQ=255)
                                                                                                                                                                                         tgtgtGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGt                                                      >  1:393852/5‑251 (MQ=255)
                                                                                                                                                                                          gtgtGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGaa                                                                                                                                                 >  1:434401/4‑159 (MQ=255)
                                                                                                                                                                                          gtgggcccaATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAGGCTGTAATCGACGCCGGaa                                                                                                                                                 <  2:434401/152‑1 (MQ=255)
                                                                                                                                                                                           tgtGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCCCCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGCCGCCGGAAGGTTTGCTCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTg                                                    >  2:692805/3‑251 (MQ=255)
                                                                                                                                                                                              gCACTATGCCGCCAATGGACGGGAATGTCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGTAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGcc                                                 <  2:476957/251‑2 (MQ=255)
                                                                                                                                                                                                  aaTGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGCCGCCGCGCTCGCGATAAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCctct                                              >  2:396057/1‑250 (MQ=255)
                                                                                                                                                                                                       cgccAATGGAAGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAAttcc                                        >  1:605351/1‑248 (MQ=255)
                                                                                                                                                                                                         ccaatgggcgggaagtccttaTAGAGCAGGCGGTGAGTGGTCGGTTCGGAGGTGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCCAACCACGACCATATCGTCAGTGGTCATGATGCTGTAATCGACGCAGGAAGGTTTGAGCACAAAGACGCCGCGATCGCGTTCATCGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAAttcgtg                                      <  2:51207/231‑6 (MQ=255)
                                                                                                                                                                                                          cAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGAGATCAACGGCGCTGACGTTGCCCCATGGGAGCGGGACCAGGTTATGATTTGGCAGCGCCAGGTAGGCCTCTCATACCTgg                                     >  2:654737/1‑251 (MQ=255)
                                                                                                                                                                                                            aTGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCg                                   <  2:12957/251‑1 (MQ=255)
                                                                                                                                                                                                             tGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTGGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAACCGACGCCGGAAGGTTTGATCACAAGGACGCCGCGCTCGCGATCAACGGCGCTGACGGTGCCCCATGTGAGCGTGGCCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTATTACCTggagt                                  >  2:203741/1‑248 (MQ=255)
                                                                                                                                                                                                              ggacgggAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGtt                                 >  1:463710/5‑251 (MQ=255)
                                                                                                                                                                                                                 cgggAATGCCTGATAGAGCAGTCGGGGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGCCGTTGCCCCATGTGAGCTTGCCCAGGTTATGATTTGTCAGCGCCAGGTTGGCCGCTAATACCTCGCGTTTg                               >  2:290081/2‑250 (MQ=255)
                                                                                                                                                                                                                    gAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTg                               >  1:642198/1‑247 (MQ=255)
                                                                                                                                                                                                                     aaTGCCTGGTAGAGCAGTCGGTGCGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTCACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACACAGACGCCGCGCTCGCGATCAACGGCGATGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTCATACCTGGCGTTTg                               >  2:31632/1‑246 (MQ=255)
                                                                                                                                                                                                                      aTGCCTGATAGAGCAGTCGGTGAGTGGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGAtct                         <  2:458499/251‑1 (MQ=255)
                                                                                                                                                                                                                         ccTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGAtcttct                      >  1:451433/1‑251 (MQ=255)
                                                                                                                                                                                                                           tGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTa                     >  1:550505/1‑250 (MQ=255)
                                                                                                                                                                                                                            gATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCTCCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTGTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTa                     <  1:274872/249‑1 (MQ=255)
                                                                                                                                                                                                                            gATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTa                     <  1:471535/249‑1 (MQ=255)
                                                                                                                                                                                                                            gATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTa                     <  1:296800/249‑1 (MQ=255)
                                                                                                                                                                                                                            gATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGGTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTaa                    <  1:32353/250‑1 (MQ=255)
                                                                                                                                                                                                                               agagCAGTCGGTGAGTTGGCGTAGCGGAGGAGGTCTTTTTCGCTCCTTCAACCACTTCGCCGGTTTCGATGCTAACCACGACCATATTGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTTACGTTTCCCCATGTGTGCGTGACCTGGTTATGTGTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATg                <  1:94045/251‑1 (MQ=255)
                                                                                                                                                                                                                                 agcagtCGGTGAGTTGGCGTAACGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGtt              <  1:11960/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cagtcGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTg             >  1:132621/5‑250 (MQ=255)
                                                                                                                                                                                                                                   cagtcGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGa            <  1:610535/247‑1 (MQ=255)
                                                                                                                                                                                                                                        ggTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATCTCGTCAGCGGTCATGATGCTGTAATCGACGCCGAAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCCGCGCCCGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCt       >  1:671791/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ggTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCt       >  1:577158/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ggTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCt       >  1:8182/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gTGAGGTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCTCTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTGTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACtactg      <  2:225871/251‑5 (MQ=255)
                                                                                                                                                                                                                                           gAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCg    <  1:671725/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gAGTTGGCGTATCGGAGGAGGGCTTTTCCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACGCCTTc     <  1:264679/250‑1 (MQ=255)
                                                                                                                                                                                                                                             gatGGCGTATCGGCGGAGGGCTTTTTCGCAACTTCAACCACTTCACCGGTTTCGATGATAACCACGACCATATCGACAGCGGTCATGCTGCTGTAATCGACGCAGGAAGGTTTGATCACAAAGACGCCGCGATCGCGCTCAACGGCGCTGACGTTGCCCCATGCGCGCGTGACCAGGTTATGTTTTGGAAGCGCCAGGTTGGCATCTAATACCTGGCGTTTGAGATCTTCTAACCTGTTGACTCCTTCGtt  >  2:467564/3‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
ACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTG  >  NC_000913/66080‑66565

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