Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,471,964 T→G 17.9% T145P (ACT→CCT)  pyrB ← aspartate carbamoyltransferase, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,471,9640TG17.9% 177.9 / 2.8 83T145P (ACT→CCT) pyrBaspartate carbamoyltransferase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (32/33);  new base G (12/6);  total (44/39)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTTCTTCAATAGAGCTGTGCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGC  >  NC_000913/4471740‑4472192
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                    
cACTTCTTCAATAGAGCTGTGCAGACTCCATGCGCTCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGGGCGGCCAAATTTCAGGCCACCAACCATTGGGACGTGG‑GAGTGTCCAGAACCCCCCGGGGTTCCCGGAAAAGGGaaa                                                                                                                                                                                                                                  >  2:709293/1‑228 (MQ=255)
         atatgggttgtgCAACCCCCATGCCATCCCTTTTTCATCGAGCATAGCCAGAATGGATTGCGGCTTTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTGCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCACCTGGGTTTCCTGAATAGTGAATAAGTCCAGc                                                                                                                                                                                                                         <  2:412575/221‑1 (MQ=255)
           tAGAGCTGTGCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGACGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCCCCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATTGGGAAAAAGTCCAGCaa                                                                                                                                                                                                                       >  2:506563/1‑229 (MQ=255)
           tAGAGCTGTGCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTACCGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGGGATTGTCCCG‑ACGCCCCCGGGTTTCCGGAATAGTGAATAAATCCAGCaa                                                                                                                                                                                                                       >  2:272830/1‑229 (MQ=255)
           tAGAGCTGTGCAGACTCCATGCGATCCCTTTTGCATCGAGCATATCCAGAATGGATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGc                                                                                                                                                                                                                         <  1:272830/227‑1 (MQ=255)
                cTGTGCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATATGGAATAAGTCCAGCAAGGttt                                                                                                                                                                                                                  >  1:375657/1‑229 (MQ=255)
                cTGTGCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGGAAATAAGCCCAGCAAGGttt                                                                                                                                                                                                                  >  1:384206/1‑229 (MQ=255)
                cTGTGCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGGGCGGCCATATTTCAGGTCACCACCCATTGCGCCGTGGGGATTGTCCAG‑ACGCCCCTGGGGTTCCTGAATAGGGAATAAGTCCAGAAAgtgtt                                                                                                                                                                                                                  >  2:608796/1‑225 (MQ=255)
                cTGTGCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGAGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACTATTGCGACGTGGGGCTTGTCCGA‑ATGCCCCTGGGGTTCCTGAATATGGAATAAGTCCAGCAAGGatt                                                                                                                                                                                                                  >  1:432505/1‑229 (MQ=255)
                   tgCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCg                                                                                                                                                                                                               >  2:78004/1‑229 (MQ=255)
                   tgCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTCCCTGAATAGGGATTAAGTCAAGCAAGGTTTGct                                                                                                                                                                                                               >  1:458015/1‑228 (MQ=255)
                   tgCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGCCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGAAAGGTTTGCg                                                                                                                                                                                                               >  1:520069/1‑229 (MQ=255)
                    ccAGACTCCATGTGATCCCTTGTTCATCGAGCACATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGt                                                                                                                                                                                                              <  1:333770/228‑1 (MQ=255)
                        aCTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGGGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAAAAGTGAATAAGTCCAGAAGAGTTTGGGCCggg                                                                                                                                                                                                          >  2:177937/1‑228 (MQ=255)
                         cTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGGCCGGCGCGGTGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGGCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCgg                                                                                                                                                                                                           <  1:550402/227‑1 (MQ=255)
                         cTCCAAGCGATCCCTTTTACATCGGGCCTATCCAGAATGTATTGCGGCATTTCCAGCGGATCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGGGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCgg                                                                                                                                                                                                           <  1:2614/227‑1 (MQ=255)
                             aTGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATgtt                                                                                                                                                                                                      >  2:88268/1‑228 (MQ=255)
                                  aTCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCCTATTTCAGTTCACCACCCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCAAGCACGGTTTGCGCGGGATGtgggttgg                                                                                                                                                                                                >  2:159592/1‑229 (MQ=255)
                                   tCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATgttggtt                                                                                                                                                                                                  <  1:88268/226‑1 (MQ=255)
                                    cccTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAAATTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGGGACTAAGTACAGCCAGGTTTTCGTCGGTTTTTGCTTGGAg                                                                                                                                                                                              >  1:145175/1‑229 (MQ=255)
                                      cTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGGCCAG‑ACGCCCCTGGGTTTCCTGAAAATTGAATAAGTCCAGCCAGGTTTGCGCCGGATGTTGGTTGGAGcc                                                                                                                                                                                            >  2:651502/1‑229 (MQ=255)
                                        ttttCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGCGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGACGTTGGTTGGAGCCa                                                                                                                                                                                           >  1:39015/1‑228 (MQ=255)
                                        ttttCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCAt                                                                                                                                                                                          <  1:75479/229‑1 (MQ=255)
                                           tCACCGGGCGTATCCAGGAGGTATTGCGGCTTTGCCCGCGCGGCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTGCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCAc                                                                                                                                                                                       <  2:618033/229‑1 (MQ=255)
                                             aTCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGGGACTAAGTCCGGCAAGGTTTGCGTCGGATGATGGTTGGAGCCATcccc                                                                                                                                                                                      >  1:738466/1‑228 (MQ=255)
                                               cGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGGGACCGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCACGGTTTGCGGCGGAGGTTGGTTGGGGCCAACACCGGc                                                                                                                                                                                   >  2:197962/1‑229 (MQ=255)
                                                gggacatgcccgAATTTTTTGCGGCATTTCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGGTGCCGGCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTCCAGGTCCCCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCTCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCa                                                                                                                                                                                  <  2:317076/218‑1 (MQ=255)
                                                gagaATATCCAGAAAGGTTTGCGGCCTTGCCAGCGCGGCCGGCGGGGGGAAGGAAAAACGGTTGCCGGCGAACTTCGCTAACGCCTGGGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCa                                                                                                                                                                                  <  2:443701/229‑1 (MQ=255)
                                                  gCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGGCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCAtt                                                                                                                                                                                <  1:506563/229‑1 (MQ=255)
                                                            aaTGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGGGGAGATTGTCCAA‑ACGCCCCTGGGTTTCCTGATTAGTGAATAAGTCCAGCAAGGTTTGCGTCGGGTGTTGGTTGGAGCCATACCCggaaatgtaccgggcc                                                                                                                                                                      >  2:175494/1‑215 (MQ=255)
                                                              ttgtttGCGGGAATGCCAGCCCGTCCGGCGCGATGGAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTa                                                                                                                                                                     <  1:389932/224‑1 (MQ=255)
                                                                tATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGACAAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGACATCACCGGCCTTCCGTACCggacca                                                                                                                                                                   >  1:80749/1‑224 (MQ=255)
                                                                 aTTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCCGGTCACCAACCATTGCGACGTGGAGATTGTCCCG‑ACGCCCCCGGGATTCCTGAATCGGGAATAAGTCCCGCAAGGCTTGCGTCTGCTGGTGGTTGGGGCCATCACCGGCATTCAGTACCGGTACaat                                                                                                                                                                 >  2:379145/1‑227 (MQ=255)
                                                                 aTTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGGTTCCTGAATAGGGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCCGGACCGGTACAt                                                                                                                                                                  >  1:15183/1‑228 (MQ=255)
                                                                 aTTGCGGCATTGCCAGCGCGTCCGGCGCGAAGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGACTTCAGGACCGGTACAt                                                                                                                                                                  >  1:266728/1‑228 (MQ=255)
                                                                 aTTGCGGCATTGCCAGCGCGTCCGGCGAGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCCGGGTTTCCGGAATAGTGAATACGGCCAGACAGGTTGTCGTCGGAAGTTGGTTGGAGCCATACCCGGCATTCAAGACCgggcaat                                                                                                                                                                  >  2:107670/1‑223 (MQ=255)
                                                                  ttGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCAAT‑‑GTACCGGTActgaatg                                                                                                                                                              <  1:175494/229‑7 (MQ=255)
                                                                          ttGCCATCACGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTAAACGGTGCGGCCATATTTCAGGTCACCAACCATTACGACGAGAAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGGGAATAAGTCCAGCAAGGTATGCGCCGGATGTGGGTTGGAGCCATCCCCGGCATTCAGGTCCGGGCCATTgagcgaaaa                                                                                                                                                        >  2:574066/1‑221 (MQ=255)
                                                                                            cgATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGgcgcg                                                                                                                                       <  1:552038/228‑1 (MQ=255)
                                                                                            cgATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGgcgcg                                                                                                                                       <  1:514298/228‑1 (MQ=255)
                                                                                            cgATGAAGTAAAAACGGTTGCCGTCCAACTTCACTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATTTTTCAGGTCCCCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGgcgcg                                                                                                                                       <  1:608797/228‑1 (MQ=255)
                                                                                                 aaGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCa                                                                                                                                  <  1:78004/228‑1 (MQ=255)
                                                                                                     aaaacCGGTTTCCGGCCAACTTCCCTAACGCCTGGGGCAGGGAGTGAACGGTGGGGCCATTTTTCAGGGCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTc                                                                                                                             <  2:375657/229‑1 (MQ=255)
                                                                                                           ggTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATACTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCTTGCg                                                                                                                        >  1:342870/1‑228 (MQ=255)
                                                                                                            gttgcgcgcGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCACTATTGCAGGTCCCCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGGCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGa                                                                                                                      <  2:534022/221‑1 (MQ=255)
                                                                                                            ggtGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGACCGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCTCCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGa                                                                                                                      <  1:160790/227‑1 (MQ=255)
                                                                                                             ctGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACATGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTTGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGGGGCCAGGCGCGCCGCACCTTCCTGCgg                                                                                                                       >  1:107615/2‑227 (MQ=255)
                                                                                                               tccGTCGTACTTCGCTAACGCTCGAGTGAGGGAGTGAACGTTGCGGCCATATTTCAGGTCACCAGCCATTGCGACGTGGAGATTGTCCAG‑ACGCCCGTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGa                                                                                                                   <  2:266728/228‑1 (MQ=255)
                                                                                                               gCCGTCCAACTTCGATAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGAt                                                                                                                     <  1:177937/227‑1 (MQ=255)
                                                                                                                    cGAACTTCCCTAACGCCTGAGTCAGGGGGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCAt                                                                                                              <  1:159592/229‑1 (MQ=255)
                                                                                                                        cTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCAc                                                                                                           <  1:197962/228‑1 (MQ=255)
                                                                                                                                   ccTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGGACCGGGACATTGCCGGAAAACTCGGGGGCCAGGCGCGCCGCCCCTTCCTGCGGAGGACGCAACACTATGCCGTCGAc                                                                                               >  2:629799/1‑229 (MQ=255)
                                                                                                                                    ctgggccAGGGGGTGACCGTTGGGGCCAGTTTTCGGGTCACCAACATGTGGGACGT‑GTGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCCCTATCGCATCGACGt                                                                                             <  1:709294/223‑1 (MQ=255)
                                                                                                                                               aGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGAt                                                                                   <  1:552208/229‑1 (MQ=255)
                                                                                                                                                 tGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGTATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGCCgaacgggcggaaaa                                                                                 >  2:220773/1‑216 (MQ=255)
                                                                                                                                                   aaCGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTCGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGACAACTCGTTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGCGGCAAGTGCTGATacc                                                                                >  1:77934/1‑226 (MQ=255)
                                                                                                                                                       gTGCGGCCATATTTCATGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTGGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTaa                                                                                           <  1:603754/213‑1 (MQ=255)
                                                                                                                                                       gTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACCTTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGAAGCCTCACTATCGCAGCGCGGTAAGTGCCGATaaaccgtc                                                                           >  2:603753/1‑223 (MQ=255)
                                                                                                                                                       gTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGGGCCTCCCCTTCCTGCGGATGTCGCATCACTATCGCACCGACGAAAATGCTGCTACACGaa                                                                            >  2:268310/1‑228 (MQ=255)
                                                                                                                                                        ttcGGCCCAATTTCAGGTCACCAACCATTGCGACGTGGAGATTGGCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTa                                                                                                                                                                     <  2:467391/136‑1 (MQ=255)
                                                                                                                                                        tGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTa                                                                                                                                                                     >  1:467391/1‑138 (MQ=255)
                                                                                                                                                                  tttCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGCTAACCGAAATGGTATAGGc                                                                 >  1:101075/1‑228 (MQ=255)
                                                                                                                                                                            ccaaccaTTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGGATCGGCCAGCGTTTCg                                                       >  1:381992/1‑228 (MQ=255)
                                                                                                                                                                                    tGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCCCGATCGCATCGACGTAAGTGCTGATAACCGAAATGGCATCGGCCAGCGTTTCGCCCtt                                                  >  1:431871/1‑225 (MQ=255)
                                                                                                                                                                                    tGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCACCGCCGTAAGGGCTGTTAACCCAAATGGTATCGCCCCACGGTTCGCCCTTTTTa                                              >  2:378339/1‑229 (MQ=255)
                                                                                                                                                                                      cGAGGTGGGGGTTGTCCGG‑ACCGCCCTGGGTTTCCTGAAAAGGGAATAAGTCCAGCAAGGTTTGCGTGGGATGTTGGGTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTAcc                                            <  2:342870/229‑1 (MQ=255)
                                                                                                                                                                                      cGACGGTGAGATTGTCCAG‑GCGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCCAGGTTTGCGTCGGATGTTGGTTGGAGCCAGCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTAcc                                            <  1:651503/229‑1 (MQ=255)
                                                                                                                                                                                       ggattgggtgagtttccggtcGCCCCTGGTTCTCCTGAATGGGCAATAAGTCCGGCTAGGGTTGCGTCGTATGTTGGTTGGCGCCATCACCGGCATTCAGTACCGGTACATTGCCTGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTAcc                                            <  2:107615/209‑1 (MQ=255)
                                                                                                                                                                                         cttGGAGATTTTCCGG‑ACGCCCCGGGGTTTCCGGAATTGGCAATAAGGCCAGGAAGGTGTGCGTCGGATGTTGGTTGGGGCCATCGCCGGCATTCAGGTCCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCCCCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAg                                         <  2:15183/227‑1 (MQ=255)
                                                                                                                                                                                          ggtgggaatttccaaggcGCCCCTGGGTTTCCTGAATGGGGAATAAGGCCAGCAAGGTTTGCGGCGGGTGTGGTTTGGAGCCATCCCCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAg                                         <  2:520069/212‑1 (MQ=255)
                                                                                                                                                                                           tGGAGATTGTCATG‑ACGGCCCTGGTTTTCCGGAATAGTTAATAAGGCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCCCCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTg                                       <  1:357143/229‑1 (MQ=255)
                                                                                                                                                                                                      cctgggcgcccgtggTTTTCGGAAAAGTGAATAAGGCCGGAAAGTTTTGGGTCGGATGTTGGTTGGGGCCATCCCCGGCATCCCGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGgc                             <  2:381992/216‑1 (MQ=255)
                                                                                                                                                                                                        g‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCCCCGGCATTCAGTACCGGTACCTTGCCGGAAAACTCGGTGGCCAGGCCCGCCGCACCTTCCTGCGGATGACGCAACACCTTCGCATCGCAGGACGTCAGGTCAACGGAAATGGAATCGGCCACCGTTTCGCCCTTTTTACCCCGAGACGAATTgccacg                          >  2:2861/1‑224 (MQ=255)
                                                                                                                                                                                                         ggctccctTGGGTTTCCTGAATGTGGAACAAGGTCAGCAATGTTTGAGTCGTGTGTTGGGTGGAGCCATCACGGGAATTCAGTGCGCGTACATTGCCGAAAAACTCGGTGGCCAGGCGGGCCGCACCTTCCTGCGGATGACGCATCACTATCGCGTCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTGCGCCCTTGTTACCCAGTGATGTATTGgcgc                           <  2:376229/224‑1 (MQ=255)
                                                                                                                                                                                                            gtcccTGGGTTTCCTGAATAGGATAGAAGTCCAGCAAGGGTTGTCTGGGATTGTTGTTGGAGCCATCACCGGCATTCAGTACCGGGACATGGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTTATGTATTGGCGCTGt                        <  2:513007/226‑1 (MQ=255)
                                                                                                                                                                                                               ccTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGCGCTGATAACCCAAATGGTTTCGGCCAACCTTTCGCCCTTTTTACCCCGTGATGTATTgggcccttggga                    >  2:339350/1‑219 (MQ=255)
                                                                                                                                                                                                                     tttCCGGAAGAGTAAATAAGTCCAGCAAGGCTGGCGTCGGATGTGAGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACACGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGcc              <  1:508247/229‑1 (MQ=255)
                                                                                                                                                                                                                        ccTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTTCCGGGGAAGCcca            >  1:43705/1‑228 (MQ=255)
                                                                                                                                                                                                                        ccTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGAAGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGAATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGAGGTATTGGCGCTGTCGGGAAAACcca            >  2:193836/1‑228 (MQ=255)
                                                                                                                                                                                                                        ccTGAATAGTGAATAAGGCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCGAGCGTTTGGCCCTTTTTACCCAGTGATGCATTGGCGCTGTCGGGGAAGCcca            >  1:688755/1‑228 (MQ=255)
                                                                                                                                                                                                                         cTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCccac           >  2:373412/1‑228 (MQ=255)
                                                                                                                                                                                                                         cTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGAATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGccccc           >  1:25255/1‑226 (MQ=255)
                                                                                                                                                                                                                             aTAGTTAATTAGTCCCGCAAGGTTGGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGc      <  1:378339/229‑1 (MQ=255)
                                                                                                                                                                                                                                 tGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGCCGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAATCCCACCACGCTGGc  >  1:619362/1‑229 (MQ=255)
                                                                                                                                                                                                                                 tGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGTCGCTGTCGGAGAAGCCCCCACCGCTGGc  >  1:55335/1‑229 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                    
CACTTCTTCAATAGAGCTGTGCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAG‑ACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGC  >  NC_000913/4471740‑4472192

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