Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,539,805 G→A 9.3% intergenic (+99/‑358) 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,8050GA 9.3% 250.0 / 12.0 85intergenic (+99/‑358)yhgF/feoAputative transcriptional accessory protein/ferrous iron transporter, protein A
Reads supporting (aligned to +/- strand):  ref base G (43/34);  new base A (4/4);  total (47/38)
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.78e-01

AGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGC  >  NC_000913/3539659‑3539950
                                                                                                                                                  |                                                                                                                                                  
aGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTg                                                                                                                                                  >  4:235859/1‑149 (MQ=255)
  cgcgATGATGGATGCGCTGGCGGCGGCAATGGGCAAACAACGTTAGACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGcc                                                                                                                                                <  1:721276/149‑1 (MQ=255)
  cgcgATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGcc                                                                                                                                                <  1:446205/149‑1 (MQ=255)
  cgcgATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTAGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGcc                                                                                                                                                >  4:787536/1‑149 (MQ=255)
     gatgatGGATGCGCTGGCGGCGGCAATGGGCAAAAAGCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGACATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGACTGACATAAAAGTGCCGga                                                                                                                                             <  2:570087/149‑1 (MQ=255)
     gatgatGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGga                                                                                                                                             >  1:689995/1‑149 (MQ=255)
     gatgatGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGga                                                                                                                                             >  1:382559/1‑149 (MQ=255)
        gatggatgCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGaa                                                                                                                                          >  2:687304/1‑149 (MQ=255)
         atggatgCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAat                                                                                                                                         <  2:419409/149‑1 (MQ=255)
            gatgCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑tc                                                                                                                                     >  3:191688/1‑149 (MQ=255)
            gatgCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑tc                                                                                                                                     >  1:191183/1‑149 (MQ=255)
               gcgcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTcc                                                                                                                                  <  1:760289/149‑1 (MQ=255)
               gcgcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTcc                                                                                                                                  <  2:245208/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                                                                        <  1:648642/149‑1 (MQ=255)
                                                                         aTTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatttcat                                                                        <  3:235859/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                                                  >  3:276465/1‑149 (MQ=255)
                                                                                               gATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaattaa                                                  <  4:77203/149‑1 (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:154670/1‑149 (MQ=255)
                                                                                                             gccgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTcc                                    >  4:371136/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                 <  4:219336/149‑1 (MQ=255)
                                                                                                                                gCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGaa                 >  1:17104/1‑149 (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     <  4:624939/149‑1 (MQ=255)
                                                                                                                                            aTAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGc     >  3:406379/1‑149 (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)
                                                                                                                                                  |                                                                                                                                                  
AGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGC  >  NC_000913/3539659‑3539950

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