Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,888,180 T→A 6.1% C451* (TGT→TGA mnmE → tRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0009133,888,1800T→A 6.1% 4.7 ‑0.5 ‑0.366C451* (TGT→TGAmnmEtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase
Reads supporting (aligned to +/- strand):  new base (2/2):  ref base (18/44):  total (20/46)
Fisher's exact test strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.57e-01

ACGCGTTAATTCGTCTCTCGGCAAGGACTGGTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATG  >  NC_000913/3887882‑3888454
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                  
aCGCGTTAATTCGTCTCTCGGCAAGGACTGGTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTAt                                                                                                                                                                                                                                                                                  <  2:246729/301‑1 (MQ=255)
aCGCGTTAATTCGTCTCTCGGCAAGGACTGGTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTAt                                                                                                                                                                                                                                                                                  <  2:378787/301‑1 (MQ=255)
aCGCGTTAATTCGTCTCTCGGCAAGGACTGGTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTAt                                                                                                                                                                                                                                                                                  <  2:136741/301‑1 (MQ=255)
aCGCGTTAATTCGTCTCTCGGCAAGGACTGGTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTAt                                                                                                                                                                                                                                                                                  <  2:861321/301‑1 (MQ=255)
aCGCGTTAATTCGTCTCTCGGCAAGGACTGGGGAAGGCGTGGACGTGCTGCGTAACCATCTAAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTAt                                                                                                                                                                                                                                                                                  <  1:775385/301‑1 (MQ=255)
  gcgTTAATTCGTCTCTCGGCAAGGACTGGTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTg                                                                                                                                                                                                                                                                                >  2:811529/1‑301 (MQ=255)
     ttAATTCGTCTCTCGGCAAGGACTGGTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGgta                                                                                                                                                                                                                                                                             <  2:716144/301‑1 (MQ=255)
     ttAAGTCGTCTTTCGGCAAGGACTGGGGAAGGCGTGGACGTGCTGGGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGgta                                                                                                                                                                                                                                                                             <  1:416787/301‑1 (MQ=255)
     ttAAGTCGTCTCTCGTCAATGACTTGGGAAGGCGTGGACGTGCTGCGTAACCATTCCAAACAGGGCAGTGGCTTTGACACCAACAGGTAAGGCGGCTTCCTGGCGCGGCGGCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGgta                                                                                                                                                                                                                                                                             <  1:775563/301‑1 (MQ=255)
       tattgttgtttcGGGAAGGAGCGGGGAAGGCTCGGACGGGTGGTGTGATCATCTCAAACAGCGCATGCGCTTTGACACTAACATGGAAGGCGGGTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAGGCGCAACTGTTGGGAGCCTGTGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGgtaag                                                                                                                                                                                                                                                                           <  1:164/291‑1 (MQ=255)
       aaTTCGTCTCTCGGCAAGGACTGGTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGgtaag                                                                                                                                                                                                                                                                           >  2:297486/1‑301 (MQ=255)
           cgtcagtcGGGAATGACTGGGGAAGTGGTGGACGTGGGTATGAACGAACTGCAACCGAGCATGCGGGTTGGAACGAACATGGAGGGGGGCGTCCGGGCGCGTCGGCGGCACCTACAGGCGCTGGAACAGGCAGAGGAACAGCTACAAGAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAAc                                                                                                                                                                                                                                                                       <  1:319177/295‑1 (MQ=255)
           cGTCTCTCGGCAAGGACTGGTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAAc                                                                                                                                                                                                                                                                       <  1:832150/301‑1 (MQ=255)
           cGTCTCTCGGCAAGGACTGGTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAAc                                                                                                                                                                                                                                                                       <  2:909701/301‑1 (MQ=255)
            ggctgtcGGCAAGGACTGGTGAAGGCGTGGACGTGCGGCGTAACCATCTCAATCAGAGCATGGGCTTGGACCCCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGTTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAAc                                                                                                                                                                                                                                                                       <  1:260932/295‑1 (MQ=255)
             tctctcGGCAAGGACTGGTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACcg                                                                                                                                                                                                                                                                     <  1:259779/301‑1 (MQ=255)
                    gCAAGGACTGGTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTa                                                                                                                                                                                                                                                               >  2:879890/1‑300 (MQ=255)
                      aaGGACTGGTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGa                                                                                                                                                                                                                                                            <  2:861214/301‑1 (MQ=255)
                         gACTGGTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCGGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGACATCACCGGGGAATTTACTTCAGACGACCTGCGAGGGCGGATTTTCGCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGCAgg                                                                                                                                                                                                                                                         >  1:770429/1‑300 (MQ=255)
                             ggTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGc                                                                                                                                                                                                                                                      >  2:641323/1‑300 (MQ=255)
                              gTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCAt                                                                                                                                                                                                                                                    >  2:199686/1‑301 (MQ=255)
                                     cggggACGGGCTGCGTAACCATCTCAAACAGAGAAGGGGCTTTGACACCAACATGGAAGGCGGCTTCCGGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGAc                                                                                                                                                                                                                                              <  1:6778/297‑1 (MQ=255)
                                       tGGACGTTCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGTGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTg                                                                                                                                                                                                                                            <  1:879861/300‑1 (MQ=255)
                                         gACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTAGAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCa                                                                                                                                                                                                                                         >  1:540104/1‑301 (MQ=255)
                                                 gcGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGgtaagtaa                                                                                                                                                                                                                                                                        <  1:521587/262‑1 (MQ=255)
                                                 gcGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGgtaagtaa                                                                                                                                                                                                                                                                        >  2:521592/1‑262 (MQ=255)
                                                         aTCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTAGTAAGTAACCGCGCTTACGAAGCCGCATTATGACTGTCAGATGCGgcttcgctt                                                                                                                                                                                                                          >  1:453452/1‑300 (MQ=255)
                                                         aTCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCCCTGGTACTGGCAGCGGAACCTCTACAACAGGTCAAAGCGCCACTGTTGGGAGCCTGTGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCCCAGCCGACCTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTCACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttagctt                                                                                                                                                                                                                          >  1:55/1‑300 (MQ=255)
                                                          tctcAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGCCCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTCACTTTCAGATGCGtcttctcttc                                                                                                                                                                                                                         >  1:939630/1‑300 (MQ=255)
                                                          tctcAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGAATTGGTAAGTAAc                                                                                                                                                                                                                                                                       <  1:614104/254‑1 (MQ=255)
                                                          tctcAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGAATTGGTAAGTAAc                                                                                                                                                                                                                                                                       >  2:614110/1‑254 (MQ=255)
                                                             aaaaCAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCGTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTg                                                                                                                                                                                                                     <  1:199684/300‑1 (MQ=255)
                                                             aaaaCAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCGTCCGGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCTAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTg                                                                                                                                                                                                                     <  1:800550/300‑1 (MQ=255)
                                                             aaaaCAGAGCATGGGCTTTGACACCAACAGGGAAGGCGGCTGCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTg                                                                                                                                                                                                                     <  1:641312/300‑1 (MQ=255)
                                                                   gagCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCa                                                                                                                                                                                                               >  2:471469/1‑301 (MQ=255)
                                                                    ggCATGGGGTTGGACACAAACATGGAAGGCGGCTTCCGGGGGGGTCGTCGCCACCTACAGGCGCGGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCa                                                                                                                                                                                                               <  1:378884/299‑1 (MQ=255)
                                                                    agCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTAGAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAc                                                                                                                                                                                                                 <  2:540109/298‑1 (MQ=255)
                                                                    agCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAcc                                                                                                                                                                                                                <  2:770449/299‑1 (MQ=255)
                                                                            ttttGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTa                                                                                                                                                                                                      <  1:811505/300‑1 (MQ=255)
                                                                            cTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTa                                                                                                                                                                                                      <  2:495844/301‑1 (MQ=255)
                                                                               tGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACCGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTc                                                                                                                                                                                                   <  1:471465/301‑1 (MQ=255)
                                                                               tGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTc                                                                                                                                                                                                   <  2:793656/301‑1 (MQ=255)
                                                                                  cacCAACATGGAAGGCGGCTGCCTGGCGCGTCGTCGCCACCGACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTc                                                                                                                                                                                                <  2:55/301‑1 (MQ=255)
                                                                                       aaaTGGAAGGCGGCTTGCTGGCGCGTCGTCGCCACCGACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAAGCACCGGGGAATTTACTGCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAccc                                                                                                                                                                                           <  1:399094/299‑1 (MQ=255)
                                                                                       aCATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAccc                                                                                                                                                                                           <  2:939671/301‑1 (MQ=255)
                                                                                        cATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTgg                                                                                                                                                                                                                                                                               >  2:875292/1‑216 (MQ=255)
                                                                                        cATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTgg                                                                                                                                                                                                                                                                               <  1:875263/216‑1 (MQ=255)
                                                                                                   gcTTCCTGGCGCGTCGTCGCCACCTACAGGCGCGGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                                                                                                                                                                <  1:297483/300‑1 (MQ=255)
                                                                                                           gcacgacctcgccagctcgcggCGCTGGTAAAGGCAGCGGAACATCGTCAACAGGCCAGAGCGCAACCGTTGGGAGCCGGCGCAAGGCAACTGCTGGCGTGAGAGTCTCCTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTAATTCAGACGACCGGCTGGGGCGGATTTTCTCCAGCTCCTGTATTGGTAAGTAACCGCCCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAAt                                                                                                                                                                        <  1:105637/280‑1 (MQ=255)
                                                                                                                        cccTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta                                                                                                                                                          <  2:453456/300‑1 (MQ=255)
                                                                                                                          cTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTtata                                                                                                                                                        >  1:642476/1‑301 (MQ=255)
                                                                                                                                gcgcTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATc                                                                                                                                                     <  2:642487/298‑1 (MQ=255)
                                                                                                                                               gCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGGACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                                                                                                                                                             <  2:91254/243‑1 (MQ=255)
                                                                                                                                                  gaagttCTACAAGAGGGCAAAGCGCAACTGTGGGAGCCTGTGCAGGTGAACTGCTGGCGGAAGAGTTGAGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAatat                                                                                                                                <  1:342934/295‑1 (MQ=255)
                                                                                                                                                     aaTCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGc                                                                                                                             <  2:332885/300‑1 (MQ=255)
                                                                                                                                                                            aCTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgattcgat                                                                                                      >  2:737635/1‑301 (MQ=255)
                                                                                                                                                                                                   aaCTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTTACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGaa                                                                                <  2:237802/300‑1 (MQ=255)
                                                                                                                                                                                                     ctgctgGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTc                                                                                                                                                                                                   <  1:221952/183‑1 (MQ=255)
                                                                                                                                                                                                     ctgctgGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTc                                                                                                                                                                                                   >  2:221954/1‑183 (MQ=255)
                                                                                                                                                                                                              gAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAAc                                                                     <  1:737617/300‑1 (MQ=255)
                                                                                                                                                                                                                             gCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTa                                                                                                                                                                                                      <  2:296111/156‑1 (MQ=255)
                                                                                                                                                                                                                             gCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTa                                                                                                                                                                                                      >  1:296108/1‑156 (MQ=255)
                                                                                                                                                                                                                                                   tCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGAATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAatat                                                                                                                                >  1:11710/1‑204 (MQ=255)
                                                                                                                                                                                                                                                   tCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGAATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAatat                                                                                                                                <  2:11710/204‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ctgctgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCCCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACatat     >  2:629464/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                 tgGGGCGGATTTTCTCCAGCTTCTGTATTTGTAAGTTACCGCGCTTACGAAGCCGCATGCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATGTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGAGTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATAtg  <  1:426909/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                  
ACGCGTTAATTCGTCTCTCGGCAAGGACTGGTGAAGGCGTGGACGTGCTGCGTAACCATCTCAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATG  >  NC_000913/3887882‑3888454

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.