Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,880,380 A→C 5.7% L281* (TTA→TGA)  recF ← gap repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,880,3800AC 5.7% 208.8 / 2.8 70L281* (TTA→TGA) recFgap repair protein
Reads supporting (aligned to +/- strand):  ref base A (33/33);  new base C (2/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

CTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCA  >  NC_000913/3880086‑3880680
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                            
ttttGATATTGTAGGTGTCATAAGAATTGGACATCATCGTTTTTCGCTCTTTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGgcgc                                                                                                                                                                                                                                                                                                        <  1:573370/300‑1 (MQ=255)
   tgttACTGGAGGGGTCATAAGAATTCGACATCAACGTTTCTCGTTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCTTACGTCTGGCGcaag                                                                                                                                                                                                                                                                                                     <  1:223161/298‑3 (MQ=255)
      tACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATAATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTAGCGCCAGACGTAAGGCCcacagcc                                                                                                                                                                                                                                                                                                  >  1:142336/1‑298 (MQ=255)
          ggaggaGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACac                                                                                                                                                                                                                                                                                              >  1:224013/1‑300 (MQ=255)
                   aTAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCAGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCACTCAGGTAGAGACACCGCCGCCCGCTTTAACGGGTGAGGACCTCTCCTTGCGCAAGACGTAAGGCAAACATACAAAGCATAAgaag                                                                                                                                                                                                                                                                                     >  2:627652/1‑298 (MQ=255)
                           tCGACATCAACGTTTCTCGTTCATTTATACTTGGGTTAATCCGTTATTTTTCCCTTTTCCACGGTAAACATCTTCGAATTTGCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCAcgcg                                                                                                                                                                                                                                                                             <  1:169893/301‑1 (MQ=255)
                                 tCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCCACAGCTAAAGCTGCCCACACGCAAAgg                                                                                                                                                                                                                                                                       >  1:623337/1‑301 (MQ=255)
                                  cAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCATTCACGGGTGAGGAACTCTCCTTGAGCCAGACGGAAGGCGCACATCAAAAGCTTAAGCTGCCCACCCAATAAggg                                                                                                                                                                                                                                                                      >  1:244687/1‑300 (MQ=255)
                                                  ttttTACTTTGTTTAATCCTTTATTTTTCCCTTTTCCACGGTATACATCTTCGAATTTTCGGCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGAAAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCGCGCGACGCTCATCATCAAGCTCAGTGGCAAAATCATCTATCGGGTAGAGACACCGCCGCCCGCTATCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATGTTCCACCGGCGCa                                                                                                                                                                                                                                                      <  2:204539/301‑1 (MQ=255)
                                                   ttttgtgggggTTGATCGGTTATTTTACCTGTTTCTACGGGAATCATCTTCGAATTTTCGTCCGACATGTATATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAGCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCAc                                                                                                                                                                                                                                                     <  1:281295/293‑1 (MQ=255)
                                                     ataCTTGGGTTAATCCGTTATTTTACTCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTCCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCAccc                                                                                                                                                                                                                                                   >  2:263416/1‑300 (MQ=255)
                                                        cTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTAAAGCTGCCCACGCGATAAGGTATCTTACACCGGCGCACCGtag                                                                                                                                                                                                                                                >  2:143981/1‑299 (MQ=255)
                                                          cagacgtgtgctcttccgatctcTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGgc                                                                                                                                                                                                                                              <  2:312050/279‑1 (MQ=255)
                                                            ggTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGACAAAAACATCAAACAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGCGGAACTCCCATTGAGCCAGACGTAAGGCGCCACTCCACAGCTTAAGCTGCCCACGCGAAAAGGTATCTACCACCGGCCCACCgaccgccc                                                                                                                                                                                                                                            >  1:478309/1‑294 (MQ=255)
                                                                   tctttATTTTTCCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGCTCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTACGCGGCTGGCAAGCAGACCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTGAAGCTGCCCACGCGATAAGGTATCGTCCACCGGCGCACCGTCGGCGCGAATGCg                                                                                                                                                                                                                                     <  1:156055/298‑1 (MQ=255)
                                                                    ggTTTTTTTACCCTTTTCTCGGGTAAACTTCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGt                                                                                                                                                                                                                                    <  1:138453/300‑1 (MQ=255)
                                                                    cGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCga                                                                                                                                                                                                                                    >  2:156206/1‑300 (MQ=255)
                                                                           ttACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGCGGAACTCTCCTTGCGCCAGACGTAAGGCGCACCTCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTACACCGGCGCACCGTCGGCCCGAAAGCGCAAGTcac                                                                                                                                                                                                                             >  1:519941/1‑299 (MQ=255)
                                                                              ccttttgtcACGGGAAACATGTTCGAAGTTTCGTGCGACAGTTCGATAACGGGTGCAGCACTGATCGCGCTGACAAAGACCGGTGATGGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAGTGCGTAAGTCCGCtt                                                                                                                                                                                                                          <  1:144693/293‑1 (MQ=255)
                                                                                cTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGgcagatcggaagagcgtcgtgacg                                                                                                                                                                                                                        >  1:311646/1‑279 (MQ=255)
                                                                                   ttCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGc                                                                                                                                                                                                                     <  2:478910/301‑1 (MQ=255)
                                                                                      cACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCggg                                                                                                                                                                                                                  >  2:351204/1‑301 (MQ=255)
                                                                                                      gAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAgg                                                                                                                                                                                                  >  2:314211/1‑301 (MQ=255)
                                                                                                        aTTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGta                                                                                                                                                                                                >  2:96724/1‑300 (MQ=255)
                                                                                                          tttCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTCAg                                                                                                                                                                                              >  1:109664/1‑301 (MQ=255)
                                                                                                                ccGACGTGTCTATACCTTGGTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGAAAGCAGCCCGCGACGCGCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGTGTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCg                                                                                                                                                                                        <  1:350741/301‑1 (MQ=255)
                                                                                                                 cGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGGTAGCTGGCGa                                                                                                                                                                                       >  2:347321/1‑301 (MQ=255)
                                                                                                                         atgtaatgtgttcggtccgggtgcgtcTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTc                                                                                                                                                                               <  1:474314/275‑1 (MQ=255)
                                                                                                                          tataACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCa                                                                                                                                                                              >  2:267813/1‑301 (MQ=255)
                                                                                                                           ataACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCaa                                                                                                                                                                             >  2:35878/1‑301 (MQ=255)
                                                                                                                               cGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAAAGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAAtt                                                                                                                                                                         <  2:244954/301‑1 (MQ=255)
                                                                                                                                   gtcagcgcgatcagtgcTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTACGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGGGTCAGTTAGCTGGCCATCGCGTTCCAAATTACGt                                                                                                                                                                     >  1:129219/16‑301 (MQ=255)
                                                                                                                                    tCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGtt                                                                                                                                                                    <  2:224256/301‑1 (MQ=255)
                                                                                                                                       gCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCGAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCa                                                                                                                                                                 >  2:387676/1‑301 (MQ=255)
                                                                                                                                        tgttcagacgtgtgctcttccgatctgATTGCTTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAg                                                                                                                                                                <  2:324579/276‑1 (MQ=255)
                                                                                                                                            gATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCAcc                                                                                                                                                            <  2:189318/301‑1 (MQ=255)
                                                                                                                                                       aaaggACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTc                                                                                                                                                 <  1:35827/297‑1 (MQ=255)
                                                                                                                                                               tgtggtgttgTCGCTTTTAATCGGCTTGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGTCAAAATCATCTAGCAGGTAGAGACACCGCCGCCCGCTTTCACTGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACCTCAACAGCTTCAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCtt                                                                                                                                         <  1:626775/292‑1 (MQ=255)
                                                                                                                                                                gtgATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGAACCTCAGCATATTCTGTCTCttt                                                                                                                                        >  2:523108/1‑301 (MQ=255)
                                                                                                                                                                  gATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATACGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGGGCGCGGTGTAGGATAGCCAGCAACCGCGATCCAAATTACGTTTCAgaaatcaggagagcgtcgtggtgggga                                                                                                                                      >  1:324175/1‑275 (MQ=255)
                                                                                                                                                                     tgctttgCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTccc                                                                                                                                   <  2:538221/295‑1 (MQ=255)
                                                                                                                                                                            cTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCCAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCcgcc                                                                                                                            >  1:440148/1‑300 (MQ=255)
                                                                                                                                                                                  aGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGaa                                                                                                                      <  2:520613/301‑1 (MQ=255)
                                                                                                                                                                                           cAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGCCCGTGCGCGGTGTAGGTTAGCTGGCGAGCGCGTTCAACATTACGTTCAACCACCTCAGCATATTCTGTCCTTTTCTCCCcgccgcccccgaccgagcaaat                                                                                                             >  1:328290/1‑285 (MQ=255)
                                                                                                                                                                                                agcCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCagag                                                                                                        <  2:109788/301‑1 (MQ=255)
                                                                                                                                                                                                     gcgACTCTTATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCACACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAAc                                                                                                   <  2:49627/301‑1 (MQ=255)
                                                                                                                                                                                                        ccGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCa                                                                                                <  2:142461/300‑1 (MQ=255)
                                                                                                                                                                                                          ccttattattcaGCTTAGAGGCAAAATCATTTATTAGGGTGAGACACCGGCGCCCGCTTTCACGGGGGAGGAAATCTTCTTGGGGCAGAAGTAAGGCGGAAATTAACAGCTTAAGCCGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAgg                                                                                              <  2:578441/290‑1 (MQ=255)
                                                                                                                                                                                                               tcatcaAGCTTATGGGCAAAAGCATCTAGCCGGTAGAGACTCCGACGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGaa                                                                                         <  1:96600/301‑1 (MQ=255)
                                                                                                                                                                                                                    ttccgatctggCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGc                                                                                    <  2:138889/292‑1 (MQ=255)
                                                                                                                                                                                                                         cAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCCCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGACAACTCAGGGAGACATTGcatcca                                                                               >  1:163019/1‑296 (MQ=255)
                                                                                                                                                                                                                             ggCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAAAAATTGcagatcggaa                                                                           >  1:138764/1‑292 (MQ=255)
                                                                                                                                                                                                                                  aaTCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCCACCCAGCATATTCTGCCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGACAAAACTCAGGAAGAAATTGCTTACAGATAACacc                                                                      >  1:561694/1‑299 (MQ=255)
                                                                                                                                                                                                                                  aaTCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGc                                                                      <  1:263148/301‑1 (MQ=255)
                                                                                                                                                                                                                                              ggTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCcgg                                                          >  2:4062/1‑300 (MQ=255)
                                                                                                                                                                                                                                                      cacCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCgg                                                  >  2:116570/1‑301 (MQ=255)
                                                                                                                                                                                                                                                            cgcccgcTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCtat                                            <  2:328694/301‑1 (MQ=255)
                                                                                                                                                                                                                                                             gcccgcTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTCCGCGGGGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCCGCATCTTCTGTCACTTTCTCCCAGCCGCGCTGGACAGGGAAAGTCCGAGAAAACACAGGGAGAAATTGCTAACAGGTATCCTCCATATCGGCCGCGATCACGGCGCtaaa                                           >  1:588041/1‑299 (MQ=255)
                                                                                                                                                                                                                                                              cccggTTTCACGGGTGCGGAACTTTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATAc                                          <  1:4062/301‑1 (MQ=255)
                                                                                                                                                                                                                                                              cccgcTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATAc                                          >  2:389256/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                  cttgtacgggggtgggtctttcctttgcccAGACTTAAGGCGCACTGCAACAGCTTCAGCTGCCCACGCGAGAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCcg                                      <  1:313807/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gggTGGGGAACTCTCCTGGGGCCAGAGGTAAGGCGCACCGCAACAGGTGAAGCTGCCCACGCGGTAAGGTATCTTCCACCGGCGGACCGTCGGGGGGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATGCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCa                               <  1:522436/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                            tGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCCTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAggg                            >  2:219694/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                             gatgAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTg                           <  2:624196/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                              tttccgatcttgtGCGCCTTACGTCTGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGc                          <  2:362923/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                               ggaatcctgtgtGCGCCAGCCGTAAGGTGAACATCATCAGCTTCTGCGGCCCACGCTGTAAGTGTGCTAGCACGGGGGCACCGTCGTGGCGAATGCGTATGTCCGCTTTGTGCGGGCCGTGCGCGTTGAAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATTTTCTGTCGCTTGCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCACGCGCCAGGTGCt                         <  1:267545/290‑1 (MQ=255)
                                                                                                                                                                                                                                                                                gctgttgatgtGCGCCTTACGTCTGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTg                        >  2:11317/11‑301 (MQ=255)
                                                                                                                                                                                                                                                                                   ctctCCTTGCGCCAGACGTAAGGCGCACATCAACAGCGTAAGCTGCCCACGCGATAAGGTATCTGCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGAtt                     <  1:116446/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       ccTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCt                 >  2:126704/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                        tgtGCGCCTTACGTCTGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGCAATTGCTTACAGGTATCAGCCATATCGGCAGCGATACCGGCCCTATACTCCGCGCACCAGgcgcagcacggaag                >  1:362439/3‑291 (MQ=255)
                                                                                                                                                                                                                                                                                                      aaGTCGCACATCCACATGTTAAGCTGCCCACGCTATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCa  <  1:388742/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                            
CTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCA  >  NC_000913/3880086‑3880680

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