Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 221,963 A→C 12.7% V228G (GTG→GGG)  metN ← DL‑methionine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913221,9630AC12.7% 229.3 / 2.8 96V228G (GTG→GGG) metNDL‑methionine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (45/34);  new base C (2/13);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 3.63e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCCACCGGCGTAATCCATCTGCGCGCTAATAATGTTGTTGTTGACGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTTCATCACACAGCAATACTTTGGGATTGCTGGCTAACGCA  >  NC_000913/221725‑222187
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                
aaCGCCACCGGCGTAATCCATCTGCGCGCTAATAATGTTGTTGTTGACGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAaca                                                                                                                                                                                                                                  <  2:834863/239‑1 (MQ=255)
                  cATCTGCGCGCTAATAATGTTGTTGTTGACGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTcc                                                                                                                                                                                                                 >  2:24618/1‑238 (MQ=255)
                    tCTGCGCGTTATTAATGTTGGTGTTGACGGTGAAACGACGCGGGTTTTCAGAAAGCAGGGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGc                                                                                                                                                                                                              <  1:131887/239‑1 (MQ=255)
                        cgcgcTAATAATGTTGTTGTTGACGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCCGTAAATGGCTCCGCTTGCAGACGTTCCTGGCAATCTTCCGGGATATCCAGATGCAGGGGCGACTGAATAAACTTTCGC‑GCCA‑GCGGCGTTTTCGGATTCGTGACCACTTCACTTACCGTGTCCTGcgct                                                                                                                                                                                                           >  1:100825/1‑235 (MQ=255)
                        cgcgcTAATAATGTTGTTGTTGACGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGCGAACACTTCACTTACCGTGTc                                                                                                                                                                                                                  >  2:124507/1‑231 (MQ=255)
                         gcgcTAATAATTTTGTTGTTGACGTTGAAACGAGCGGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGAt                                                                                                                                                                                                         <  1:683456/239‑1 (MQ=255)
                         ccggctaatatttttttgcttgcGTTGAACGCACGCGCGGTTTCAAAAAGCAGTGGGGCATCGCCCGATTGACCGGTAAACCCCAGGCGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGAt                                                                                                                                                                                                         <  1:912647/217‑1 (MQ=255)
                         cccttatatatttgttttgtgcggTTAAAACGCGGCGCGGGTTCAGAAAGCAGGGTGGCATCGACAGATTGACCGGTAAACCCAAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGAt                                                                                                                                                                                                         <  1:340613/216‑1 (MQ=255)
                          cgcTAATAATGTTGTTGTTGACGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGACCACTTCACTTACCGTGTCCGGCTCGATc                                                                                                                                                                                                        >  2:399136/1‑239 (MQ=255)
                          cgcTAATAATGTTGTTGTTGACGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATc                                                                                                                                                                                                        >  2:446515/1‑239 (MQ=255)
                          cgataataaTGTTGTTTTGGGCGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGAt                                                                                                                                                                                                         <  2:933428/235‑1 (MQ=255)
                             taataaTGTTGTTGTGGACGTTGAAACGACGCGCGGTTTCAGAAAGCCGTGGGGCATCGACCGATTGACCGGTAAACTCCAGGCGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAg                                                                                                                                                                                                      <  2:1024762/238‑1 (MQ=255)
                                   ttttgttgtggACGTTGCAACGACGCGGGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGGCGGGTAAACTCCAGACGCAGGATCGGCACGCAGTCAGGAAATGGCTCCGCTTGCAGACGGTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCa                                                                                                                                                                                               <  1:1021001/237‑1 (MQ=255)
                                       gttgttgACGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCCGATGCAGGGCCGCCTGAATAAACTTCTCC‑GCCC‑GCGGCGTTTTCGGATGCGAGAACAATTCCCTTACCGGTCCCTGCCGGATCAGTTCCCCATTGc                                                                                                                                                                                           >  1:694933/1‑239 (MQ=255)
                                       gttgttgACGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCTAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCCtt                                                                                                                                                                                             >  2:56727/1‑237 (MQ=255)
                                       gttgttgACGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTc                                                                                                                                                                                                                  <  1:808798/216‑1 (MQ=255)
                                       gttgttgACGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTc                                                                                                                                                                                                                  >  2:808798/1‑216 (MQ=255)
                                       gttgttgACGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCg                                                                                                                                                                                                           <  2:195691/223‑1 (MQ=255)
                                       gttgttgACGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAAAAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCg                                                                                                                                                                                                           >  1:195691/1‑223 (MQ=255)
                                       gttgttgACGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACCGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCCCTTACCGTGTCCTGCCCGATCAGTTCTCCAtt                                                                                                                                                                                             >  1:179720/1‑237 (MQ=255)
                                             gACGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGCCGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGTAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGACAACCCTTCACTTACCGTGTCCTGTTCGATCAGTTTTCCATTgagtata                                                                                                                                                                                      >  2:596085/1‑232 (MQ=255)
                                               cGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGGGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATAAc                                                                                                                                                                                    >  2:389274/1‑238 (MQ=255)
                                                          cggggggTTTGAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAAt                                                                                                                                                                         <  2:927376/233‑1 (MQ=255)
                                                                 tttCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGAt                                                                                                                                                                                                         <  2:905563/199‑1 (MQ=255)
                                                                           cAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGGGAACACTTCACTTACCGTGTCCTGCTCGCTCAGTTCTCCATTGCTGATGACCCCCACGCAACCACAAATGCGCTTccaa                                                                                                                                                        >  1:453534/1‑235 (MQ=255)
                                                                           cAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGACACCTTCACTTACCGTTTCCTCCCCGATCATTTCCCCATTGCTGGGGACCGCCACGCAATCCAAAAttccgctccaca                                                                                                                                                       >  1:924690/1‑228 (MQ=255)
                                                                            aGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑ACCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACaa                                                                                                                                                       <  2:918628/238‑1 (MQ=255)
                                                                              tGGGGCATCGACCGGTTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCCCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGCGAACACTTCACTTACCGTGTCCTGCTCGCTCAGTTCTCCCTTCCTACTGACCGCCACGCAGTCCCACAAGCCCTTcaaacac                                                                                                                                                     >  1:669515/1‑233 (MQ=255)
                                                                              tGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCCGTTGACCGCCACGCAATCACAAATGCGCTTACCAACg                                                                                                                                                     >  2:774169/1‑238 (MQ=255)
                                                                               ggggCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTAGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCa                                                                                                                                                                                               >  2:396397/1‑195 (MQ=255)
                                                                                gggctgcgcccGATTGGCCGGTAAACGCCAGCCGGCGCAGCGGCACCCAGTCCGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTc                                                                                                                                                   <  1:24618/229‑1 (MQ=255)
                                                                                        aCCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGGGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGGCCTTCACAACGCCCATTTCggg                                                                                                                                          >  1:473734/1‑237 (MQ=255)
                                                                                         ccGGTTGACCGGTAACCCCCAGGCGCAGCCTCGGCAGGCAGTCCGTAAATGGCTCCGCTTGCAGACGTTCCTGGAAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTgg                                                                                                                                         <  2:69186/239‑1 (MQ=255)
                                                                                            aTTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACAGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTg                                                                                                                                      >  2:253871/1‑239 (MQ=255)
                                                                                             ttGCCCGGTAAACTCCAGACGAAGCATCGGCAGCCAGTCCGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTg                                                                                                                                      <  2:453534/238‑1 (MQ=255)
                                                                                             tggctcGGTAAGCTGCAGGCGACGCATCGTCACGCAGCCAGTAAAGTGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTGACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTg                                                                                                                                      <  2:200064/233‑1 (MQ=255)
                                                                                              tGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGAt                                                                                                                                    <  2:136160/239‑1 (MQ=255)
                                                                                                   ggTAAACTCCAGACGCAGACTCGGCACGGAGTCAGTAAATGGCTCCGTTTGCAGACGTTCCGGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACGTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACa                                                                                                                               <  2:100825/239‑1 (MQ=255)
                                                                                                     tAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGGGACCACTTCACTTACCGTGTCATGCTCGATCAGTTGTCCATTGATGATGACCGCCACGCAATCCCAAATTGCCTTCACAACGCCCTTTTCGGGGGGGATCACCaaa                                                                                                                             >  1:223263/1‑237 (MQ=255)
                                                                                                          tcccgtcgcttgctCGGCACGCAGTCAGAAATTGGATCTGGTGGCAGATGTTCCTGGTAATCTTACGGGATATCCAGTTGCAGGGTCGACTGAATAAACCTCTGC‑GCCA‑GCGGGCGTTTCGGATGCGAGAACCCTTCACATACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGCTGACCGATACGAAATAACAAAAGCGCTTCACAACGTACTTTTCGTTGGTGATAAACAGAATCGt                                                                                                                        <  1:183483/226‑1 (MQ=255)
                                                                                                          tCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGGTATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCGCGATCCGTTCTCCATTGCTGAGGACCGACACGCAATCAAAAACACGCTTGACAACGCCCCTTTCggtggcagtaacaagaaacc                                                                                                                         >  1:414239/1‑219 (MQ=255)
                                                                                                               aCGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACACCGTCCATTTCGTGGGTGATCAACAGAATCGTCAAc                                                                                                                    >  2:724951/1‑238 (MQ=255)
                                                                                                                    gcaTCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTCAGTGGGTGATCAACAGAATCGCAAACCCCaca                                                                                                              >  1:193050/1‑237 (MQ=255)
                                                                                                                         ggCACCGAGTCAGTAAATGGCCCCGGTTGCAGACCTTCCTGGTAATCTTCCGGGATAACCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGAcggcg                                                                                                         <  2:906197/239‑1 (MQ=255)
                                                                                                                              gcggtcggtAACTGGTTCCGTTTTCAAGCGTTGCTGGTACTCCTCGGGGTTGACCAGATGCGGGGGCGGCTTAATAAACTGTTGC‑GCCA‑GCGGCGTTTTCGGAGGCGAGAACACTTCCCTTACCGTTTCCTGCTCGGTCAGTTCTCCATTGCTGATGACCGCCACGCAATCACGAATGCGCTTCACAACGTCGATTTCGTGGGGGATCAAGAGAATCGTCAACCCCAGACGGCGGTTGa                                                                                                    <  1:317562/236‑1 (MQ=255)
                                                                                                                              gCAGTCAGGAAATGGCCCCGTTGGCGGACGTTCCGGGTAGACTCCCGGGATATCCATAAGCGGGGTGGCCTGAATAAACTTCTGC‑GCCG‑TGGGCGTTTTCGGAAGCGAGAACGCTTCCCTTGCCGGGTCCTGCTCGATCAGTTCTCCAGTGGTGATGGCCGCCACACAATCACAAATGCGCGTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGa                                                                                                    <  1:492403/239‑1 (MQ=255)
                                                                                                                               cagtcagtAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCCCAAATGCGCTTCACAACGTCCATTTCGTGGGGGATCACCAGAATCGTCAACCCCAGACGGCGGTTGag                                                                                                   >  2:689529/1‑238 (MQ=255)
                                                                                                                                      caatgggCCCCGTTTGCGGGCGTTCCTTGTAATCGTCCGGGATACCCAGATGCAGGGTCGGCTGAATAAACTTCTGC‑GCCC‑GCGGGTTTTTCGGATGCGGGAACCCTTCAATTACCGTTTCCTGATCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAGCGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTc                                                                                            <  1:446515/234‑1 (MQ=255)
                                                                                                                                         attgCTCCGCTTGCAGACGTTCCTTGTAATCTTCCGGGATAGCCAGAGTCAGGGCCGACCAAATAAACTTGTGC‑GCCA‑GCGGGCTTTTTGGTTGCGAGAACCTTTCACTTACCGTTTCCTGCTCGATCAGTTCTCCATTGCTGATGACGGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTcagc                                                                                         <  1:615208/236‑1 (MQ=255)
                                                                                                                                          tGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCCATCACAAATTCGCTTCACACAGTCCATTTCGTGGGTGATCAACAGAATCGTCCACCCCAGACGGCGGTTGATGTCTTTcagag                                                                                        >  1:297261/1‑237 (MQ=255)
                                                                                                                                          tGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACCGAATCGTAAACCCCAAACGGCGGTTGTTGTCTTTcggcc                                                                                        >  2:317562/1‑238 (MQ=255)
                                                                                                                                            gCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCACCAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGt                                                                                      >  2:979244/1‑239 (MQ=255)
                                                                                                                                              tCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCa                                                                                                                                                                                               <  2:901641/132‑1 (MQ=255)
                                                                                                                                              tCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCa                                                                                                                                                                                               >  1:901641/1‑132 (MQ=255)
                                                                                                                                              tCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGCCGGCGGTTGATGTCTTTCAGCCGTTc                                                                                    >  2:614347/1‑239 (MQ=255)
                                                                                                                                              tCCGCTTGCAGACGTTCCTGGTAATCCGCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTTC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTc                                                                                    <  2:317361/239‑1 (MQ=255)
                                                                                                                                                cGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCCCACATGCGCTTACCACCGTCCATTTCTGGGGTGTCCAACAGAATGGTAAACCCCAGACGGCGGTGGATGGCTTT‑AGGAGATCgag                                                                                 >  1:735690/1‑239 (MQ=255)
                                                                                                                                                cGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCga                                                                                  >  2:353526/1‑239 (MQ=255)
                                                                                                                                                  cTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCCCAAATGCGCTTCCCAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGAGTACTTTCAGCAGTTCcaga                                                                                >  2:249735/1‑239 (MQ=255)
                                                                                                                                                      cAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGGGGGTGATCACCCGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCCGTTCGAGAATAg                                                                            >  2:611292/1‑239 (MQ=255)
                                                                                                                                                          cGTTCCGTGTTATCTTCCGGGAATCCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCCCTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAAc                                                                         <  1:979244/238‑1 (MQ=255)
                                                                                                                                                             tCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACaa                                                                                                                                                       >  1:678190/1‑157 (MQ=255)
                                                                                                                                                             tCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACaa                                                                                                                                                       <  2:678190/157‑1 (MQ=255)
                                                                                                                                                                     aTCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGGGATCAACAGAATCGTCCACCCCAGACGGCGGGTGATGTCTTTCACAAGTTCGAGAACAAAACgtggagtgcgcg                                                             >  1:352907/1‑230 (MQ=255)
                                                                                                                                                                       cTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGCCGGCGGTTGATGTCTTTCAGCAGTTCGAGAATGAAACAGGTCGTTGCCggg                                                           >  1:672966/1‑239 (MQ=255)
                                                                                                                                                                       cTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGGGATCAACAGAATCGTCAACCCCAACCGGCGGTTGGTGTCTTTCAGCCGTTCGGGAATAGAACGTGGCGTTGCgggg                                                           >  1:935359/1‑239 (MQ=255)
                                                                                                                                                                         tCCGGGTTATCCGGGGGGAGGGTGGACTTAATAAACTTCTGG‑CCCA‑GCGGGGTTTTGGGATGCGAAAACCCTTCCCTTACCGTGTCCTGCCCGACCAGTTCCCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTc                                                         <  1:614347/239‑1 (MQ=255)
                                                                                                                                                                                  tCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGGGAATAGAACGGGTCGTTGCCGGGTCCAgccc                                                   >  2:791764/1‑234 (MQ=255)
                                                                                                                                                                                   ccgggtGCGGGGGCGACTACATAACCTTCTGC‑CCCA‑TGCGCGTTTACGGATGCGAGAACACTTCACTTCCCGGGTCCTGCTCGATCATGTCTCCATTGCTGATGACCGGCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCtg                                               <  2:428176/234‑1 (MQ=255)
                                                                                                                                                                                   ccAGATGCAGGGTCGACTTAATAAACTTCTGC‑GCCA‑GCGGCTTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCGGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCtg                                               <  2:352907/239‑1 (MQ=255)
                                                                                                                                                                                       aTGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGc                                                                                                                                                                                                              >  1:252880/1‑76 (MQ=255)
                                                                                                                                                                                       aTGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGc                                                                                                                                                                                                              <  2:252880/76‑1 (MQ=255)
                                                                                                                                                                                        tGCAGGTGCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACCCTTCCCTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCtggt                                             <  2:223263/236‑1 (MQ=255)
                                                                                                                                                                                            gggTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTGGTTGCCGGGTCCACCGCGCGTGTGGCTtc                                       >  2:700837/1‑238 (MQ=255)
                                                                                                                                                                                            gggTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTtc                                       >  2:437307/1‑238 (MQ=255)
                                                                                                                                                                                             ggTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGGGAATAGAACGTGTCGTTGCCGGGTCCAGCGCCCTGGGGGCTtca                                      >  2:7855/1‑238 (MQ=255)
                                                                                                                                                                                             ggTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGCGACTAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTtca                                      >  2:430609/1‑238 (MQ=255)
                                                                                                                                                                                               agccaggAATAAACTTCTGG‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTtcatca                                   <  2:832062/233‑1 (MQ=255)
                                                                                                                                                                                               acGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCCTTTCGTGCGGGATCAACTGAATCGTCAACCCCAGACGGCGGTTGATGTCTTGCAGAAGTTCGAAAATAGAACATGGCGTTGTCCGGCCCAgacccgtgggggctcccgac                                   >  1:214041/2‑220 (MQ=255)
                                                                                                                                                                                                tgccttaaTAAACTTCTGCGGCCG‑GCGG‑GTTTTCGGATGCGAGAACCCTTCACTTACCGTGTCCTGCTCGATCAGTTCCCCAGTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTtcatca                                   <  1:689529/232‑1 (MQ=255)
                                                                                                                                                                                                cgtcttaataatCTTCTTC‑GCCC‑GCGGGGTTTTCGGATGGGAAACCACTTGCCTTACCGTGCACTTGTCGACCAGTTCACCATGGCTGATGCCCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAAGCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTtcatca                                   <  1:399136/232‑1 (MQ=255)
                                                                                                                                                                                                cGACTTAATAAACTTCTGG‑CCCA‑GCGGCTTTTTCGGATGGGAGAACCCTTCACCTACCGTGTCCCGCTCGATCAGTTCCCCTTTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTtcatca                                   <  1:7855/238‑1 (MQ=255)
                                                                                                                                                                                                cGACTGATTTAACTTCTGT‑CCCA‑GCGGCGTTTTCGGGTGCGAAAACCCTTCCCTTCCCGGTGCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTtcatca                                   <  1:56727/238‑1 (MQ=255)
                                                                                                                                                                                                cGACTAAATAAACCTCTGC‑GCCA‑CCGGCGTTTTCGGTTGCGAGAACCCTTCACTTACCGTGCCCTGCCCGAACAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTtcatca                                   <  2:669515/238‑1 (MQ=255)
                                                                                                                                                                                                 gACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCCCTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTTCATcaca                                 <  2:907841/239‑1 (MQ=255)
                                                                                                                                                                                                   cttaaTAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCCCTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTTCATcacac                                <  2:694933/235‑1 (MQ=255)
                                                                                                                                                                                                      aaaaaaCTTCTGC‑CCCACGGGGCGTTTT‑GGGTGGCAGAACCCTTTCATTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTTCATCACACAGCa                            <  2:662899/236‑1 (MQ=255)
                                                                                                                                                                                                       agaaccgtCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGGGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACACCGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTTCATCACACAGCa                            <  2:600120/231‑1 (MQ=255)
                                                                                                                                                                                                       aTACACTTGTGC‑CCCA‑GCGGAGTTTTCGGATGCGAGAACCCTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCAGTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTTATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTTCATCACACAGCa                            <  2:851492/238‑1 (MQ=255)
                                                                                                                                                                                                       aTAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACCTTCCGTTGCCGGGGCAGGCCCCCTGGGGGCTtcaccaaaaagaaa                           >  1:396389/1‑231 (MQ=255)
                                                                                                                                                                                                        ataaagttcgtgcccA‑GGGGCGTTTTCGGATGCGGGAACCCTTCACTTACCGGGTCCTGCTCGATCAGTTCTCCATTGGTGATGACCGCCACGCAATCACAAATGCGCCTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTTCATCACACAGCaa                           <  2:28961/226‑1 (MQ=255)
                                                                                                                                                                                                                tGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGCCGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAAAACGGGTCGTTGCCGGGGCCGGCCCGCTGGGGGCTtcatcaaaaagaaaaaccggta                   >  1:832062/1‑222 (MQ=255)
                                                                                                                                                                                                                tGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTTCATACCACAGCAATACTTTggg                  >  2:309183/1‑239 (MQ=255)
                                                                                                                                                                                                                   ggccccgcggcgTTTTCGGAAGCGAGAACACTTCACTTCCCGTGCCCTGCGCGATAAGTTCTCCATTGGTGATGACCGCCCAGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTTCATCACACAGCAATACTTTGGGa                 <  1:791764/233‑1 (MQ=255)
                                                                                                                                                                                                                     ccA‑GCCGCTGTTGCGGATGCGAAAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTTCATCACACAGCAATACTTTGGGAtt               <  2:297261/238‑1 (MQ=255)
                                                                                                                                                                                                                             tgTTTTCGGATGCGAGAACCTTTCCCTTTCCTTGTCCGCCTCGATCAGTTCCCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTTCATCACACAGCAATACTTTGGGATTGCTGGCt        <  1:437307/237‑1 (MQ=255)
                                                                                                                                                                                                                                  tCGGATGCGAGAACTCTTCACTTACCGGGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTTCATCACACAGCAATACTTTGGGATTGCTGGCTAACGCa  <  2:179720/239‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                
AACGCCACCGGCGTAATCCATCTGCGCGCTAATAATGTTGTTGTTGACGTTGAAACGACGCGCGGTTTCAGAAAGCAGTGGGGCATCGACCGATTGACCGGTAAACTCCAGACGCAGCATCGGCACGCAGTCAGTAAATGGCTCCGCTTGCAGACGTTCCTGGTAATCTTCCGGGATATCCAGATGCAGGGTCGACTGAATAAACTTCTGC‑GCCA‑GCGGCGTTTTCGGATGCGAGAACACTTCACTTACCGTGTCCTGCTCGATCAGTTCTCCATTGCTGATGACCGCCACGCAATCACAAATGCGCTTCACAACGTCCATTTCGTGGGTGATCAACAGAATCGTCAACCCCAGACGGCGGTTGATGTCTTTCAGCAGTTCGAGAATAGAACGTGTCGTTGCCGGGTCCAGCGCGCTGGTGGCTTCATCACACAGCAATACTTTGGGATTGCTGGCTAACGCA  >  NC_000913/221725‑222187

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: ‑