Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,435,725 C→A 18.0% intergenic (‑106/+35) uspF ← / ← ompN stress‑induced protein, ATP‑binding protein/outer membrane pore protein N, non‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,435,7250CA18.0% 188.6 / 5.5 77intergenic (‑106/+35)uspF/ompNstress‑induced protein, ATP‑binding protein/outer membrane pore protein N, non‑specific
Reads supporting (aligned to +/- strand):  ref base C (57/5);  new base A (2/13);  total (59/18)
Fisher's exact test for biased strand distribution p-value = 4.21e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGT  >  NC_000913/1435488‑1435953
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                     
cAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCa                                                                                                                                                                                                                                      >  2:263386/1‑239 (MQ=255)
   tGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGaa                                                                                                                                                                                                                                   >  2:1112888/1‑239 (MQ=255)
     aaGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGaaaa                                                                                                                                                                                                                                 >  2:867975/1‑239 (MQ=255)
     aaGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAAAAAGCAGaaaa                                                                                                                                                                                                                                 >  1:148608/1‑239 (MQ=255)
           ttACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGccc                                                                                                                                                                                                                           >  2:703732/1‑238 (MQ=255)
             aCCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCCTAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGACCGTTGACACTTTTTGTGTAATGACGGGGTTCCCATTTTTAATTTAGATCAAAGGAGGAAGAAAAAACAGAAAAAGcccccc                                                                                                                                                                                                                        >  1:603374/1‑239 (MQ=255)
             aCCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTGTCCGAATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCCAAAGAGGAAGAATAAGCAGAAAAAGCCCGcc                                                                                                                                                                                                                        >  1:756970/1‑239 (MQ=255)
                 acAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACAATTTTTGTGTACTGACGGAGTTCCAGTTTTTAATTTAGAACAAAGGAGGAAGAATAAAAAGaaaaaaacaccaatac                                                                                                                                                                                                                     >  1:1144330/2‑228 (MQ=255)
                   aGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACa                                                                                                                                                                                                                  >  2:211675/1‑239 (MQ=255)
                                     ccaccatttttggCCTCTTCCTCATGGTGGTTAACCACGTGTGGAGTTAATTCTGAACCGAAAATATCGATAGGGAAAAGATTGGTTCTTTAAAAAAACCCTCTTTTTTTAGTACCCGAAAAAAGTGTACCGCAAGATGACTTTTTTGGGTAAATGGGGGAGTTCACATTTTTAATTAAAAAAAAAGGAGGAAGAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAg                                                                                                                                                                                                <  2:11755/227‑1 (MQ=255)
                                        tcaATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGa                                                                                                                                                                                               >  2:883427/1‑237 (MQ=255)
                                           aTCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAAAAGCGGGCAGGAGGATTTAGAAATGa                                                                                                                                                                                          >  1:11755/1‑239 (MQ=255)
                                                tGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGCATAAGCAGAAAAAACCCCACAAAAAAGCCGGCAGAGGGATTTaaaaagaaaaaaa                                                                                                                                                                                     >  1:912271/1‑230 (MQ=255)
                                                cgcttctctcTCAACGTGGCTATCCACGTGTTGATTTAATTTTGAATGGAAATTATGGAGAGGGAAAAAAATGGTTCTGTTAAAAAACCCTCTTTTTTTAGTCTCCGAAAAAGGTTCAACCCAAAAAGACATTTTTAGATTAAGGAGGGAGTCAACATTTTTAATTAAGAGAAAGGGAGAAAAAAAAACAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAc                                                                                                                                                                                     <  1:766441/230‑1 (MQ=255)
                                                            aaCGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGTAACTTTTTGTGTAATGCCGGGGTTCCCATTTTTAATTTAGATCAAAGGAGGAAAAATAAGCAGAAAAAGCCCGCCATAAACACGGGTGGGCGGATTTAGAACTTATAAACCACACCaaaaac                                                                                                                                                                         >  1:668631/1‑237 (MQ=255)
                                                                tGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACa                                                                                                                                                                     >  2:962956/1‑239 (MQ=255)
                                                                tGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAAAAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAACCGACa                                                                                                                                                                     >  2:981491/1‑239 (MQ=255)
                                                                  gCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAATAATCGTTCTGTTAATAAAACCTCCTGTTTTAGAATCAGCATAAAGTGGAACGCCAGATGACACTTTTTGTGTAATGCATGCGTTCACATTTTTAATTTAATTCCAAAGAGGAGAAATAAGAGAAAAAAGTCCGCCATAACAGAGGGCAGAAAGATGTAGAACTGATAAACCAGAACCaaaacacaaca                                                                                                                                                                   >  1:747165/1‑231 (MQ=255)
                                                                  gCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAAACCTAAAGCGACAat                                                                                                                                                                   >  2:1280502/1‑239 (MQ=255)
                                                                   caaaccACGCGTGGATTTAATCCAGAATCGGAAATACCGATAGGGACAAGAATCGTTCTGTTCATAAACCCTTGTTTTTTATTACCCGCATAAAGTGTAACGCAAGAGGACACTTTTGGTGTAATGAGGGATTTCACATTTTAAATTAAGACAAAAGGAGGAAGAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAat                                                                                                                                                                   <  2:527543/233‑1 (MQ=255)
                                                                    taatccggtgTTGATTTAATTCAGAATGGGAAATATCGATAGGGAAAAGAATCGTTTTTTTAAAAAACCCTCCTTTTTTAGTACCCGAAAAAGGTGTAACGCAAAATGACATTTTTTGGTAAAGGACGGAGTTCACATTTTTAATTTAGACAAAAGGAGGAAGAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAatat                                                                                                                                                                 <  2:1092016/230‑1 (MQ=255)
                                                                            cgTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAAACCCAAAGCGAAAAAATAATCGGTa                                                                                                                                                         >  1:1092016/1‑239 (MQ=255)
                                                                                           gAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCa                                                                                                                                          >  2:902395/1‑239 (MQ=255)
                                                                                            aaTCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAAATGATAAACCAAACCTAAAACGACAAAATAATCGGAAAAAATGCAAATGGAAg                                                                                                                                         >  1:67846/1‑239 (MQ=255)
                                                                                            aaTCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAAAAGCGGGCAGGAGGATTTAGAACTGATAAAACAGACCTAAAGCGACAAAAACATGGGTaactgtcccttaaaaaa                                                                                                                                         >  1:635705/1‑223 (MQ=255)
                                                                                                 gAAATATCGTTAGGGACAAGAATCGTTCTGTTCATAAAACCTACTGTTTTAGTATCCGCATAAAGCGTAACGCCAGATGTCACCTTTTTTGTAATGACGGTGTACACATTTTTAATTTAGATCAAAGTAGGAAGAATCAGCAGAAAAAGCCCGCATTAACAGCGGCAAGGAGGATTTAGAACTGCTAAACCAGACCTAAACCGCCAATACCAACGGTAAAGATTCAAATagcagcagaa                                                                                                                                    >  1:518394/1‑235 (MQ=255)
                                                                                                      atCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAAATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAACGACAAAACAATCGGGAAAGATTCCAAAGGCAACGAAAAAgag                                                                                                                               >  1:508539/1‑237 (MQ=255)
                                                                                                            aGGGACAAGACTCGTTCGGTTTATACCACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCACATGACCCTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCACAGGCGGAAGACTAAGCAGAAAAAGCCCGCCATAACAACGGGCAGGAGGATTTAGAACTGAAAAACCAGACCTAAAGCGACAATATCCTCGGTAGAAATGCCATAGGCACGGTAAAAGcaggaaaa                                                                                                                         >  1:673833/1‑232 (MQ=255)
                                                                                                            aGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGCCAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCa                                                                                                                           >  2:1238811/1‑237 (MQ=255)
                                                                                                                aCAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGCTGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAACGGAGGACGAATAAGCAGAAAAACCCCGCCATAACAGCGGGACGGAGGATTTAGAACATATAAAACAGACATAAAACGACAAAATCGACGGTAGAGAAGCAATTGCCaagcacaaactgggacacttca                                                                                                                     >  1:1302953/1‑218 (MQ=255)
                                                                                                                aCAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCa                                                                                                                     >  2:1228476/1‑239 (MQ=255)
                                                                                                                 cAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCACCTTAGAAGCTGTAACCTTct                                                                                                                     >  2:176085/1‑237 (MQ=255)
                                                                                                                         gTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAgg                                                                                                            >  2:1277367/1‑239 (MQ=255)
                                                                                                                          ttttGTTTATAAAACCTTTTTTTTTATTACCCAAAAAAAGGTAACGCCAAAGTGACATTTTGTGAGAAAGGAGGGACTTCCATTTATAAATTAAAAAAAAGGGAGAAAGAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGt                                                                                                           <  1:981491/236‑1 (MQ=255)
                                                                                                                            cTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACATCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCGGGCATCTTCATCCAACAGGTTg                                                                                                         >  2:359270/1‑239 (MQ=255)
                                                                                                                               ttCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAAAAAATGTAATCTTCATCAAACAGGGTGAtt                                                                                                      >  1:977890/1‑239 (MQ=255)
                                                                                                                                    aaccctttttttttttttATCCTCAAAAAGGGTACCGCAGGAGGATATTTTGTGATTAATGACGGATTCCTCTTTTTAAATATAGAACAAGGGAGGAATAAGAACAAAAAAAGCCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtata                                                                                                 <  1:1219598/222‑1 (MQ=255)
                                                                                                                                    aaaaCCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGAAGGAGGATTTAGAACTGACAAACAAGACCTAAAACGCCCATATCATCGGTAGAGGTGCCATTGGCAGCGTAGAAACTGGCACCTTAATCCAACAgggtactttataa                                                                                                 >  1:327259/1‑228 (MQ=255)
                                                                                                                                       ccctcatTTTTTAGTCCCCGTAAAAAGTGAAACCCAAAGTGACTTTTTTTGTGAAGTGGCGGATTCCACTTTTTAAATAAAAACAAAAGGAGAAAGAAAAAGCAAAAAAAGACCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTc                                                                                              <  2:259498/238‑1 (MQ=255)
                                                                                                                                         ctcctGTTTTAGTATCCGCATCAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGCAGAATAAGCAGAAAAAGCCCGCCCTAACACCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTTGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTAATCCAAAAGGTGGacatttaatcaaa                                                                                            >  1:259498/1‑227 (MQ=255)
                                                                                                                                         ctcctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCaa                                                                                            >  2:487451/1‑239 (MQ=255)
                                                                                                                                         ctcctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTACTCTTAATCCAAAAGGTTGATTTtaaaataaa                                                                                            >  1:36824/1‑232 (MQ=255)
                                                                                                                                           cctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGCCAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGCTTTTATAGTCAACg                                                                                          >  2:92408/1‑239 (MQ=255)
                                                                                                                                           cctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCAACGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATTTTCATCCAACAGGCTGACTTTATAGTCAACg                                                                                          >  2:663277/1‑239 (MQ=255)
                                                                                                                                                      aTCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGAGCAGGAGGATTTAGAACTGCTAAACCAGACCTAAAGAGACAATATCATCGTGAGAGATGCCCTTGGCAACGTAGAAGCGGTCACCTTCATCCAACAGGTTGATTTTATAGTGaaaaaagagggaaag                                                                               >  1:691616/1‑226 (MQ=255)
                                                                                                                                                         cGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGtt                                                                            >  2:93879/1‑239 (MQ=255)
                                                                                                                                                         cGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGtt                                                                            >  2:16498/1‑239 (MQ=255)
                                                                                                                                                         cGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGtt                                                                            >  2:203435/1‑239 (MQ=255)
                                                                                                                                                         cGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAACCGACACTATCATCGGTAGAGAAGCCATTTGCACCCTACAAGCTGTCAACTTCATCCAAAAGGTTGATTTTATAATCAACCAACGGGGACAtttt                                                                            >  1:376421/1‑239 (MQ=255)
                                                                                                                                                         cGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGttt                                                                           >  2:1200866/1‑239 (MQ=255)
                                                                                                                                                         cGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACAAGACCTAAAACGACAATATCATCGGTAGAGATGCCCATGGCAGCGTAGAAACTTGCATCTTCATCCAAAAGGTTGAATTTAAAGTCAACGAAGGGGGAACTGtt                                                                            >  1:996726/1‑239 (MQ=255)
                                                                                                                                                             aaaaGTTAAACGCCAGAGAACATTTTGTGTGAAAGGGCAGAGTTAACTTTTTTAATTAAGATAAAAGAAAAAAAAAAAAGAAGAAAAGGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTa                                                                        <  2:318592/238‑1 (MQ=255)
                                                                                                                                                             aaaaGGGAAACGACAGAGGACTCTTTTGGTGAAAGGACGGAGTTCAAATTAGTAATGAAGATAAAGGGAGAAAGAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTa                                                                        <  2:112920/238‑1 (MQ=255)
                                                                                                                                                                  gaaaaCCCCAGATGACTTTTTTGTGTAAATGGCGGAGTTACAATTTTAAATTAAAAACAAAGGAGGAAGAAAAAGAAAAAAAGGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTg                                                                     <  2:508539/234‑1 (MQ=255)
                                                                                                                                                                    tAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGAACTGTTTTTATTGAAAt                                                                 >  2:306574/1‑239 (MQ=255)
                                                                                                                                                                      aCGCCAGATGACAATTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCCGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGAAGAAGCTGTCATCTTCATCCACCAGGTTGCTTTTAAAGTCACCGTAGGTGGACAAGTTTTTTTTGAAATAg                                                               >  1:1067460/1‑239 (MQ=255)
                                                                                                                                                                           aGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTAAAGTCAACGTAGGGGGACATGTTTTTATTGAAATAGTAAg                                                           >  2:596488/1‑238 (MQ=255)
                                                                                                                                                                            gATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAATCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTCGAGATGCCATTGGCAACGTAGAAGCTGTCATCTTCCTCCAACACGTTGATTTATAAGTCAACGAAGGTGGACAGATTTTTATTCAAATAATAAggt                                                         >  1:453442/1‑237 (MQ=255)
                                                                                                                                                                            gATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGTGGGCAGGAGGATTTAGAACTGTTAAACCAGACCTAAACAGACAATATCATCGGTAGAGTTGCCATTGGCAGCGTCGAACCGGTAATCTCACTCCAAAAGGTTGATTTAATAGTAAACGAAGGTGGACAGGTTTTTATTTAAAAAGAAAGta                                                         >  1:858304/1‑238 (MQ=255)
                                                                                                                                                                                      ttttGTGTAAGAGCGGATTTCATATTTTTAATTTAGACAAAAGGAAGAAAATAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAata                                                <  2:912271/239‑1 (MQ=255)
                                                                                                                                                                                       tttGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAatat                                               >  2:1210015/1‑239 (MQ=255)
                                                                                                                                                                                            ttaaatagacgAGTTCACATTTATAATGTGATATAAAAGAGAGAAAAAAAACAAAAAAATCCCCGCCATAACTGCGGGCAGGAGGATTTAAAACTGAGAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCgg                                            <  1:388107/229‑1 (MQ=255)
                                                                                                                                                                                                atagtggATTTCATTTTTTTAATTAAATAAAAGGGAGAAAAAAAAAAAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTAt                                       <  1:204728/234‑1 (MQ=255)
                                                                                                                                                                                                aagACGGATTTCACTTTTTAAATTTAGATAAAAGGAGGAAAAAAAAGCAAAAAACGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTAt                                       <  2:36824/237‑1 (MQ=255)
                                                                                                                                                                                                      gAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAAcc                                 >  2:204728/1‑239 (MQ=255)
                                                                                                                                                                                                        gTTTACATTTTTAAATTTAGAAAAAAGGAGAAAAAAAAAGAAAAAAGAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAAcc                                 <  2:481592/238‑1 (MQ=255)
                                                                                                                                                                                                         ttCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGCTTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCACGCCAATATCGGCGTATTTAAcccaa                              >  2:698634/1‑236 (MQ=255)
                                                                                                                                                                                                         ttCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGa                              >  2:1203870/1‑239 (MQ=255)
                                                                                                                                                                                                                 tttAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAACTAGTAAGTTGCGCCAATATAGGCGAATTTAACCAGATCTTTatc                      >  2:1174844/1‑239 (MQ=255)
                                                                                                                                                                                                                 tttAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGAATTTAACCAGATCTTTgtc                      >  2:680864/1‑239 (MQ=255)
                                                                                                                                                                                                                      aTTTAGAAAAAAGGTGGAAAAAAAACAAAAAAACCCCCGTCAATAAAGCGGGCAGGAAGATTTAGAAGTGAAAAACCAGCCCTAAAGCGACAATATCATCGGTAGAGATGACATTGGCAGCGTAGAAGCTGTCATTTTCATGCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTAtcatca                  <  1:962956/239‑1 (MQ=255)
                                                                                                                                                                                                                        ttaaaaaaaaaGGGGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCaa                 <  2:67846/231‑1 (MQ=255)
                                                                                                                                                                                                                            gATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGAGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAAATTTTTATAATAAaaac              >  1:408133/1‑238 (MQ=255)
                                                                                                                                                                                                                                     gaggaAAAATAAGAAAAAAACGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCgg       <  2:977890/237‑1 (MQ=255)
                                                                                                                                                                                                                                        gaagaaTAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGTGCCAATATCGGCGAATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGt  >  2:1168158/1‑239 (MQ=255)
                                                                                                                                                                                                                                        gaagaaTAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTg   >  2:943356/1‑238 (MQ=255)
                                                                                                                                                                                                                                        gaagaaTAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTACAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTCTAGCCAACGTAGGCGGACATGTTTTTTTTGAAATAGTAAGTCGCGCCAACATCGGCGTATTTAACCAGATCTTTATCAAAAAAACCAACCGGGtggc  >  1:656546/1‑236 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                     
CAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGT  >  NC_000913/1435488‑1435953

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