Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,728,949 T→C 15.0% E141G (GAG→GGG)  xylB ← xylulokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,728,9490TC15.0% 201.6 / 3.7 84E141G (GAG→GGG) xylBxylulokinase
Reads supporting (aligned to +/- strand):  ref base T (48/17);  new base C (3/13);  total (52/32)
Fisher's exact test for biased strand distribution p-value = 8.00e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCAATGCTTTAACGTCCT  >  NC_000913/3728732‑3729174
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                   
aaaaaaaGCCCCAGGAAGTGCGCTGCCTTGTGAATAGGGCGGACTCTGGTCACGAGGAAGGTCGCAAGCCCGAAGTCGGACGTCACTCCAGTCACGCTTTGCAACATCCAGCCAAATGGTGCCAGCTGCGTCAGACATAGCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATg                                                                                                                                                                                                                          <  1:504821/224‑1 (MQ=255)
aaaCAAAGCACCACTCATTCCCCTGCATACTTATAATGCGGGAAACTGGTCACGAGATGAGCCGCAACCCTGCAGCTTCCGGTCACTCCAATCACTCTTTGCGACATCCAGCCACCTGGTGCCAGTTGCGTCAGACATATCGCTGGCATACTCCCCCGTCATACGCAGACGCAAGTAGTCTTTCGGTAATAATCCTTTGTCGATTTGACGGAATATCTCCGGCTCATg                                                                                                                                                                                                                          <  1:204755/227‑1 (MQ=255)
              tAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCCCTCTAGTCACGCTTTGCGACATCAAGCCACATGGTGCCAGCTGCGTCCGACATCTCGCCGGCAAACTCCCCCGTCATACGCAGACGCAAGAAATCTTTCGGTAATAATACTTTGCCGCTTTGCCGACACACGCCCGGCTCAGGCCACAGAACCCATa                                                                                                                                                                                                            >  1:297605/1‑228 (MQ=255)
                        gCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCACCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCACACGCCAGCAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCCTGCCCCTAAACCCATAGCAATTTAga                                                                                                                                                                                                  >  1:137259/1‑227 (MQ=255)
                        gCCTTCGTATAATGCGGGAATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGTCGTAACTCCAGTCACGCATTGCGACATCCAGACACATGGTGCCAGCTGCGTCAGACATACAGCTGGCAAAATCCCCATTCATACGCAGACGCAAGTAATCCTTCGGTACTACTACTTTGTGTGTTTGACGTAATAGCTCCGACTCAGTCCGCTTGACCCTTAGCAACTTAg                                                                                                                                                                                                   >  2:319176/1‑227 (MQ=255)
                        gCCTTAGTATAATGCGGGCATCTGGTCACGAAATAAGTCGCAAGCCTGCAGCATGCCGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAgg                                                                                                                                                                                                  <  2:75147/228‑1 (MQ=255)
                             ttggataACGCGGGCATCTGGTACCGAGATAAGTCGCAAGCCTGCAGCATGACGTCCCCCCAGTCACGCTTTGCGACATCCAGCCACAGGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCGTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCa                                                                                                                                                                                              <  2:592557/224‑1 (MQ=255)
                                     tGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGAAGACGCAAGGAATCTTTCGGTAATAATACTTTGTCGATTGGACGGAATATCCCCGGTCACTGCCCCTGAACCCTTACCAATTACGGCGCAGAAAATcc                                                                                                                                                                                      >  1:75147/1‑228 (MQ=255)
                                         ggCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCTACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATCCGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAACATCTCCGGCACATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCggg                                                                                                                                                                                   >  2:90089/1‑227 (MQ=255)
                                         ggCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAGTATCTCCGGCTCATGCCGCTGAACCCATCACAATTTAGGCGCAGTAAATCCggg                                                                                                                                                                                   >  1:613599/1‑227 (MQ=255)
                                         ggCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGCAAATCCGGGc                                                                                                                                                                                  >  2:513278/1‑228 (MQ=255)
                                         ggCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCACCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATCCGCAGACGCAAGGAATCTTTCGGGAATAATACTTTGTCGATTTGACGGAAAAACACCGACTCATGCCCCTGAACCCTCAGCACTTTAGGGCCAGTAAACCCcgg                                                                                                                                                                                   >  1:62226/1‑227 (MQ=255)
                                          gCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGCTTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGc                                                                                                                                                                                  >  2:304624/1‑227 (MQ=255)
                                           cATCTGGTAACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGCAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCGGCTGAACCCATAGCAATTTAGGCGCAGTACATCCGGGca                                                                                                                                                                                 >  2:203898/1‑227 (MQ=255)
                                                   cccGAGATAAGCCGCAAGCCTGCAGCTGTACGTCCTCCCAGTCACCTTTTGCGACTACCAGCCCCATGGTGCCAGCTGCGTTCGACATAGCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGTAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGt                                                                                                                                                                        <  1:501604/226‑1 (MQ=255)
                                                      gagaTAAGTCGCAAGCCTGCATCATGACGTCACTCCAGTCACGCTTTGCGCCATCCATCCACATGGTGCCAGCTGCGTCAGACATATCTCTGGCATCCTCCCCCGTCATACCCAGCCGCAAGTAATCTTTCGGTAATAATACTTGCTAGCTTTGACGGAATATCCACGCATCATGCCGCTGAACCCAAAGCAATTATCGCGCAAAAAAACCggcagtacacaggatcc                                                                                                                                                                     >  1:337387/1‑213 (MQ=255)
                                                       agaTAAGTCGCAAGCCTGAAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGc                                                                                                                                                                     <  2:652792/227‑1 (MQ=255)
                                                                       tGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTACTAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGGGCAGTAAATCCTGGCACCATCCGGACGCCGGTAAACACCCGCGaa                                                                                                                                                    >  1:115180/1‑227 (MQ=255)
                                                                       tGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCCCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGa                                                                                                                                                     >  2:124658/1‑227 (MQ=255)
                                                                        gcagcaTGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGAt                                                                                                                                                    >  2:240657/1‑227 (MQ=255)
                                                                         cagcaTGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGGAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCcgggctt                                                                                                                                                   >  2:231076/1‑222 (MQ=255)
                                                                                      cTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCTTACCAATTTAGCCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCt                                                                                                                                     >  1:8710/1‑228 (MQ=255)
                                                                                       tCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGCCGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAgg                                                                                                                                                                                                  >  1:173076/1‑166 (MQ=255)
                                                                                       tCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAgg                                                                                                                                                                                                  <  2:173076/166‑1 (MQ=255)
                                                                                             ccccgtttggcccctCCAGCCCCATGTGGCCAAGTGCGTCGAACTTATCGCTGGAAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCa                                                                                                                              <  1:474867/214‑1 (MQ=255)
                                                                                                 cTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTCGGGAACTCGCGCTTCCAGCAAAg                                                                                                                           >  1:462190/1‑227 (MQ=255)
                                                                                                   ttGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGGGCCGTAAATCCGGGCATCCTCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCaaaagga                                                                                                                        >  1:45911/1‑224 (MQ=255)
                                                                                                    tGCCACACCCCGCCACATGGTGCCAGGTGCGTCAGACATTTCGCTGGAAAACTCCCCCGTCATTCGCAGACGCAAGTAATCTTCCGGTAATAATACTTTGTCGATTTGACGGAATATCCCCGGCTCATGCCGCTGATCCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCa                                                                                                                       <  1:240657/228‑1 (MQ=255)
                                                                                                      cGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCATTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGGGCAcc                                                                                                                     >  1:549594/1‑227 (MQ=255)
                                                                                                              cgCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATGCGCACACTCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCAAAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGca                                                                                                             <  2:401020/227‑1 (MQ=255)
                                                                                                              aGCCACCTGGTTCCAGCTGCGTCAGACCTATTGCTGGCAAACTCCCCCGTCATCAGCAGACGAAAGTAATCTTTCGGTAATAATACTTTGGCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGca                                                                                                             <  2:337387/228‑1 (MQ=255)
                                                                                                                ccACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGCCGATTTGACGGAATATCTCCGGCTCATGTCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATACTCCGGTGTCCGGTAATAACCCGGGATTGCGGAACTCGCACTTCCAGCAAAAGGCACCCTTgtgc                                                                                                              >  1:404453/1‑225 (MQ=255)
                                                                                                                        tGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCTTCATCAGATTGCCGGTAATCACCCGCGATTGCAGAACTCGCACTTCCAGCAAAAGGCACTCTTGCCCCCAGCCCCcgc                                                                                                   >  1:679765/1‑227 (MQ=255)
                                                                                                                          ccAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCTATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAATGCACTCTTGCTCCCAGCGCCcgtcg                                                                                                 >  2:19374/1‑228 (MQ=255)
                                                                                                                             gCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCTGGCATCATCAGGTTGCCGGAAATAACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGGGGCCTCTGCCGCCCAGCGCCCGTCGTTc                                                                                              >  1:13841/1‑228 (MQ=255)
                                                                                                                                 cGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGGCGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCcaca                                                                                          <  2:137259/228‑1 (MQ=255)
                                                                                                                                  gTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTCCGCACAGCGCCCGTCGTTCCACaa                                                                                         >  2:111247/1‑228 (MQ=255)
                                                                                                                                        cATACCGCTGGCAAATCACCCCGTCATACGGAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGc                                                                                  <  2:9361/228‑1 (MQ=255)
                                                                                                                                         ttttctctgcaaaaCTCCCCCGTCATCCGCAGACGCAATTAAGCTTTCGGTAATAATACTTTGGCGATTTGACGGAATATCCCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCa                                                                                 <  1:111247/217‑1 (MQ=255)
                                                                                                                                         atatCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATCCTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCa                                                                                 <  2:449129/228‑1 (MQ=255)
                                                                                                                                          tatggagtgCAAACCCCCCCGTTAAACGCCAACGCAAGTAATCTTTCGGTAATAATACTTTGCCGATTTGGCGGAATATCGCCCGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAg                                                                                <  1:513278/220‑1 (MQ=255)
                                                                                                                                           atCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACACCGCCCGTCGTTCCCCAAAATGGCAgg                                                                               >  2:420884/1‑228 (MQ=255)
                                                                                                                                           atCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCCTCCCGCAAAGTGCACTCTTGCGCACAGCCCCCGTCGTTCCCCAAAATGCCAgg                                                                               >  2:327457/1‑228 (MQ=255)
                                                                                                                                                 ggCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCg                                                                            >  2:173574/1‑225 (MQ=255)
                                                                                                                                                     taCTCCCCCGTCATACGCAGACGCAAGGAAGCGTCCGGTAATAATACGTTGGCGATGTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGGAAATCCGGGCATCATCAGGTTGCCGGTAATCCCCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACAcc                                                                     <  2:418317/227‑1 (MQ=255)
                                                                                                                                                               tCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCACCAAAGTGACCTCTTCCGCCCAGCGCCCGCCGTTCCACAAAATGGCAGGCCGAAACAcctctcccttaa                                                           >  1:451470/1‑218 (MQ=255)
                                                                                                                                                                cATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCa                                                                                 <  2:451470/205‑1 (MQ=255)
                                                                                                                                                                      cAGACGGAAGTCATCTTTCCGTTAAATACCTTTGTCGATTTGACGGAATATCCCCGGCTCATGCCCCTGGACCCATAGCAATTTAGGCGCAGTAAATCCTGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAg                                                    <  1:203898/228‑1 (MQ=255)
                                                                                                                                                                      aaagcGCAAGTAATCTGTCGGTAATAATATTTTGTCGATTTGACGGAATATCCCCGGCTCATGCCGCTGAACCCAGAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCCCCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAg                                                    <  1:173574/224‑1 (MQ=255)
                                                                                                                                                                        gaagCAAGTTATTTTTTGGTAATAATTCTTTTTCAATTTTACGGAATTTCCCCGGCTCATTCCGCTGGACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCa                                                  <  1:420884/225‑1 (MQ=255)
                                                                                                                                                                         aCGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCACTTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGCCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGAATCCAGcca                                                 >  2:92515/1‑226 (MQ=255)
                                                                                                                                                                         aCGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGCTTGCGGAACTCGCGCTTCCAGCAAAGTGACCTCTTGCGCCCAGCGCCCGTCGTTCACAACAAAGGACGGAGGTACCACCCGTTGCtcgccatcccccaa                                                 >  1:5050/1‑215 (MQ=255)
                                                                                                                                                                         aCGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCCCAAAATGGCAGGGCGTACACCCCGTTGCTGGCCCTCCACCaa                                                 >  1:123368/1‑228 (MQ=255)
                                                                                                                                                                         aCGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTCCGCACAGCGCCCGACGTTCCACAAAATGGACGGGCGTAACACCGCTTTCTGacataccaccaa                                                 >  2:330447/1‑217 (MQ=255)
                                                                                                                                                                         aCGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGAGCATTAAATCCGGGCATCATCAGGTTGCCGGTAATCAACCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCCCCTGGCCCGCCGTTCCCCAAATTGACAGGGCGAAACcccccctaccaggcaaccaccca                                                 >  1:293575/1‑209 (MQ=255)
                                                                                                                                                                                   tCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCCCAAAAGGGAAGGCCGTAACACCCCTTGCTCACCACCCacaaacgagtctccg                                       >  1:544590/1‑214 (MQ=255)
                                                                                                                                                                                       tCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCCCAAAATGGCAGGGCGTAACACCCGTTGCTGCGCATCCCGCAAGGTTGCTCCGTGCa                                   >  2:342521/1‑228 (MQ=255)
                                                                                                                                                                                       tCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTcctggca                                   >  2:514009/1‑223 (MQ=255)
                                                                                                                                                                                       tCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAAAACCCGTTGCTTAGCATCCAGCAAGGTTGCTCCGTGCa                                   >  2:621837/1‑228 (MQ=255)
                                                                                                                                                                                          cttattaaatcTTTGTCGATTTTACGGAATTACGCCCGGTCAAGCCGCTGGACCCATAGCCATTTAGGCGCAGAAAATCCGGGCATCGTCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCCTGCCCACAGTGCCCGTCGTTCCACAAAATGGCAGGGCGTAACCCCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCAt                                  <  1:305983/217‑1 (MQ=255)
                                                                                                                                                                                           taataaTACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTACCACCCGTTTCTGAGCATCCAGCAAGGTTGCTCCGCGCa                                   >  2:305983/1‑224 (MQ=255)
                                                                                                                                                                                             acaattCTTTGGCGATTTGACGGAAATCCCCCGGCTAAGGCCGCGGAACCCATAGCAATTTAGGCGCAGTAAAACCCGGCATCATCAGGTTGCCGGTAATCCCCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACCCCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGc                             <  1:514009/226‑1 (MQ=255)
                                                                                                                                                                                              ataaccatttGTCGATTTTGCGGAATTTTCCCGGGCAACGCCGCTGAACCCCAAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGc                             <  1:265356/221‑1 (MQ=255)
                                                                                                                                                                                                  aCTTTTTCGTTTTGGCGGAATATCCCCGGCTCATGCCGTTGAACCCAGAGCAATTTGGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAACCACCCGCGATTGCGGAACCCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACCCCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCa                        <  1:327457/228‑1 (MQ=255)
                                                                                                                                                                                                   ttttttTCGATTTGCAGGAATTTCCCCGGCACATGCCGCTGAACCCAAAGCAATTGGGGCGCAGTAAATCCGGGCATCACCAGGTTGCCCGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACGCCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCa                        <  1:342521/226‑1 (MQ=255)
                                                                                                                                                                                                      ttGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACTCATAGCAATTTAGGCTCAGTAAATCCAGGACACATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTACCAGCAAAATGCACTCTTGCGCCAACTGCCCGACGCTCCCGAAAATAGCCAGGCGACAGACCCGTTCCGTGGACTCCAGCAAAGATGCTCCAGCCACATGGCCGGCacaa                     >  1:671067/1‑225 (MQ=255)
                                                                                                                                                                                                      ttGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGCTTGCGGAACTCGCGCTTCCCGCAACGAGCACCCTTGCGCCCCGCGCCCGTCGTTCCGCCAAATGGCAGGGCGTACAACCCGTGGCTCAGCATCCCGCAACGTTGCTCCTGGGCTCTGGCCGCCatac                     >  1:132321/1‑225 (MQ=255)
                                                                                                                                                                                                      ttGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCCGCAAAGTGCACTCTTGCGCACAGCGCCcgtc                                                                                                  >  1:268590/1‑152 (MQ=255)
                                                                                                                                                                                                      ttGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCCAAGTGCACTCTTCCGCACAACGCCCGTCGTTCACCAAAATGCCAGGGCTCACCACCCGTGCCCGAGCATCCAGACAGGCTGCTCCGCGAAACTGGCCGCCaaaa                     >  1:323558/1‑226 (MQ=255)
                                                                                                                                                                                                      ttGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGCTGTGCCGTGCATCTGCCCGGCAATa                     >  2:560574/1‑228 (MQ=255)
                                                                                                                                                                                                      ttGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCAAGCAAGGTTGCTCCGTGCACCTGACCGGCCATa                     >  2:352608/1‑228 (MQ=255)
                                                                                                                                                                                                      ttGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGAAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGATTCCAGCAAAGTTCACTCTTGCGCCCAGCGCCCGCCGTTCCACAAAATGGACGGGCGAACCACCCGTTGCTGACCACCAGGCAAGGATCCTCCGAGCATGTGCCCAGCaaac                     >  1:199800/1‑226 (MQ=255)
                                                                                                                                                                                                      gtgccGATTTGACGGAATACGTCTGGCTCAAGACGCTGAAACCTTAACAATTTAGGGGCAGTAAATCAGGGCATCATCAGGTTGCCGGTAATCAACCGCGATTGCGGAACTCGCGCTTCCAGCAAAGAGCACACTTGCGCACAGCGCACGTCGTTCCCCAAAATGGCAGTGCGTACCACCCGTTGCTCGGAACCCAGCACGACGGTCCCGGGCTTCTGGCCGGCaaaa                     >  1:643560/5‑226 (MQ=255)
                                                                                                                                                                                                       tGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCCTAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGCTTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACCAACTGGCAGGGCGTAACACCCGTCGCTGAGCCTCCAGCACGCTGGCTCCCTGCCCCAGTCCGGCaacaa                    >  1:535020/1‑225 (MQ=255)
                                                                                                                                                                                                       tGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCcgtcg                                                                                                 <  2:268590/152‑1 (MQ=255)
                                                                                                                                                                                                        ttccattttccgaaaTATCTCCGGATCATGCCGCTGAACCCCTAGCAAATCAGGCGCAGTAAAGCCGGGCATCACCAGGTTGCCGGTAATCCCCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATTGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATAcc                   <  1:330447/215‑1 (MQ=255)
                                                                                                                                                                                                         tCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCCCAAAATGGCAGGGCGTAACACCCGTTCCTGCGCATCCACCAAGGTAGCTCCGGGCATCTGCCCGCCAAAAcca                  >  1:110993/1‑227 (MQ=255)
                                                                                                                                                                                                          ggttttaccGGGATACCCCCGGCACGTGCCGCTGCACCCACAGCAATTTGGGCGCAGTACATCCGGGACTCCTCAGGTTGCCGGTAATCCCCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGAACAGCGCCCGTCGTTCCACCAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCa                 <  1:525962/220‑1 (MQ=255)
                                                                                                                                                                                                          cGATTGGACGGAAATTCCCCGGCTCCTGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCTTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCa                 <  2:671067/228‑1 (MQ=255)
                                                                                                                                                                                                           gATTTGGCGGAATTTCCCCGGCCCATTCCGCTGAACCCATACAAATTTAGGCGACGTAAATCCGGGCACCATCAGGTTGCGCGTAATCCCCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCaa                <  1:621837/228‑1 (MQ=255)
                                                                                                                                                                                                           caTTTGACGGAATATCCCCGGCCACTGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCCCCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCATAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCaa                <  2:13841/227‑1 (MQ=255)
                                                                                                                                                                                                            aTTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAACATGGCAGGGCGTACCACCCGTTGCTGACCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCaa                >  2:601350/1‑227 (MQ=255)
                                                                                                                                                                                                              ttGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGATGCTCCGTGCATCTGGCCGGCAATACCaaa                >  2:467780/1‑225 (MQ=255)
                                                                                                                                                                                                              ttGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCACCAAGGTTGCTCCGTCCATCTGGCCGGCATTACCCAATg              >  2:225181/1‑227 (MQ=255)
                                                                                                                                                                                                                        ttctcCGGCTCATGCCGCTAAACCCATAGCAATTTAGGCGCAGGAAATCCGGTCATCATCAGGTTGCCGGTAATCCCCCGCGATTGCGGAACTCGCTCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCAATGCTTTAACGTc    <  2:404453/226‑1 (MQ=255)
                                                                                                                                                                                                                         cctcCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCAATGCTTTAACGTCCt  <  2:613599/227‑1 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                   
TAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGACGGAATATCTCCGGCTCATGCCGCTGAACCCATAGCAATTTAGGCGCAGTAAATCCGGGCATCATCAGGTTGCCGGTAATCACCCGCGATTGCGGAACTCGCGCTTCCAGCAAAGTGCACTCTTGCGCACAGCGCCCGTCGTTCCACAAAATGGCAGGGCGTAACACCCGTTGCTGAGCATCCAGCAAGGTTGCTCCGTGCATCTGGCCGGCAATACCCAATGCTTTAACGTCCT  >  NC_000913/3728732‑3729174

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