Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,142,657 C→G 16.7% G43G (GGC→GGG frwC → putative enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,142,6570CG16.7% 202.1 / 6.9 86G43G (GGC→GGGfrwCputative enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base C (14/56);  new base G (14/2);  total (28/58)
Fisher's exact test for biased strand distribution p-value = 7.09e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATCGGCGGCATCGTGGTGCATTACCTGAAGAAA  >  NC_000913/4142420‑4142894
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
aaaGTCCATTAAATGGCCTTTTTTTCCCCTGTTTGCTCCGCCTTGGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGgcg                                                                                                                                                                                                                                                        <  2:789675/239‑1 (MQ=255)
aaaGTCCATTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGgcg                                                                                                                                                                                                                                                        <  2:835620/239‑1 (MQ=255)
        tttaaTGGCTTTTTTTTCCCTTGTTTTCCCCCGCTGCGATTGCCTTTTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCg                                                                                                                                                                                                                                               <  2:271413/236‑1 (MQ=255)
         tAAATGGCCTGGTTATCCACTGTTTGCTCCGGCTGCGATTGCCTATTGTTTCCTCCATCAAAAAGATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCgg                                                                                                                                                                                                                                              <  1:613787/239‑1 (MQ=255)
          aattGGCCTTTTTTACCCTTGTTTGCCCCGCCTGCGATTGCCTATTGTTCCCCCCATCAAAAATATGGGCCTGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGa                                                                                                                                                                                                                                             <  2:76262/236‑1 (MQ=255)
          aaaTGGCCTTTTTTTCCCTTTTTTGCTCCGCCTGCGTTTGCTTTTTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGa                                                                                                                                                                                                                                             <  2:465002/239‑1 (MQ=255)
          aaaTGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCgg                                                                                                                                                                                                                                              <  1:881509/238‑1 (MQ=255)
          aaaTGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCgg                                                                                                                                                                                                                                              <  1:556073/238‑1 (MQ=255)
              ggCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATAGGGCGGGATTTTGCTGGCGGGTTTCGGTAATGTGGTATGGAAAGGGGCacgggcgggaggc                                                                                                                                                                                                                                          >  2:505422/1‑227 (MQ=255)
                       atccagttgttGCTCCGCCTGCGATTGCCTATTGTTTCCCCCATCAAAAATATGGACATTGGCTTTCGCAAGCCCGGTAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGAt                                                                                                                                                                                                                             <  1:1054563/230‑1 (MQ=255)
                            cTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTTCGTCATGTTGTATGGCCAAGGGGCAGTGCCGGGTGGCGTAACCGGACCaaa                                                                                                                                                                                                                      >  1:640353/1‑239 (MQ=255)
                            cTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGCTTTCGTCTTTTTTGATGGGCAAGGGGCGGGGTCGGATGCCGTAGCCGGTCCcaa                                                                                                                                                                                                                      >  1:477208/1‑239 (MQ=255)
                              tttttccccGCCTGCGATTTCCTTTTTTTTCCTCCATAAAAAATTGGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATc                                                                                                                                                                                                                    <  2:541791/232‑1 (MQ=255)
                              gTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAACAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACCTGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGGTTCCGCTCAGTTTGTTTGGAAATGCGGCAGTCCGGGATGCCGAAGCCGACCCaaac                                                                                                                                                                                                                     >  2:762043/1‑238 (MQ=255)
                               tttGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGGTTCCGTCCTGTTGTATGGGAAAGGGGCCGTGCCGGATGCCGTAGCCGCACCAAATCt                                                                                                                                                                                                                  >  1:138697/1‑239 (MQ=255)
                                  tccccGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACGTTGGCTTTCGCAAGCCCGGGAACAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGaa                                                                                                                                                                                                               <  1:505422/236‑1 (MQ=255)
                                       gCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTTCGTCACGTTGTATGGCCAAGGCGCAGTGCCGGGTGCCGTAGCCGACCCCAAACTTGaaaaaa                                                                                                                                                                                                          >  1:50101/1‑239 (MQ=255)
                                                  ccTTTTTTTTCCTCCATAAAAAATTGGGGCATGGGTTTTCGAAAGCCCGGGAGAAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAAGCCCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATc                                                                                                                                                                                              <  2:374145/239‑1 (MQ=255)
                                                  acTATTGGTTCCTCCATCAAAAATATGGACCTGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAAGGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGaaaaaa                                                                                                                                                                                                          <  1:509697/226‑1 (MQ=255)
                                                         tttCCTCCATTAAAAATTTGGGCATGGGGTTTCGCAAACCCGGGCGCAGCTTATGAAGGAGTGGGTGCAGATCCTGAAAATCTCCCGTCAGCATTTAATGCCGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTg                                                                                                                                                                                       <  2:439345/239‑1 (MQ=255)
                                                          ttCCTCCATAAAAAATATGGCAATGGGCTTTGGCAAGCCCGGGAGCACGTTTTGAATGAGTTGGGGCAGACCCTGAAAAATCCCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGc                                                                                                                                                                                      <  2:82404/239‑1 (MQ=255)
                                                           tcctccATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGGTTTCGGCCATTTTTATGGCACAAGGGGAGTGGCGGGTGCCCGTGGCGATCCCAATATTAAAAAAAAGtttttaagatggggtggg                                                                                                                                                                                     >  2:1050847/1‑224 (MQ=255)
                                                           tcctccATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCCGCCTTTTGTATGGCTAAGGGGCAGTGCCGGGTGCCGTAGCCGATCCAAATCTTAAAAAAACTTTTTGattccgggtgtcg                                                                                                                                                                                      >  1:411405/1‑228 (MQ=255)
                                                           tcctccATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGAATCGGGCGGTATTTTGCTGGCGGTTTTCTCTCTGTTTTAGTGGAAAAGCGGCGGGGCCGGTGCCGTGGCCGCTCACAAATTTGAAAAAAATGttcaatttggggggtgg                                                                                                                                                                                      >  2:116076/1‑224 (MQ=255)
                                                            cctccATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTTCGTCCTTTTTTTTGGGCAAGGCTGAGGTGCGGATGGCGTAACCGGGCCCAAACTTGAAAAAATTTTTTTATTTCgcgtttgag                                                                                                                                                                                     >  1:885533/1‑231 (MQ=255)
                                                                cATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATTTTGTTTGGCAAAGGGGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCtt                                                                                                                                                                                >  1:340330/1‑239 (MQ=255)
                                                                   cAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTAGTTTGCTGGCGGTTTACGTCATGTTGGTTGGAAAAGGGGCAGTGCCGGAGGGCGCAGGCGATCCAAACATGTAAAAACCATTTGATATCGGAGTTGGGGGCGTGAc                                                                                                                                                                             >  1:505302/1‑239 (MQ=255)
                                                                    aaaaaTTTGGGCATTGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGATTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACg                                                                                                                                                                            <  2:1081847/239‑1 (MQ=255)
                                                                                   ggCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTTCGTCCTTGTTTTTGGCAAAGGGGCAGGGCCGGGGGCCCGTGCCGGTCCAAAATCGGAAAAAACTTTTGTTAACGGCGGTGggggggttagcgggggtggggcctt                                                                                                                                                              >  2:824273/1‑217 (MQ=255)
                                                                                     ttttCGAAAGCCGGGGAGCTGCTTTTGAATGAGTTGGTTCAGATCCTGAAAAATTCCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGGGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTcc                                                                                                                                                            <  2:796201/237‑1 (MQ=255)
                                                                                       ttCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCCTGTTGTTTGGCAAAGGCGCCGGGCCGGATGCCGTAGCCGGTCCACATCTTCAAAAAATGTTTGGTATCTGGGTTTCGGGGtttgcctgggtggggccttttcctc                                                                                                                                                         >  1:1019140/1‑216 (MQ=255)
                                                                                       ttCCCAAGCCCGGGGGCAGCTTATGAATGAGTTGGTGCAGTTCCTGAAAAATCCCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCg                                                                                                                                                         <  2:477208/239‑1 (MQ=255)
                                                                                        tCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTTCGGCCTGTTGTATGGCAAAGGGGCCGTGCCCGGTGCCGTACGCCGTCCACATCTGAAAAAACTTGTTGTTATCGGGGTTTCGGGGTTTGCGCTGATGGGGCCTTTtctccc                                                                                                                                                        >  1:139202/1‑237 (MQ=255)
                                                                                         gcaaaGCCCGGGAGCAGGTTTGGAATGGGTTGGTGCCGGTCCTGAAAAATTCCCGTCAGCATTTAATGCCGGGCGTTGCACACATGATTCCCTTCGTGGTTTCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCa                                                                                                                                                       <  2:623649/236‑1 (MQ=255)
                                                                                                ccggggggtgtttattaATGGGTTGGTGCAGATCCTGAAAAATCCCCGTCAGCATTTAAGGACGGGCTTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACa                                                                                                                                               <  2:546236/228‑1 (MQ=255)
                                                                                                ccgggacccttttttgattggGTTGGTTCTGTTCCTGAAAATTACCCTGATGCATTGAATGCCGGGCGTTTCACACATTATTCCCTTCGTGGTGTCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACa                                                                                                                                               <  2:485920/225‑1 (MQ=255)
                                                                                                ccGGGAGCAGTTTTTGAATGGGTTGGTGCAGATCCTGAAAAATCCCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACa                                                                                                                                               <  2:340330/239‑1 (MQ=255)
                                                                                                ccGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAGTACCCGTCAGCATTTAATGACGGGCGTTTCACACTTGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACa                                                                                                                                               <  1:76984/239‑1 (MQ=255)
                                                                                                ccGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATTCCCGTCAGCATTTAATGGCGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACa                                                                                                                                               <  1:133546/239‑1 (MQ=255)
                                                                                                      gcagttttttaattggTTGGGCCAGTTCCTGAAAAATGCCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGt                                                                                                                                          <  2:885533/223‑1 (MQ=255)
                                                                                                      gcagtTTTTGAATGAGTTGGTGCAGATCCTGAAAAATCCCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGt                                                                                                                                          <  2:1073446/238‑1 (MQ=255)
                                                                                                      gcagtTTATGAATGAGTTTGTTCAGATCCTTAAAAATTCCCGTCCGCATTTAAGGGCGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATTTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGt                                                                                                                                          <  2:505302/238‑1 (MQ=255)
                                                                                                      gaaggttATGAATGAGTTGTGGCAGTTCCGGAAAAATGCCCGGCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGt                                                                                                                                          <  1:1050847/233‑1 (MQ=255)
                                                                                                           ttttgaatgaGTTGGTGCAGATCCTGAAAAATCCCCGTCAACATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTc                                                                                                                                     <  2:681870/235‑1 (MQ=255)
                                                                                                           ttATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTAtt                                                                                                                                      <  1:824273/237‑1 (MQ=255)
                                                                                                            tATTAATGGGTTGGGTCAGAACCTAAAAAATTCCCGTCAGGATTTAATGACGGGCGTTTCCCACATGATTCCCTTCGGGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCa                                                                                                                                   <  2:640353/239‑1 (MQ=255)
                                                                                                            tATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCa                                                                                                                                   <  1:33867/239‑1 (MQ=255)
                                                                                                                aatgaGTTGGTGCAGATCCTGAAAAATCCCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGc                                                                                                                                                                                      <  1:243707/185‑1 (MQ=255)
                                                                                                                   ggtgtggtggCGATCCTGAAAAATTCCCGGCAGCCTTTAAGGGCCGGCGTTTCCCCACGGATTCCCGTCGCGGTATCGGGCGGTTGTTTGCTGGCGGTTGCCGTCAGGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCaga                                                                                                                            <  2:689553/235‑1 (MQ=255)
                                                                                                                    aGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGGGCCGGGTGCCGTAGCCGATACAAATCTGCAAAAACTGTTTGTTATCGGCGTTGCGGGCTTTGCGCTGCTGGTGCCCTTCCTCGCCGCGTACATCGGTTATTCCCTTTCcgag                                                                                                                           >  1:72655/1‑239 (MQ=255)
                                                                                                                     gTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGGGCAGGGCCGGGGGCCGGAGCCGGTCCCAAACTGAAAAAACTGTTTGATATCCGCGTTGCGGGCTTGTGGCTGATGGTGCCTTTCCTCCCCGCGGACATCGGTGATTTCAtttctaggc                                                                                                                          >  1:321852/1‑232 (MQ=255)
                                                                                                                                   tGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGATTTCGTCATGTTTTATTGGCAAGGCGCACTGCCCGGATCCCTAGGCGACCCAAATATGAACAAAACTGTTTGTAGTGGCGTTTCGGGGTTTTCGGCGCTGGGGCCCTTCTTCCCACCGTACCACGGGTATTCCCTTTCCGAACGttttcggccggc                                                                                                             >  2:688821/1‑228 (MQ=255)
                                                                                                                                      aaatttcccctccgtttttAATGACGGGCGTTTCACAACTGATTCCCTTGGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGGGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCg                                                                                                         <  2:349330/223‑1 (MQ=255)
                                                                                                                                      aaaaTTCCCGTCAGCATTTAATGGCGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCg                                                                                                         <  1:881254/239‑1 (MQ=255)
                                                                                                                                               cgttcggatttaattgggggGTTTTACATCATGATTTCTGTCGGGGTATCGGGCGGGATTTGGGTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCGGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCGCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTAt                                                                                                 <  2:513234/223‑1 (MQ=255)
                                                                                                                                                      tttAATGGCGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCt                                                                                         <  1:204918/239‑1 (MQ=255)
                                                                                                                                                      tttAATGCCGGGGTTTTCACACAGTTTTCCTTTGGTGGTATCGGGCGGTATTTGGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCt                                                                                         <  2:1019140/239‑1 (MQ=255)
                                                                                                                                                               gggCTTTGCACACATGATTCCCTTCGTGGTATCGGGCTGTATTTTGCTGGCGGTTTCCGGCATGGTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTGCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGt                                                                                <  1:903247/238‑1 (MQ=255)
                                                                                                                                                                 gctttttccacaTTTTTCCCTTCGTGGTATCGGGCGGTTTTTTGCTGGCGGTTTCCGTCATGTTTGATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTaa                                                                              <  2:138697/230‑1 (MQ=255)
                                                                                                                                                                 gCGTTTCCCACATGTTCCCCGTCGTGTGATGGGGTGGTTTTTGTCTGGCGGGTTCGTTCATGTTGTTAGGCACATGCGAAGTGCCGGATTCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGTTGCCGTGCCTGGCAGCGTACAGCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTaa                                                                              <  2:968973/238‑1 (MQ=255)
                                                                                                                                                                             tggATTCCCGTCGGGGTACCGGTCGGTATTGTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATAGCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTg                                                                   <  1:326070/236‑1 (MQ=255)
                                                                                                                                                                             aTGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGcc                                                                                                                                                                 <  1:331435/146‑1 (MQ=255)
                                                                                                                                                                              tttttCCCTTCGTGGTTTCGGGCGGTTTTTTGCTGGCGGTTTCCGTCATGTTGTATGGTAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCg                                                                 <  2:139202/236‑1 (MQ=255)
                                                                                                                                                                              tgtttCCCTTCGTGGTTTCGGGGGTTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCg                                                                 <  2:321852/236‑1 (MQ=255)
                                                                                                                                                                               gATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATTGCAAAGGCGCAGTGCCCGAAGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCTCAGGGTACATCGGGTATTCCCTTGCAGAACGTTCTGCGCTGGCTCCCTGCGCTATCGGGGCCTGGGTTGGGAACAGCTttgttggtg                                                                >  1:237210/1‑236 (MQ=255)
                                                                                                                                                                                  tCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCaga                                                                                                                            <  1:348654/177‑1 (MQ=255)
                                                                                                                                                                                   cccTGCGGGGTATCGGGCGGTAGTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGctt                                                            <  1:366082/239‑1 (MQ=255)
                                                                                                                                                                                     cTTGGTGGTTGCGGGGGGTTTTTTGGTGGCGGTTTCCGTCAGGTGGTATGGCAAAGGCGCAGGGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGcttct                                                          <  2:72655/239‑1 (MQ=255)
                                                                                                                                                                                      ttCGTGGTTTCGGGGGTTTTTTTGCGGGCGGTTTCCGTCATTTTGTAGGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGcttctt                                                         <  2:493954/239‑1 (MQ=255)
                                                                                                                                                                                      ttCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGGCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGTTCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGcttctt                                                         <  1:417592/239‑1 (MQ=255)
                                                                                                                                                                                       tCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATTTTTTATGGGAAAGGGGCGGGGCCGGATGCCGTAGCCGCGCCAAAACTTAAAAAACTGGTTTTTAACCGGGTTTCGGGGTTGGCGCCGATGGTGCCTTTCCTCCCCGCGGCACCCCGGGATTTCCTTTGCGCGCGTTCTGCTGTGGGTCCCTGCGCTATTGGGGGCTCGGGGGGGGAACAACTTTGGGGGGGGgttttt                                                         >  2:937579/1‑239 (MQ=255)
                                                                                                                                                                                               tCGGGCGGTATTTTGCTGGCGGTTTCCCTCATGTTTTATGGCAAAGGCGCAGTGCCGGATGCCGGAGCCCATCCCAATCTGAAAAAAACGTTTGGTATCGGCGGTGCGGGCTTGACCCGGATGGGGCCTTTCCTCGCGACGTAACTCGGGTGTTTCCATGCCGAGGGTTTTTCGCTGGCTCCGCGGGCTCTCGGGGCCCGGGTGGGTAAAAGCGTTGGTGCGGGGTGCTTTggggg                                                   >  2:264895/1‑233 (MQ=255)
                                                                                                                                                                                                 gggCCGTATTTTGCTGGCGTTTTCCGTCATGTTGTAAGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAGTCTTGAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGa                                              <  1:535693/239‑1 (MQ=255)
                                                                                                                                                                                                       tATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCa                                                                                                                              <  1:597190/154‑1 (MQ=255)
                                                                                                                                                                                                                  cGGTTTCCGTTATGTTGTATGGCAAAGGCGCAGGCCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATcgg                              <  1:1040138/238‑1 (MQ=255)
                                                                                                                                                                                                                   gggTTTCCGTCTTTTTGTTTGGAAAAGGGGCCGTGCCGGATGCCGGAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATcgg                              <  2:483344/237‑1 (MQ=255)
                                                                                                                                                                                                                   cgttttCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTTTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATcgg                              <  1:916674/235‑1 (MQ=255)
                                                                                                                                                                                                                       ttCCGTCATGTTGTATGGCAAAGGGGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATCGGCGGCa                         <  1:883930/239‑1 (MQ=255)
                                                                                                                                                                                                                       ttCCGTCATGTTGGATGGCAAAGGCGCAGTGCCGGATGCCGTAGACGATCCAAATCTGAAAAAACTTTTTGATATCGGCGTTGCGGGCTTGCCGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATCGGCGGCa                         <  1:721426/239‑1 (MQ=255)
                                                                                                                                                                                                                          cGTCTTGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTGGCCGATCCAAATCTGAAAAAACTGTTTGATTTCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATCGGCGGCATCg                      <  2:942671/239‑1 (MQ=255)
                                                                                                                                                                                                                              tcttttgttagggaaagggCGCGGTTCCGGAGGCCGTGGCCGAACCAAATCGTAAAAACCTGTTGGTTGTCGGCGGTGCGGGCTTGCCGCTGATGGTGCCTTTCCTGGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGGGCGCTAGCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATCGGCGGCATCgtc                    <  2:696423/222‑2 (MQ=255)
                                                                                                                                                                                                                                     tATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATCGGCGGCATCGTGGTGCAtt             >  1:785831/1‑238 (MQ=255)
                                                                                                                                                                                                                                        tGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTTTCGGCGGCACCGGGGGGCATTccc          >  2:1046209/1‑239 (MQ=255)
                                                                                                                                                                                                                                           cAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATGGGTGCCTGGGTTGGTAAAAGATTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATCGGCGGCATCGTGGGGGATTACCTga       >  1:506006/1‑239 (MQ=255)
                                                                                                                                                                                                                                             aaGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCCCTCGTCGCCGGGATTATCGGCGGCCTCGTGGGGCCTTACCCgaag     >  2:884345/1‑239 (MQ=255)
                                                                                                                                                                                                                                              aGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTTGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCCGGCTTCTTTGGTGCACGTTTCGCCGGGTTTATCGGCGGCTTCGTGGTGCTATTCCTgaaaa    >  2:159010/1‑237 (MQ=255)
                                                                                                                                                                                                                                              aGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATCGGCGGGATCGGGGGGCATAACCGgaaga    >  2:379034/1‑239 (MQ=255)
                                                                                                                                                                                                                                              aGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGAAAACGCTTTGGTGGGGGCTTCTTTGGTGCCCTGATCGCCGGGATTATCGGGGGCATCGGGGGGCATTACCTgagga    >  2:1052438/1‑239 (MQ=255)
                                                                                                                                                                                                                                                gcgcAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATCGGCGGCATCGTGGTGCATTACCTGAAGaaa  >  1:748340/1‑239 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
AAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATCGGCGGCATCGTGGTGCATTACCTGAAGAAA  >  NC_000913/4142420‑4142894

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: