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

CACTTCTTCAATAGAGCTGTGCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGC  >  NC_000913/4471740‑4472192
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                    
cACTTCTTCAATAGAGCTGTGCAGACTCCATGCGCTCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGGGCGGCCAAATTTCAGGCCACCAACCATTGGGACGTGGGAGTGTCCAGAACCCCCCGGGGTTCCCGGAAAAGGGaaa                                                                                                                                                                                                                                  >  2:709293/1‑228 (MQ=255)
         atatgggttgtgCAACCCCCATGCCATCCCTTTTTCATCGAGCATAGCCAGAATGGATTGCGGCTTTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTGCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCACCTGGGTTTCCTGAATAGTGAATAAGTCCAGc                                                                                                                                                                                                                         <  2:412575/221‑1 (MQ=255)
           tAGAGCTGTGCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGACGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCCCCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATTGGGAAAAAGTCCAGCaa                                                                                                                                                                                                                       >  2:506563/1‑229 (MQ=255)
           tAGAGCTGTGCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTACCGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGGGATTGTCCCGACGCCCCCGGGTTTCCGGAATAGTGAATAAATCCAGCaa                                                                                                                                                                                                                       >  2:272830/1‑229 (MQ=255)
           tAGAGCTGTGCAGACTCCATGCGATCCCTTTTGCATCGAGCATATCCAGAATGGATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGc                                                                                                                                                                                                                         <  1:272830/227‑1 (MQ=255)
                cTGTGCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATATGGAATAAGTCCAGCAAGGttt                                                                                                                                                                                                                  >  1:375657/1‑229 (MQ=255)
                cTGTGCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGGAAATAAGCCCAGCAAGGttt                                                                                                                                                                                                                  >  1:384206/1‑229 (MQ=255)
                cTGTGCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGGGCGGCCATATTTCAGGTCACCACCCATTGCGCCGTGGGGATTGTCCAGACGCCCCTGGGGTTCCTGAATAGGGAATAAGTCCAGAAAgtgtt                                                                                                                                                                                                                  >  2:608796/1‑225 (MQ=255)
                cTGTGCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGAGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACTATTGCGACGTGGGGCTTGTCCGAATGCCCCTGGGGTTCCTGAATATGGAATAAGTCCAGCAAGGatt                                                                                                                                                                                                                  >  1:432505/1‑229 (MQ=255)
                   tgCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCg                                                                                                                                                                                                               >  2:78004/1‑229 (MQ=255)
                   tgCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTCCCTGAATAGGGATTAAGTCAAGCAAGGTTTGct                                                                                                                                                                                                               >  1:458015/1‑228 (MQ=255)
                   tgCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGCCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGAAAGGTTTGCg                                                                                                                                                                                                               >  1:520069/1‑229 (MQ=255)
                    ccAGACTCCATGTGATCCCTTGTTCATCGAGCACATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGt                                                                                                                                                                                                              <  1:333770/228‑1 (MQ=255)
                        aCTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGGGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAAAAGTGAATAAGTCCAGAAGAGTTTGGGCCggg                                                                                                                                                                                                          >  2:177937/1‑228 (MQ=255)
                         cTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGGCCGGCGCGGTGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGGCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCgg                                                                                                                                                                                                           <  1:550402/227‑1 (MQ=255)
                         cTCCAAGCGATCCCTTTTACATCGGGCCTATCCAGAATGTATTGCGGCATTTCCAGCGGATCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGGGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCgg                                                                                                                                                                                                           <  1:2614/227‑1 (MQ=255)
                             aTGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATgtt                                                                                                                                                                                                      >  2:88268/1‑228 (MQ=255)
                                  aTCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCCTATTTCAGTTCACCACCCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCAAGCACGGTTTGCGCGGGATGtgggttgg                                                                                                                                                                                                >  2:159592/1‑229 (MQ=255)
                                   tCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATgttggtt                                                                                                                                                                                                  <  1:88268/226‑1 (MQ=255)
                                    cccTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAAATTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGGGACTAAGTACAGCCAGGTTTTCGTCGGTTTTTGCTTGGAg                                                                                                                                                                                              >  1:145175/1‑229 (MQ=255)
                                      cTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGGCCAGACGCCCCTGGGTTTCCTGAAAATTGAATAAGTCCAGCCAGGTTTGCGCCGGATGTTGGTTGGAGcc                                                                                                                                                                                            >  2:651502/1‑229 (MQ=255)
                                        ttttCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGCGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGACGTTGGTTGGAGCCa                                                                                                                                                                                           >  1:39015/1‑228 (MQ=255)
                                        ttttCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCAt                                                                                                                                                                                          <  1:75479/229‑1 (MQ=255)
                                           tCACCGGGCGTATCCAGGAGGTATTGCGGCTTTGCCCGCGCGGCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTGCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCAc                                                                                                                                                                                       <  2:618033/229‑1 (MQ=255)
                                             aTCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGGGACTAAGTCCGGCAAGGTTTGCGTCGGATGATGGTTGGAGCCATcccc                                                                                                                                                                                      >  1:738466/1‑228 (MQ=255)
                                               cGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGGGACCGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCACGGTTTGCGGCGGAGGTTGGTTGGGGCCAACACCGGc                                                                                                                                                                                   >  2:197962/1‑229 (MQ=255)
                                                gggacatgcccgAATTTTTTGCGGCATTTCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGGTGCCGGCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTCCAGGTCCCCAACCATTGCGACGTGGAGATTGTCCAGACGCCTCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCa                                                                                                                                                                                  <  2:317076/218‑1 (MQ=255)
                                                gagaATATCCAGAAAGGTTTGCGGCCTTGCCAGCGCGGCCGGCGGGGGGAAGGAAAAACGGTTGCCGGCGAACTTCGCTAACGCCTGGGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCa                                                                                                                                                                                  <  2:443701/229‑1 (MQ=255)
                                                  gCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGGCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCAtt                                                                                                                                                                                <  1:506563/229‑1 (MQ=255)
                                                            aaTGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGGGGAGATTGTCCAAACGCCCCTGGGTTTCCTGATTAGTGAATAAGTCCAGCAAGGTTTGCGTCGGGTGTTGGTTGGAGCCATACCCggaaatgtaccgggcc                                                                                                                                                                      >  2:175494/1‑215 (MQ=255)
                                                              ttgtttGCGGGAATGCCAGCCCGTCCGGCGCGATGGAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTa                                                                                                                                                                     <  1:389932/224‑1 (MQ=255)
                                                                tATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGACAAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGACATCACCGGCCTTCCGTACCggacca                                                                                                                                                                   >  1:80749/1‑224 (MQ=255)
                                                                 aTTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCCGGTCACCAACCATTGCGACGTGGAGATTGTCCCGACGCCCCCGGGATTCCTGAATCGGGAATAAGTCCCGCAAGGCTTGCGTCTGCTGGTGGTTGGGGCCATCACCGGCATTCAGTACCGGTACaat                                                                                                                                                                 >  2:379145/1‑227 (MQ=255)
                                                                 aTTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGGTTCCTGAATAGGGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCCGGACCGGTACAt                                                                                                                                                                  >  1:15183/1‑228 (MQ=255)
                                                                 aTTGCGGCATTGCCAGCGCGTCCGGCGCGAAGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGACTTCAGGACCGGTACAt                                                                                                                                                                  >  1:266728/1‑228 (MQ=255)
                                                                 aTTGCGGCATTGCCAGCGCGTCCGGCGAGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCCGGGTTTCCGGAATAGTGAATACGGCCAGACAGGTTGTCGTCGGAAGTTGGTTGGAGCCATACCCGGCATTCAAGACCgggcaat                                                                                                                                                                  >  2:107670/1‑223 (MQ=255)
                                                                  ttGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCAAT‑‑GTACCGGTActgaatg                                                                                                                                                              <  1:175494/229‑7 (MQ=255)
                                                                          ttGCCATCACGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTAAACGGTGCGGCCATATTTCAGGTCACCAACCATTACGACGAGAAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGGGAATAAGTCCAGCAAGGTATGCGCCGGATGTGGGTTGGAGCCATCCCCGGCATTCAGGTCCGGGCCATTgagcgaaaa                                                                                                                                                        >  2:574066/1‑221 (MQ=255)
                                                                                            cgATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGgcgcg                                                                                                                                       <  1:552038/228‑1 (MQ=255)
                                                                                            cgATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGgcgcg                                                                                                                                       <  1:514298/228‑1 (MQ=255)
                                                                                            cgATGAAGTAAAAACGGTTGCCGTCCAACTTCACTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATTTTTCAGGTCCCCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGgcgcg                                                                                                                                       <  1:608797/228‑1 (MQ=255)
                                                                                                 aaGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCa                                                                                                                                  <  1:78004/228‑1 (MQ=255)
                                                                                                     aaaacCGGTTTCCGGCCAACTTCCCTAACGCCTGGGGCAGGGAGTGAACGGTGGGGCCATTTTTCAGGGCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTc                                                                                                                             <  2:375657/229‑1 (MQ=255)
                                                                                                           ggTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATACTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCTTGCg                                                                                                                        >  1:342870/1‑228 (MQ=255)
                                                                                                            gttgcgcgcGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCACTATTGCAGGTCCCCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGGCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGa                                                                                                                      <  2:534022/221‑1 (MQ=255)
                                                                                                            ggtGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGACCGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCTCCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGa                                                                                                                      <  1:160790/227‑1 (MQ=255)
                                                                                                             ctGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACATGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTTGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGGGGCCAGGCGCGCCGCACCTTCCTGCgg                                                                                                                       >  1:107615/2‑227 (MQ=255)
                                                                                                               tccGTCGTACTTCGCTAACGCTCGAGTGAGGGAGTGAACGTTGCGGCCATATTTCAGGTCACCAGCCATTGCGACGTGGAGATTGTCCAGACGCCCGTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGa                                                                                                                   <  2:266728/228‑1 (MQ=255)
                                                                                                               gCCGTCCAACTTCGATAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGAt                                                                                                                     <  1:177937/227‑1 (MQ=255)
                                                                                                                    cGAACTTCCCTAACGCCTGAGTCAGGGGGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCAt                                                                                                              <  1:159592/229‑1 (MQ=255)
                                                                                                                        cTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCAc                                                                                                           <  1:197962/228‑1 (MQ=255)
                                                                                                                                   ccTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGGACCGGGACATTGCCGGAAAACTCGGGGGCCAGGCGCGCCGCCCCTTCCTGCGGAGGACGCAACACTATGCCGTCGAc                                                                                               >  2:629799/1‑229 (MQ=255)
                                                                                                                                    ctgggccAGGGGGTGACCGTTGGGGCCAGTTTTCGGGTCACCAACATGTGGGACGTGTGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCCCTATCGCATCGACGt                                                                                             <  1:709294/223‑1 (MQ=255)
                                                                                                                                               aGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGAt                                                                                   <  1:552208/229‑1 (MQ=255)
                                                                                                                                                 tGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGTATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGCCgaacgggcggaaaa                                                                                 >  2:220773/1‑216 (MQ=255)
                                                                                                                                                   aaCGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTCGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGACAACTCGTTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGCGGCAAGTGCTGATacc                                                                                >  1:77934/1‑226 (MQ=255)
                                                                                                                                                       gTGCGGCCATATTTCATGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTGGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTaa                                                                                           <  1:603754/213‑1 (MQ=255)
                                                                                                                                                       gTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACCTTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGAAGCCTCACTATCGCAGCGCGGTAAGTGCCGATaaaccgtc                                                                           >  2:603753/1‑223 (MQ=255)
                                                                                                                                                       gTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGGGCCTCCCCTTCCTGCGGATGTCGCATCACTATCGCACCGACGAAAATGCTGCTACACGaa                                                                            >  2:268310/1‑228 (MQ=255)
                                                                                                                                                        ttcGGCCCAATTTCAGGTCACCAACCATTGCGACGTGGAGATTGGCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTa                                                                                                                                                                     <  2:467391/136‑1 (MQ=255)
                                                                                                                                                        tGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTa                                                                                                                                                                     >  1:467391/1‑138 (MQ=255)
                                                                                                                                                                  tttCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGCTAACCGAAATGGTATAGGc                                                                 >  1:101075/1‑228 (MQ=255)
                                                                                                                                                                            ccaaccaTTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGGATCGGCCAGCGTTTCg                                                       >  1:381992/1‑228 (MQ=255)
                                                                                                                                                                                    tGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCCCGATCGCATCGACGTAAGTGCTGATAACCGAAATGGCATCGGCCAGCGTTTCGCCCtt                                                  >  1:431871/1‑225 (MQ=255)
                                                                                                                                                                                    tGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCACCGCCGTAAGGGCTGTTAACCCAAATGGTATCGCCCCACGGTTCGCCCTTTTTa                                              >  2:378339/1‑229 (MQ=255)
                                                                                                                                                                                      cGAGGTGGGGGTTGTCCGGACCGCCCTGGGTTTCCTGAAAAGGGAATAAGTCCAGCAAGGTTTGCGTGGGATGTTGGGTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTAcc                                            <  2:342870/229‑1 (MQ=255)
                                                                                                                                                                                      cGACGGTGAGATTGTCCAGGCGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCCAGGTTTGCGTCGGATGTTGGTTGGAGCCAGCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTAcc                                            <  1:651503/229‑1 (MQ=255)
                                                                                                                                                                                       ggattgggtgagtttccggtcGCCCCTGGTTCTCCTGAATGGGCAATAAGTCCGGCTAGGGTTGCGTCGTATGTTGGTTGGCGCCATCACCGGCATTCAGTACCGGTACATTGCCTGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTAcc                                            <  2:107615/209‑1 (MQ=255)
                                                                                                                                                                                         cttGGAGATTTTCCGGACGCCCCGGGGTTTCCGGAATTGGCAATAAGGCCAGGAAGGTGTGCGTCGGATGTTGGTTGGGGCCATCGCCGGCATTCAGGTCCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCCCCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAg                                         <  2:15183/227‑1 (MQ=255)
                                                                                                                                                                                          ggtgggaatttccaaggcGCCCCTGGGTTTCCTGAATGGGGAATAAGGCCAGCAAGGTTTGCGGCGGGTGTGGTTTGGAGCCATCCCCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAg                                         <  2:520069/212‑1 (MQ=255)
                                                                                                                                                                                           tGGAGATTGTCATGACGGCCCTGGTTTTCCGGAATAGTTAATAAGGCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCCCCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTg                                       <  1:357143/229‑1 (MQ=255)
                                                                                                                                                                                                      cctgggcgcccgtggTTTTCGGAAAAGTGAATAAGGCCGGAAAGTTTTGGGTCGGATGTTGGTTGGGGCCATCCCCGGCATCCCGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGgc                             <  2:381992/216‑1 (MQ=255)
                                                                                                                                                                                                        gACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCCCCGGCATTCAGTACCGGTACCTTGCCGGAAAACTCGGTGGCCAGGCCCGCCGCACCTTCCTGCGGATGACGCAACACCTTCGCATCGCAGGACGTCAGGTCAACGGAAATGGAATCGGCCACCGTTTCGCCCTTTTTACCCCGAGACGAATTgccacg                          >  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)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                    
CACTTCTTCAATAGAGCTGTGCAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGC  >  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: