Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,446,970 T→G 18.8% V953G (GTG→GGG)  tamB → translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,446,9700TG18.8% 128.5 / 3.7 78V953G (GTG→GGG) tamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (26/34);  new base G (11/7);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGATGTGCAAAGCCCGCAGTTGTTTGGTCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTC  >  NC_000913/4446733‑4447208
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
gggttgtgAAAAGCCCGGTTTTGTTTGGGCAGTTTCCGGGTAGGGGTGGGGATATGGATGGCAACGTTATGCCGTTTGATATGACACCGAGCCAGCTTGCGGTAAACGTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTg                                                                                                                                                                                                                                                     <  1:78651/235‑1 (MQ=255)
      tgCAAAGCCCGCAGTTGTTTGGTCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGGCCGGGTAACGGCGAAAGGCGGGAAAGGCCGGAt                                                                                                                                                                                                                                                >  2:418156/1‑238 (MQ=255)
      tgCAAAGCCCGCAGTTGTTTGGTCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGGAACGGCGAAAGGCCGTAAAGGGCGGagc                                                                                                                                                                                                                                               >  2:302092/1‑237 (MQ=255)
           aGCCCGCAGTTGTTTGGTCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGAGCAAAGCGGTAAAGGGCGGGTAAccat                                                                                                                                                                                                                                          >  2:442683/1‑237 (MQ=255)
                gCAGTTGTTTGGTCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGGGAAAGGCCGTAAAAGGCGGATCACCGTgccg                                                                                                                                                                                                                                      >  2:174464/1‑238 (MQ=255)
                       tttGGTGTGTTTGAGGTTGCGGGTTGGGATTTCGAGGGCACATTTATGCCGTTTGAATTGAAACCGTGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGCCGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTGCCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGt                                                                                                                                                                                                                              <  1:562946/239‑1 (MQ=255)
                       tttGGTCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGGGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGAGGATCACCGGGCCGCCGAGGGt                                                                                                                                                                                                                              >  2:489619/1‑239 (MQ=255)
                            tCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGGAACGGCGAAAGGCAGTAAAGTGGGGGTCACCGTGCCGCCGATGGTACGaa                                                                                                                                                                                                                         >  2:849588/1‑239 (MQ=255)
                              aGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCCATGGGACGAATg                                                                                                                                                                                                                      >  2:628201/1‑239 (MQ=255)
                                cTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGGAAAGGCGAAAGGCAGTAAAGTGCGGGACACCGTGCCGCCGATGGGACGGATGGa                                                                                                                                                                                                                    >  2:144187/1‑239 (MQ=255)
                                     ggTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAAGGCGGGACACCGTgccccccgacggagggatgggtgta                                                                                                                                                                                                                 >  2:444576/1‑219 (MQ=255)
                                        tgcGGGTGGGGATAACGGTGGCAACATTATGCCGTTTGATATGCAACCGGGCCAACTTGCGGTCAACTTTAACGGGTGGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCg                                                                                                                                                                                                             <  2:73983/236‑1 (MQ=255)
                                         aCGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGGAACGGGGAAAAGCAGTAAAATGGGGAACACCGTGccgccggagggagggaaggggggatcgcc                                                                                                                                                                                                            >  1:374012/1‑216 (MQ=255)
                                          cGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGGAAGGCAGGAAAGTTCGGATCACCTGGACGCCGGTGGGACGatggggtgaattcaca                                                                                                                                                                                                           >  2:470295/1‑224 (MQ=255)
                                          cGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAATGCGGATCACCGGGCCGCCGATGGTACGAATGGGTGTATCGCCa                                                                                                                                                                                                          >  2:452702/1‑239 (MQ=255)
                                             tggtgGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGGGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGa                                                                                                                                                                                                        <  2:488690/236‑1 (MQ=255)
                                                tgGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGGTCACCGTGCCGCCGGTGGTCCGAATGGTGGTATCGCCCGAtgtt                                                                                                                                                                                                   >  2:74788/1‑239 (MQ=255)
                                                 ggATATCGATGTAAACTTTATGCCGTTTGATATGCAACCGAGCCAGGTTGCGGTCAACTTTAACGGTATGCGCGCGACGCTGTCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGAtgtt                                                                                                                                                                                                   <  2:506134/238‑1 (MQ=255)
                                                  gttttctttGGCAACTTTATGCCGTGTGATATGCACCCGGGCCAGGTTGGTGTCAACTTTAACGGTAGGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGAtgttg                                                                                                                                                                                                  <  1:630951/230‑1 (MQ=255)
                                                           ggaaattttATGCCGTTTGATATGCAACCGAGCCAGGTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGc                                                                                                                                                                                       <  2:626674/233‑1 (MQ=255)
                                                           gccaccTTTATGCCGTTTGATTTGAAACCGAGCCAACTTGCGGCAAACTTTAACGGTATGCCCCGGACGCTTGCCGGTAAAGTCCGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGc                                                                                                                                                                                       <  1:235225/234‑1 (MQ=255)
                                                                tttttgtccGTTTGTTAGGAACCCGAGCCAGCTTGCGGTCACTTTGAACGGTATGCGCTCGGCGGTTGCCGGTACAGTACGGACCCAGAAGGGTGAAATCTCCCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTacac                                                                                                                                                                                  <  1:418156/232‑1 (MQ=255)
                                                                cTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGc                                                                                                                                                                                       <  2:592878/234‑1 (MQ=255)
                                                                   tATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCaa                                                                                                                                                                               <  2:233425/239‑1 (MQ=255)
                                                                          tttggattggaaCCGGGCCAGGTTGGGGACAACTTTAAGGGTATGGGCTCGATGCTGGCCGGTACAGTACGGACCCAGGAGGGTGAAAACTACCTGAACGGTGATGCGGACTGGAGCCAAATTGAAAACGGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGtt                                                                                                                                                                        <  1:634049/229‑1 (MQ=255)
                                                                           ttGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCCCCGGTAAAGGCGAAAGACAGAAAAAGGCGGCACCCCGTGCCGCCGGAGGCACGAAGTGGGTTATCCCCACATGTTAGTTTCGGGCCTAAACAAAAcacattt                                                                                                                                                                       >  1:248687/1‑233 (MQ=255)
                                                                                     cccGGGCCAGCTTGCTGTCCAGTTTACCGGTAGGCGCCCGACGCTTGCCGGGCAAGTGCGGACCCGGCAGGGTGAAGTGTGCCTGAACGGTGATGCCGACTGGAACCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATg                                                                                                                                                             <  1:63565/238‑1 (MQ=255)
                                                                                      ccGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGGGATGCCGACTGGAGCCAACTTGAAAACTGGCGGGCGCGGGGAACGGCGAAACGCGGGAAAGTGCGGATTCCCGTGCCGCCGGTGGGACGAACGGGTGTATAGCCCGCTGTTGTATTTGAGGCACTACCAAAACTGTATAaccccccggg                                                                                                                                                            >  1:594265/1‑232 (MQ=255)
                                                                                       cggccccttttGCGTTCAACGTATACTTTAGGCGCTAGACGCTTGCCGGTCCAGTACGGACCCAGCGGGGTAAAATCTCCCCGAACGGGGAGGCCGACTGGAGCCAAATTGAAAACGGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGt                                                                                                                                                           <  1:470295/230‑1 (MQ=255)
                                                                                       cGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAAGCAGGAAAGGGCGGATGACCGGGCCGCCCATGGTACGAAAGGGAGTATTGCCCGATGGTTGTTTCGAGGCTACACCAAAACCGtttcaaccccagggg                                                                                                                                                           >  1:146530/1‑227 (MQ=255)
                                                                                         aGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAACCCTGTTTACCCTCCATGGTc                                                                                                                                                          >  2:322604/1‑238 (MQ=255)
                                                                                         aGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGGACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTcg                                                                                                                                                         >  2:428502/1‑239 (MQ=255)
                                                                                          gCCAGGTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTcg                                                                                                                                                         <  2:179860/238‑1 (MQ=255)
                                                                                              gCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGGAAAGGCGCAAGGCAGGAAAAGGAGGGTAACCCTGGCGCCCGTGGTCCGCAAGGATGTCTCGCCCAATAGTGGATTCGAGGCGACACCCAAACGCTTAACCCCCGAGGGGcgggggg                                                                                                                                                    >  1:605181/1‑234 (MQ=255)
                                                                                               gttgcggcacacgttACCGGTAGGGGCTCGGCGCTTGCGGGTACAGTGCGGACCCAGCGGGGTGAAATCTACCTGAACGGTGATGCCGACGGGAGCCAAAGTGAAAACGGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGGGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACGCCAAACCTGTTTACCCTCGATGGTCGCGTGGa                                                                                                                                                 <  1:538354/226‑1 (MQ=255)
                                                                                                ttGCGGTCAACTGTAACGGTATGCTCCCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGa                                                                                                                                                 <  2:453220/238‑1 (MQ=255)
                                                                                                tggCGGTCAACTTTAACGGTATGCCCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGa                                                                                                                                                 <  2:196247/236‑1 (MQ=255)
                                                                                                  gCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGAGCACCGGGCCGCCGGTGGTACGAATGGATGTATCGCCAGATGTTGGATTCGAGGCTACACCAAACCCGTTTACCCTCGTTGGTCGCGTGGATg                                                                                                                                               >  2:761240/1‑238 (MQ=255)
                                                                                                           tttAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGTCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGgc                                                                                                                                     <  2:603059/239‑1 (MQ=255)
                                                                                                           tttAACGGTAGGCGCTCGACGCTTGCCGGGACAGGACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACGGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGgc                                                                                                                                     <  1:404145/239‑1 (MQ=255)
                                                                                                                cGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAATGCGGAACACCCTGCCGCCCGTGGGACGAAATGATGATTCGCCCGAGGTGGTATTCGCGGCTACACCAAAACTGATTAACTCCGTGGGGCGCGTGGGATTTCCTTGGGgcggcc                                                                                                                                 >  1:215863/1‑236 (MQ=255)
                                                                                                                    aTGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCg                                                                                                                            <  2:515984/238‑1 (MQ=255)
                                                                                                                    aTGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGGACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCCCCATGGTCGCGTGGGTGTCCCTTGGGGGCGCCTCgt                                                                                                                           >  2:540032/1‑239 (MQ=255)
                                                                                                                     tGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGGGCGGGTCACCCGGCCGCCGAGGGTACGAAAGGATGTATCGCCAAATTTTGTATTTCAGGGCTAACCCAAACTGTTTTACCTCCGAGGTCGCTGGGGTGTCCCCTggggcccgccaggg                                                                                                                           >  1:515984/1‑228 (MQ=255)
                                                                                                                          tCGACGCTTGCCGGGACAGTACGGACCCAGCAGGGTGAACTCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCa                                                                                                                     <  2:82886/239‑1 (MQ=255)
                                                                                                                            gcccgttGCCGGTACAGGACGGACCCGGCAGGGTGAAATCTACCTGAACGGTGATGGCGACTGGAGGCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAAGGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACg                                                                                                                   <  2:681904/234‑1 (MQ=255)
                                                                                                                             aCGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCCGATGTTGTATTCGAGGCTACACCAAACCCGTTTACCCTCGATGGTCTCGTGGATGTCCCTTTGGGGCGCATCGTGGTGCACGa                                                                                                                  >  2:49758/1‑239 (MQ=255)
                                                                                                                             aCGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGGACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGa                                                                                                                  >  2:871576/1‑239 (MQ=255)
                                                                                                                                         aCAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGGGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAg                                                                                                      <  2:104004/239‑1 (MQ=255)
                                                                                                                                                ggACCCCGTGGGGTTAAATCTACCTGAACGGGGATGCCGATGGGGGCCAAATTGAAAACTGGGGGGCGGGGGTAACGGGGAAAGGGAGTAAAGGGCGGATCGCCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTa                                                                                               <  1:174464/239‑1 (MQ=255)
                                                                                                                                                  aCCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGt                                                                                                                                                                                                                              <  1:210564/116‑1 (MQ=255)
                                                                                                                                                  aCCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGt                                                                                                                                                                                                                              >  2:210564/1‑116 (MQ=255)
                                                                                                                                                  aCCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAgg                                                                                             >  2:708490/1‑239 (MQ=255)
                                                                                                                                                  aCCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGGGACCGATCTGCCGGAAAGCCCCGTAg                                                                                              >  2:710205/1‑238 (MQ=255)
                                                                                                                                                      agcagGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGGGCGGATCACCGTGCCGCCGATGGGACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGGGCACCGTGGGGCACGATCTCCCGGGAAGCGCAGATGGCgc                                                                                          >  2:592720/1‑237 (MQ=255)
                                                                                                                                                       gcagGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGGGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACATGTTTACCCTCGCTGGTCGCGTGGAAGTCCCCTGGGCGCGCATCGGGGGGCCCAATCTGCCGGAAAggcgcaaagggggc                                                                                         >  1:170233/1‑225 (MQ=255)
                                                                                                                                                       gcagGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGGAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGGCTGGAGGTCCCTTGGGCGCGCATCGTGGTGCCCGAATTGCCGGAAAGCGCAGTAGGCGtc                                                                                         >  2:455617/1‑238 (MQ=255)
                                                                                                                                                          gggTAAAATCTACCCGAACGGTGATGCCGGCTGGAACCAAATTGAAAACTGGGGGGCGCGGGTAACTGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTcc                                                                                      <  2:276817/238‑1 (MQ=255)
                                                                                                                                                           ggTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGa                                                                                                                                                 <  2:852457/179‑1 (MQ=255)
                                                                                                                                                              gaaagctgccgggggggTGAGGGCGTCGGGGGCCAAATTGAAAACTGGCGGGCGGGGAGAAGGGCGAAAGGCAATAAAGTGCGGATCACCGTGCCGCCGATGGATCGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGCCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCg                                                                                  <  1:322604/223‑1 (MQ=255)
                                                                                                                                                              gAAGTCTACCGGAACGGGGATGCCGACGGGAGACAAATTGAAAACTGGCGGGGGCGGGTAACGGCGAAAGGCAGTAAAGGGCGGATCACCGGGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGa                                                                                 <  2:441036/239‑1 (MQ=255)
                                                                                                                                                                aaTCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGAtg                                                                               <  2:170233/239‑1 (MQ=255)
                                                                                                                                                                     aCCTGAACGGTGATGCCGACTGGAGCCAAATTGAAACCTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATgtggtg                                                                          <  2:374012/239‑1 (MQ=255)
                                                                                                                                                                     aCCTGAAAGGTGAGGCCGAAGGGGGCCAAATTGAAAACGGGCGGGCCGGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCCCCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATgtggtg                                                                          <  2:324653/239‑1 (MQ=255)
                                                                                                                                                                          aaCGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAATGCGGATCACCGTGCCGCCGATGGTACGAAAGGATGGATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCCCCGTGTGGCGGGGGAGTTCCCTTTGGCGCCCATCTGGGGGGCAGATCTCCCGGAAAGCGCGATAGGGGCCTCCaccacaggggggatatt                                                                     >  1:270427/1‑223 (MQ=255)
                                                                                                                                                                                      ggatggtggcAAAGTGAAAACGGGGAGGGGGGGGTAACGGGGAAAAGGAGTAAAGGGCGGCTCCCGGTGCCGACGGTGGTACGAATGGAATGAACGGCAGATGTTGGATTCGAGGCTACACCAAACCTGTTTACCCGCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCTATGTGGTGATGCTTAACGATAACCt                                                         <  1:733073/230‑1 (MQ=255)
                                                                                                                                                                                        agggAGCAAAATTGAAAACTGGCGGGCGCGGGTAACGGCAAAAGGCAGTAAAGTGCGGATCCCCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTg                                                        <  2:605181/236‑1 (MQ=255)
                                                                                                                                                                                         ttgggCCAAATTGAAAATGGGCGGGGGCGGGTAACGGCGAAAGGGAGTAAAGGGCGGATCACCGGGCCGCCGATGGTACGAATGGATGGTTCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGc                                                       <  1:761240/234‑1 (MQ=255)
                                                                                                                                                                                         tGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGAAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGCGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCCGCGCTGTGGTGATGCTTAACGGTAAACTGCa                                                      >  2:303416/1‑239 (MQ=255)
                                                                                                                                                                                         gggtgCCAATTTAAAAACGGGCGGGGGCGGGTAAGGGCGAAAGGCGGTAAGGTGCGGATCCCCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCa                                                      <  1:708490/235‑1 (MQ=255)
                                                                                                                                                                                         ggaaGCCAAATTGAAAACTGGCGGGAGGGGGAAACGGCGGAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGGACGAATGGATGTAACGCCAGATGTTGTATTCGAGGCTTCACCAAACCTGTTTAACCTCCATGGTCGCGGGGATGTGCCTTTGGCGCGCACCGTGGGGCCCGAATTGCCGGGAAACCCAGTGAGGGACTCCCCCGAAGCGGGAATGCTTAAAGAAAACCCGcc                                                      >  1:39761/4‑238 (MQ=255)
                                                                                                                                                                                                 aaTTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGGAAAGGGGGGATCAACGTGCCGCCGAGAGTACGAATGGATGTATCGGCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGGTGGTAGCGTGGATGTTCCCTGGGAGCGCATAGTGGTGCCCGATCTGCCCGATAACGCAGTAGGCGTCCCCACGGACGTGGTGATACTTAACGAAAACCAGCAACCGaaaa                                              >  1:352586/1‑238 (MQ=255)
                                                                                                                                                                                                             gcgggcgCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGCGGATGTCCCTTGGGCGCGCATCGTGGGGGACGGTCTGCCGGAAAGGGCAGTAGGCGTCTCCCGCGATGTGGGGGTGCTTAACGAGAACCCCCAACCGGGAGAAGCGAAAAAgg                                  >  1:485532/1‑239 (MQ=255)
                                                                                                                                                                                                                          aCGGCGAAAGGCAGTAAAGTGCGGATCCCCGGGCCGCCGAGGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCGCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGAt                     <  2:709207/239‑1 (MQ=255)
                                                                                                                                                                                                                               aaaggggtttAAAGGGGGGTTCACCGGGCCGCCGATGGTACGAATGGATTTATCGCCAGGTGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAAc                 <  1:452702/229‑1 (MQ=255)
                                                                                                                                                                                                                                 aaGGCATTAAAGGGGGGTTCCCCGTGCCGCTGATGGTACAAATGGATGTATCCGCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAg               <  1:144187/238‑1 (MQ=255)
                                                                                                                                                                                                                                        tAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGGGATGCTTAAAGAGAACCTGCAACCGGAAGAGCCGAAAAAGGCGGGGATTCCGTTTAAAAGaaaacgga       >  1:743183/1‑234 (MQ=255)
                                                                                                                                                                                                                                         aaaGTGCGGATCCCCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCGCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGa       <  2:39761/238‑1 (MQ=255)
                                                                                                                                                                                                                                             tGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGt   >  2:758274/1‑238 (MQ=255)
                                                                                                                                                                                                                                             tGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGGGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGCCGTTTCCGATTAAAAGAAACCTGATTgac  >  1:33437/1‑237 (MQ=255)
                                                                                                                                                                                                                                             tGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCCCCCGCGATGTGGTGTTGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCTTCAATTCCGATTAAAAGAAAACTAATTGt   >  1:416274/1‑238 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
GTGATGTGCAAAGCCCGCAGTTGTTTGGTCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTC  >  NC_000913/4446733‑4447208

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