Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,539,816 T→A 7.3% 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 7.3% 267.1 / 8.0 81intergenic (+110/‑347)yhgF/feoAputative transcriptional accessory protein/ferrous iron transporter, protein A
Reads supporting (aligned to +/- strand):  ref base T (43/32);  new base A (3/3);  total (46/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGAG  >  NC_000913/3539668‑3539964
                                                                                                                                                    |                                                                                                                                                     
atggatgCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAat                                                                                                                                                       <  2:419409/149‑1 (MQ=255)
   gatgCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑tc                                                                                                                                                   >  1:191183/1‑149 (MQ=255)
   gatgCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑tc                                                                                                                                                   >  3:191688/1‑149 (MQ=255)
      gcgcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTcc                                                                                                                                                <  2:245208/149‑1 (MQ=255)
      gcgcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTcc                                                                                                                                                <  1:760289/149‑1 (MQ=255)
        gcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAGATAT‑TCTCCgg                                                                                                                                              <  1:81987/149‑1 (MQ=255)
            gcggcggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCAtt                                                                                                                                          <  4:585468/149‑1 (MQ=255)
            gcggcggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAG‑ATATTCTCCGGCact                                                                                                                                          <  3:595608/149‑3 (MQ=255)
                  gcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTAtt                                                                                                                                    <  1:115016/149‑1 (MQ=255)
                    aaTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTcc                                                                                                                                  >  2:458405/1‑149 (MQ=255)
                     aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAGATAT‑TCTCCGGCACTTTTAtgtca                                                                                                                                 <  2:812291/149‑5 (MQ=38)
                       gggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCcaca                                                                                                                               >  1:144598/1‑149 (MQ=255)
                         gCAAAAAACGTTAAACGCCCGTAATGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGTACGCCTTGCCTGACTTAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGc                                                                                                                             >  2:590467/1‑149 (MQ=255)
                           aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCa                                                                                                                           >  2:76473/1‑149 (MQ=255)
                            aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCaa                                                                                                                          >  4:866713/1‑149 (MQ=255)
                             aaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATGAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCaaa                                                                                                                         <  1:397455/149‑1 (MQ=255)
                             aaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGCGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCaaa                                                                                                                         <  4:731770/149‑1 (MQ=255)
                                aCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTc                                                                                                                      >  4:696853/1‑149 (MQ=255)
                                aCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTc                                                                                                                      >  2:147695/1‑149 (MQ=255)
                                     aaaCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCataat                                                                                                                 >  1:219878/1‑149 (MQ=255)
                                         gCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAatatat                                                                                                             >  3:165853/1‑149 (MQ=255)
                                           ccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTc                                                                                                           >  3:219336/1‑149 (MQ=255)
                                              tACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCgg                                                                                                        <  3:618933/149‑1 (MQ=255)
                                               aCTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGc                                                                                                       <  2:485887/149‑1 (MQ=255)
                                                 tGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCaa                                                                                                     <  1:891621/149‑1 (MQ=255)
                                                   gCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAata                                                                                                   >  1:522/1‑149 (MQ=255)
                                                   gCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAata                                                                                                   >  4:751490/1‑149 (MQ=255)
                                                        cGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAt                                                                                              <  4:445519/149‑1 (MQ=255)
                                                          gTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtca                                                                                            >  3:250469/1‑149 (MQ=255)
                                                           ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcat                                                                                           >  1:677218/1‑149 (MQ=255)
                                                           ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcat                                                                                           >  4:646230/1‑149 (MQ=255)
                                                               aaTTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatttca                                                                                       <  3:105591/149‑1 (MQ=255)
                                                                aTTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatttcat                                                                                      <  3:235859/149‑1 (MQ=255)
                                                                aTTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatttcat                                                                                      <  1:648642/149‑1 (MQ=255)
                                                                    gCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAGATAT‑TCTCCGGCac                                                                                                                                           <  4:328887/92‑2 (MQ=37)
                                                                    gCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAG‑ATATTCTCCGGCac                                                                                                                                           >  3:328887/1‑91 (MQ=37)
                                                                      aCGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTaacaa                                                                                >  4:562998/1‑149 (MQ=255)
                                                                        gAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACt                                                                              >  2:786891/1‑149 (MQ=255)
                                                                             cGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAAc                                                                         <  2:120963/149‑1 (MQ=255)
                                                                                 gggCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCtta                                                                     >  4:649944/1‑149 (MQ=255)
                                                                                  ggCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaa                                                                    <  4:177716/149‑1 (MQ=255)
                                                                                   gCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaat                                                                   >  1:458425/1‑149 (MQ=255)
                                                                                   gCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaat                                                                   <  3:633196/149‑1 (MQ=255)
                                                                                      gATGAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaattaa                                                                >  4:463355/1‑149 (MQ=255)
                                                                                      gATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaattaa                                                                <  4:77203/149‑1 (MQ=255)
                                                                                      gATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaattaa                                                                >  3:276465/1‑149 (MQ=255)
                                                                                       aTAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACGTTAATTaaa                                                               <  4:649710/149‑1 (MQ=255)
                                                                                       aTAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa                                                               <  4:191688/149‑1 (MQ=255)
                                                                                       aTAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCCGCAATATTTATCATTTCATTAACAACTGAAACCttaattaag                                                               <  3:696853/149‑2 (MQ=255)
                                                                                       aTAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACCCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa                                                               <  2:191183/149‑1 (MQ=255)
                                                                                        tAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAAc                                                              <  2:144598/149‑1 (MQ=255)
                                                                                            gCGTTCACGCCGCATCGGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTa                                                          >  4:600838/1‑149 (MQ=255)
                                                                                              gTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGc                                                        >  4:256639/1‑149 (MQ=255)
                                                                                                  aCGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGt                                                    >  4:314514/1‑149 (MQ=255)
                                                                                                  aCGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATCTATCATTTCATTAACAACTGAAACCTTAATTAAACATGAGCCAGt                                                    >  3:760147/1‑149 (MQ=255)
                                                                                                    gccgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTcc                                                  >  4:371136/1‑149 (MQ=255)
                                                                                                    gccgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTcc                                                  >  4:154670/1‑149 (MQ=255)
                                                                                                            tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAAtt                                          <  1:76473/149‑1 (MQ=255)
                                                                                                             ggCAACGAACGCCTTGCCTGACATAAAAATGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTc                                         >  3:86002/1‑149 (MQ=255)
                                                                                                                aacgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACt                                      >  3:137781/1‑149 (MQ=255)
                                                                                                                 acgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTa                                     <  3:226060/149‑1 (MQ=255)
                                                                                                                   gaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCACACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTAtt                                   >  4:841561/1‑149 (MQ=255)
                                                                                                                   gaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTAtt                                   >  3:824185/1‑149 (MQ=255)
                                                                                                                   gaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACCCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTAtt                                   <  1:262715/149‑1 (MQ=255)
                                                                                                                     acgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCg                                 >  3:341791/1‑149 (MQ=255)
                                                                                                                       gCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAGCTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGaa                               <  3:202085/149‑1 (MQ=255)
                                                                                                                       gCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGaa                               >  1:17104/1‑149 (MQ=255)
                                                                                                                       gCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGaa                               <  4:219336/149‑1 (MQ=255)
                                                                                                                       gCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGaa                               <  2:17104/149‑1 (MQ=255)
                                                                                                                       gCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCAGGTAATTCACTATTCGaa                               >  1:611196/1‑149 (MQ=255)
                                                                                                                          ttGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATta                            >  1:574536/1‑149 (MQ=255)
                                                                                                                               tGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAGCTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATAttt                       >  2:58775/1‑149 (MQ=255)
                                                                                                                                   aTAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGc                   >  3:406379/1‑149 (MQ=255)
                                                                                                                                   aTAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGc                   <  4:624939/149‑1 (MQ=255)
                                                                                                                                   aTAAAAATGCCGGAGATAT‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTCATTCACTATTCGAATTATATTTTCGc                   >  1:601065/1‑149 (MQ=255)
                                                                                                                                   aTAAAAATGCCGGAGATAT‑TCTCCGGCACTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGc                   >  2:93561/1‑149 (MQ=39)
                                                                                                                                      aaaGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGc                <  3:34902/149‑1 (MQ=255)
                                                                                                                                        aatGCCGGAGATAT‑TCTCCGGCACTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGa              >  1:400306/3‑149 (MQ=255)
                                                                                                                                             cGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATaa         <  3:463355/149‑1 (MQ=255)
                                                                                                                                              ggAGAATA‑TCTCCGTCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAAc        >  2:325174/1‑149 (MQ=255)
                                                                                                                                               gagaATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAAcc       >  3:375005/1‑149 (MQ=255)
                                                                                                                                               gagaATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAAcc       >  3:217058/1‑149 (MQ=255)
                                                                                                                                                   ata‑tCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGa   <  3:520683/149‑1 (MQ=255)
                                                                                                                                                   ata‑tCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGa   <  2:522/149‑1 (MQ=255)
                                                                                                                                                    ta‑tCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGAg  <  1:61242/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGAG  >  NC_000913/3539668‑3539964

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑