Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,539,816 T→A 5.1% intergenic (+110/‑347) yhgF → / → feoA putative transcriptional accessory protein/ferrous iron transporter, protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,539,8160TA 5.1% 287.4 / 4.1 107intergenic (+110/‑347)yhgF/feoAputative transcriptional accessory protein/ferrous iron transporter, protein A
Reads supporting (aligned to +/- strand):  ref base T (57/44);  new base A (4/2);  total (61/46)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 110 total aligned reads displayed.

ATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGAG  >  NC_000913/3539668‑3539964
                                                                                                                                                    |                                                                                                                                                    
ctggatgCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTCCTGGCCTACGGTTCGAATATGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAat                                                                                                                                                      <  1:233679/148‑1 (MQ=255)
  ggatgCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAatat                                                                                                                                                    <  2:624268/149‑1 (MQ=255)
   gatgCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtc                                                                                                                                                   >  1:499002/1‑149 (MQ=255)
      gcgcTGGCGGCGGAAATTGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                                                                                                                                                <  1:113887/149‑1 (MQ=255)
        gcgggcggcggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAGATATTCTCCgg                                                                                                                                              <  2:223450/146‑1 (MQ=39)
        gcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGGAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGCACGCCTAGCCTGACATAAAAGTGCAGGAGAATATCTCCgg                                                                                                                                              >  1:453074/1‑149 (MQ=255)
               gcgggAATGGGATAAAACCGTGAAACGCCCGTACTGGCCAAAGGTTCTAATTTGCACGAAATCGTAGGGCAGATATGGCGTTCACGACTAAACTGGCAACGAACGCCTTTCCTGACATAAAAGTGGCGGAGAATATCTCCGGAAttttt                                                                                                                                       <  1:869784/149‑1 (MQ=255)
               gcggcAATGGGCAAAAAACGTTTAACGCCCGTACTGGCGTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGACGGAGAATATCTCCGGCAttttt                                                                                                                                       <  2:1173173/149‑1 (MQ=255)
                  gcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAtt                                                                                                                                    <  2:135766/149‑1 (MQ=255)
                  gcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAtt                                                                                                                                    >  2:390971/1‑149 (MQ=255)
                  gcAATGGGAAAACAACCTTAAACGCACGTACTGGCCGTCGGTTCGAATGTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGTACGCCTTTCCTTACATAAAAGTGCCGGAGAATAGCTCCGGCATTTTTAtt                                                                                                                                    <  1:273829/149‑1 (MQ=255)
                    aaTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGCACGCCTTGACTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTcc                                                                                                                                  >  2:655884/1‑149 (MQ=255)
                     aTGTGCAAAAAACGTTAAACGCCAGTACTGGCCTACGGGTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACACCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCca                                                                                                                                 <  2:152457/149‑1 (MQ=255)
                     aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTccc                                                                                                                                 >  2:65497/1‑148 (MQ=255)
                     aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAGATATTCTCCGGCACTTTTAtgtca                                                                                                                                 <  1:583976/149‑5 (MQ=38)
                      tGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTTGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGCACGCCTTGCCTGACATAAAAGTGCCGGAGCATATCTCCGGCATTTTTATTCcac                                                                                                                                >  2:175203/1‑149 (MQ=255)
                        ggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAg                                                                                                                              <  2:100560/149‑1 (MQ=255)
                         gCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGc                                                                                                                             >  1:1117977/1‑149 (MQ=255)
                          cAAAAAACGTTAAACGCCCGTACTGGCCTACGCTTCGAATTTGCACGAAATCGAAGGGCAGAAAAGGCGATCACGCCGCATCTGGCAACGAACGCCTTGCCTCACATAAAAGTGCCGGAGAATATCACCGGCATTTTTATTCCACCGcc                                                                                                                            >  2:437227/1‑149 (MQ=255)
                           aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCa                                                                                                                           >  1:862356/1‑149 (MQ=255)
                            aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCaa                                                                                                                          >  2:450599/1‑149 (MQ=255)
                              aaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAAc                                                                                                                        <  2:302120/149‑1 (MQ=255)
                              aaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAAACGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTCACTTAAAAGTGCCGGTGAATATCTCCGGCATTTTTATTCCACAGCCAAAc                                                                                                                        >  1:151368/1‑149 (MQ=255)
                               aaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACt                                                                                                                       <  2:968618/149‑1 (MQ=255)
                                aCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTc                                                                                                                      <  1:184862/149‑1 (MQ=255)
                                   ttAAAGGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCata                                                                                                                   >  1:575527/1‑149 (MQ=255)
                                   ttAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCata                                                                                                                   <  1:125128/149‑1 (MQ=255)
                                     aaaCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCataat                                                                                                                 >  2:575235/1‑149 (MQ=255)
                                       aCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAatat                                                                                                               >  1:451817/1‑149 (MQ=255)
                                       aCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAatat                                                                                                               >  1:459510/1‑149 (MQ=255)
                                           ccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATTAGGCGTTCACGCCGCATCTGGCAACGAAGGCCTTGCGTGACTTAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTc                                                                                                           <  1:39775/149‑1 (MQ=255)
                                           ccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCTGCATTTTTATTCCACAGCCAAACTCATAATATAtta                                                                                                           >  2:214994/1‑148 (MQ=255)
                                           ccGTACTGGCCTACGGTTCGAATTTGAACGAAAACTTAGTGCAGATAAGTCGTTCACGCCGCATCTGGCAAAGAACGACTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTc                                                                                                           <  1:1091610/149‑1 (MQ=255)
                                               aCTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGc                                                                                                       >  2:42612/1‑149 (MQ=255)
                                                 tGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTTTTCCACAGCCAAACTCATAATATATTCCGGCaa                                                                                                     >  1:799640/1‑149 (MQ=255)
                                                   gCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAata                                                                                                   >  1:617813/1‑149 (MQ=255)
                                                   gCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAata                                                                                                   >  1:232449/1‑149 (MQ=255)
                                                   gCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAata                                                                                                   >  2:541281/1‑149 (MQ=255)
                                                    ccTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAatat                                                                                                  >  2:869432/1‑149 (MQ=255)
                                                    ccTACGGTTCGAATTTGCACGAAATCCTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAatat                                                                                                  <  2:450363/149‑1 (MQ=255)
                                                     cTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATAtt                                                                                                 >  1:575371/1‑149 (MQ=255)
                                                      tACGGTTCGAATTTTCACTAAATCTTAGGTCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAGCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATAttt                                                                                                <  2:559096/149‑1 (MQ=255)
                                                      tACGGATCGAATTTGCTCGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGACCGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCAGAGCCAAACTCATAATATATTCCGGCAATAttt                                                                                                >  2:1069611/1‑149 (MQ=255)
                                                          gTTCGAATTAGCACGAAATCGGAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtca                                                                                            <  2:818462/149‑1 (MQ=255)
                                                          gTTCGAAATTGCACGAAATCGTAGGGCAGATAAGGCGATCACGACGCATCTGTCAACGAACGCCTTGCCTGACATTAAAGTGCCGTAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCAGGCAATATTTAtca                                                                                            <  2:1095272/149‑1 (MQ=255)
                                                           ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcat                                                                                           <  1:199381/149‑1 (MQ=255)
                                                           ttCGAATTTGCAAGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGCCAACCAACGCCTTAGCTGACATAAAAGTGCCGGAGAAAATATCCCTCAATTATATACCAGAGCCAAACTCAAAATATAAACCGGCAATATTTAacat                                                                                           >  2:366237/1‑149 (MQ=255)
                                                              gAATTTGCGCGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatttc                                                                                        >  2:524196/1‑149 (MQ=255)
                                                              gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCCCAGCCAAACTCATAATATATTCCGGCAATATTTAtcatttc                                                                                        >  2:766092/1‑149 (MQ=255)
                                                              gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatttc                                                                                        <  2:283122/149‑1 (MQ=255)
                                                              gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCAAAACTCATAATATATTCCGGCAATATTTAtcatttc                                                                                        >  1:110409/1‑149 (MQ=255)
                                                                      aCGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTaacaa                                                                                <  2:866633/149‑1 (MQ=255)
                                                                       cGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTaacaac                                                                               <  2:114320/149‑1 (MQ=255)
                                                                                aGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCtt                                                                      <  1:385809/149‑1 (MQ=255)
                                                                                 gggCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCtta                                                                     <  2:1077173/149‑1 (MQ=255)
                                                                                  ggCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaa                                                                    <  1:1125946/149‑1 (MQ=255)
                                                                                    cAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaatt                                                                  <  1:648583/149‑1 (MQ=255)
                                                                                       aTAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATTAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCCCAGCCAAACTCAGAATATATTTCGTCAATATTTATCATTTCATTACCAACTGAAACCGTATTTaaa                                                               >  1:231458/1‑149 (MQ=255)
                                                                                       aTAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa                                                               >  1:389519/1‑149 (MQ=255)
                                                                                       aTAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa                                                               >  1:963745/1‑149 (MQ=255)
                                                                                       aTAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa                                                               >  2:590943/1‑149 (MQ=255)
                                                                                       aTAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGACTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa                                                               <  2:796171/149‑1 (MQ=255)
                                                                                         aaGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCTAACTCATAAAATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAAGTAAACa                                                             >  1:1101215/1‑149 (MQ=255)
                                                                                             cGTTCACGCCGCATCTGGAAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAg                                                         <  2:247961/149‑1 (MQ=255)
                                                                                             cGTTCACCCCGCATCTCGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAATTTTATTCCACAGCCAAACACATAATATAATCAGGCTATATTAATCCTTTCATAAACAACAGAATCGTAAATTAAACaaaag                                                         >  2:828322/1‑145 (MQ=255)
                                                                                                  aagccgcATCTGGCAACTAACGCATTTCCTGACATAAAAGTGCCGTAGAATATCTCCTGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATAATTTCATTAAAAACTGAAACCTTAATTAAACATTAGACAGt                                                    >  2:402168/3‑149 (MQ=255)
                                                                                                   agccgcATCTGGCAAAGAACGCATTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCAAAGACATGCTCCTAATATATTCAGGCAATAATTAACATTTCATTAACAACTTAAACCTTAATTAAACATTAGCCAGTc                                                   >  2:877723/2‑149 (MQ=255)
                                                                                                     ccgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCg                                                 >  2:553996/1‑149 (MQ=255)
                                                                                                       gcATCTGGCAACGAACGCCTTGCCTGACATAAAAGGGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCggg                                               <  2:204840/149‑1 (MQ=255)
                                                                                                            tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGCATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAAtt                                          >  1:744553/1‑149 (MQ=255)
                                                                                                            tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAAtt                                          >  1:233267/1‑149 (MQ=255)
                                                                                                            tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCAGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAAtt                                          <  1:101936/149‑1 (MQ=255)
                                                                                                             ggCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTc                                         >  2:1004505/1‑149 (MQ=255)
                                                                                                              gCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCa                                        <  2:453074/149‑1 (MQ=255)
                                                                                                               cAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTATCAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCAc                                       >  1:217653/1‑149 (MQ=255)
                                                                                                                aacgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTACACAGCCAAACTGATAATATATTCCGGCAATATTTATAATTTCATTAACAACTGAAACCTTAATTAAACATTAGGCAGTACGGGAAATTCACt                                      >  1:450889/1‑149 (MQ=255)
                                                                                                                   gaacgCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTAtt                                   <  2:363368/149‑1 (MQ=255)
                                                                                                                     acgCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCg                                 >  2:539043/1‑149 (MQ=255)
                                                                                                                      cgCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGa                                <  1:655884/149‑1 (MQ=255)
                                                                                                                       gCCTTGCCTGACATAAAAGTTCCGGAGAATATCTCCGGCATTTTTTTTCCACAGCCAAACTCATAATATATTCGGGCTATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGaa                               <  2:1171550/149‑1 (MQ=255)
                                                                                                                       gCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGaa                               <  1:42612/149‑1 (MQ=255)
                                                                                                                        ccTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAAt                              <  2:1123782/149‑1 (MQ=255)
                                                                                                                        ccTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGGAATTCACTATTCGAAt                              >  1:273834/1‑149 (MQ=255)
                                                                                                                             ccTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTatat                         >  1:623001/1‑149 (MQ=255)
                                                                                                                             ccTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTatat                         >  1:695685/1‑149 (MQ=255)
                                                                                                                                gACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAATCATTAGCCAGTCCGGGTAATTCACTATTCGAATTATAtttt                      >  1:341856/1‑149 (MQ=255)
                                                                                                                                  cATAAAAGTGCCGGAGAATATCTCCGGCATTTTTTTTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCg                    <  2:799640/149‑1 (MQ=255)
                                                                                                                                  cATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCg                    <  2:862356/149‑1 (MQ=255)
                                                                                                                                   aTAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTACCAACTGAAACCTTAAGAAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGc                   >  2:1036272/1‑149 (MQ=255)
                                                                                                                                     aaaaGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTg                 >  2:493643/1‑149 (MQ=255)
                                                                                                                                       aaatGCCGGAGATATTCTCCGGCATTTTTATTCCACAGCCAATCTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCg               >  1:298618/4‑149 (MQ=255)
                                                                                                                                       aaatGCCGGAGATATTCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACAGTAGCCAGTCCGGGTAATTCACTATTAGAATTATATTTTCGCTGCg               >  1:614890/4‑149 (MQ=255)
                                                                                                                                        aatGCCGGAGATATTCTCCGGCACTTTTATTCCACAGCCAAACTCATAAGATATTCCGGCAATATATATCATTTCATTAACAACGGAAACCTTAATTAAACATTAGCCCGTCCGGGTAATTCACTATTCGAATTATATTTTAGCTGCGa              >  1:528283/3‑149 (MQ=255)
                                                                                                                                        aGTTCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAAGTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTTATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGa              <  2:713608/149‑1 (MQ=255)
                                                                                                                                          tGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGata            >  2:488079/1‑149 (MQ=255)
                                                                                                                                           gCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGatat           <  1:979392/149‑1 (MQ=255)
                                                                                                                                              ggAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAAc        <  1:65497/149‑1 (MQ=255)
                                                                                                                                              ggAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAAc        >  2:64730/1‑149 (MQ=255)
                                                                                                                                                agaATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCt      <  2:1034862/149‑1 (MQ=255)
                                                                                                                                                    tatCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGAg  <  2:959892/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGAG  >  NC_000913/3539668‑3539964

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