Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,539,805 G→A 5.4% 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 5.4% 362.5 / 7.2 110intergenic (+99/‑358)yhgF/feoAputative transcriptional accessory protein/ferrous iron transporter, protein A
Reads supporting (aligned to +/- strand):  ref base G (57/47);  new base A (4/2);  total (61/49)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

Only 100 of 110 total aligned reads displayed.

GCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGA  >  NC_000913/3539660‑3539952
                                                                                                                                                 |                                                                                                                                                   
gcgcGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGc                                                                                                                                                  >  1:677354/1‑149 (MQ=255)
gcgcGATGATGGATGCGATGGCGGCGGCAATGGGCAAAAAAGGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGGAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGTCATAAAAGTGc                                                                                                                                                  >  4:668369/1‑149 (MQ=255)
    gatgatGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGc                                                                                                                                                  >  1:1033814/1‑145 (MQ=255)
    gatgatGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGc                                                                                                                                                  <  2:1033814/145‑1 (MQ=255)
    gatgatGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGga                                                                                                                                              <  4:76737/149‑1 (MQ=255)
      tgatgGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                                                                                                                            <  2:1029572/149‑1 (MQ=255)
         tggatgCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAata                                                                                                                                         >  4:656841/1‑149 (MQ=255)
          ggatgCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAatat                                                                                                                                        >  1:201429/1‑149 (MQ=255)
          ggatgCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAatat                                                                                                                                        <  3:746437/149‑1 (MQ=255)
                   ggcggcggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAt                                                                                                                               <  1:135689/149‑1 (MQ=255)
                    gcggcggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtt                                                                                                                              <  3:1054401/149‑1 (MQ=255)
                       gcggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAat                                                                                                                                          >  3:71370/1‑134 (MQ=255)
                         ggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAt                                                                                                                         >  3:542891/1‑149 (MQ=255)
                          gcAATGGGCAAAAAACGTTAAACGGCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAGCGCCTTGCCTGACATAAAAGTGCCGCAGAATATCTCCGGGATTTTTAtt                                                                                                                        >  3:676471/1‑149 (MQ=255)
                          gcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAtt                                                                                                                        <  2:843802/149‑1 (MQ=255)
                            aaTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTcc                                                                                                                      >  3:1029873/1‑149 (MQ=255)
                               gggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcaca                                                                                                                   <  4:872653/149‑1 (MQ=255)
                                 gCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGc                                                                                                                 >  3:790750/1‑149 (MQ=255)
                                      aaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAAc                                                                                                            >  2:950859/1‑149 (MQ=255)
                                          gTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCat                                                                                                        <  4:173324/149‑1 (MQ=255)
                                          gTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCat                                                                                                        >  3:173324/1‑149 (MQ=255)
                                              aaCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGc                                                                                                                                                  <  1:795950/103‑1 (MQ=255)
                                              aaCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGc                                                                                                                                                  >  2:795950/1‑103 (MQ=255)
                                                 gCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAatatat                                                                                                 >  2:22708/1‑149 (MQ=255)
                                                   ccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTc                                                                                               <  1:108993/149‑1 (MQ=255)
                                                      tACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATTAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCgg                                                                                            <  4:1029873/149‑1 (MQ=255)
                                                        cTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCa                                                                                          >  1:89207/1‑149 (MQ=255)
                                                         tGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCaa                                                                                         >  2:338944/1‑149 (MQ=255)
                                                            ccTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTGTATTCCACAGCCAAACTCATAATATATTCCGGCAatat                                                                                      >  4:355309/1‑149 (MQ=255)
                                                              tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCg                                                                                             <  2:770632/140‑1 (MQ=255)
                                                              tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCg                                                                                             >  1:770632/1‑140 (MQ=255)
                                                               aCGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTa                                                                                   <  4:403232/149‑1 (MQ=255)
                                                                  gTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGGAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCataa                                                                                                      <  3:254051/127‑1 (MQ=255)
                                                                  gTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCataa                                                                                                      >  4:254051/1‑127 (MQ=255)
                                                                  gTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtca                                                                                >  1:152033/1‑149 (MQ=255)
                                                                   ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcat                                                                               >  1:687303/1‑149 (MQ=255)
                                                                   ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcat                                                                               >  3:227025/1‑149 (MQ=255)
                                                                   ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcat                                                                               >  4:814582/1‑149 (MQ=255)
                                                                   ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcat                                                                               >  4:757752/1‑149 (MQ=255)
                                                                   ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcat                                                                               >  2:221409/1‑149 (MQ=255)
                                                                   ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAGATATTCTCCGGCACTTTTAt                                                                                                                         >  4:248211/1‑107 (MQ=38)
                                                                   ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAGATATTCTCCGGCACTTTTAt                                                                                                                         <  3:248211/107‑1 (MQ=38)
                                                                      taaTTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAATTCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatttc                                                                            <  4:855387/148‑1 (MQ=255)
                                                                      gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatttc                                                                            >  4:27676/1‑149 (MQ=255)
                                                                      gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatttc                                                                            <  4:949405/149‑1 (MQ=255)
                                                                       aaTTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatttca                                                                           <  4:497987/149‑1 (MQ=255)
                                                                       aaTTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatttca                                                                           >  4:604082/1‑149 (MQ=255)
                                                                        aTTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatttcat                                                                          <  1:961035/149‑1 (MQ=255)
                                                                              aCGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTaacaa                                                                    >  1:261879/1‑149 (MQ=255)
                                                                              aCGAAATCGTAGGGCAGATAAGGCGTTCACGCAGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATCCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTaacac                                                                    <  2:1035042/149‑2 (MQ=255)
                                                                               cGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTaacaac                                                                   <  2:200447/149‑1 (MQ=255)
                                                                                gAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtt                                                                                                                              >  4:173703/1‑89 (MQ=255)
                                                                                gAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACt                                                                  >  4:312911/1‑149 (MQ=255)
                                                                                gAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACt                                                                  >  2:186727/1‑149 (MQ=255)
                                                                                    tCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGaaa                                                              <  1:1040811/149‑1 (MQ=255)
                                                                                    tCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGaaa                                                              <  1:295891/149‑1 (MQ=255)
                                                                                      gTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAAcc                                                            >  3:27111/1‑149 (MQ=255)
                                                                                        aGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCtt                                                          <  1:777564/149‑1 (MQ=255)
                                                                                          ggCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaa                                                        <  3:604082/149‑1 (MQ=255)
                                                                                          ggCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAATTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATGATGTCATTAACAACTGAAACCtgaa                                                        <  4:672919/149‑1 (MQ=255)
                                                                                           gCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaat                                                       >  1:124864/1‑149 (MQ=255)
                                                                                              gATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaattaa                                                    >  1:357889/1‑149 (MQ=255)
                                                                                              gATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaattaa                                                    <  2:261879/149‑1 (MQ=255)
                                                                                                    gCGTTCACGCCGCATCTGGCAGCGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTaacaa                                                                    >  2:31689/1‑127 (MQ=255)
                                                                                                    gCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTaacaa                                                                    <  1:31689/127‑1 (MQ=255)
                                                                                                    gCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTa                                              >  3:339643/1‑149 (MQ=255)
                                                                                                     cGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAg                                             <  1:427732/149‑1 (MQ=255)
                                                                                                      gTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGc                                            <  1:950859/149‑1 (MQ=255)
                                                                                                      gTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGc                                            <  3:20491/149‑1 (MQ=255)
                                                                                                          aCGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGt                                        >  3:971464/1‑149 (MQ=255)
                                                                                                              cgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCgg                                    <  4:944983/149‑1 (MQ=255)
                                                                                                                    tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAAc                                                             <  3:219247/118‑1 (MQ=255)
                                                                                                                    tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAAc                                                             >  4:219247/1‑118 (MQ=255)
                                                                                                                    tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAAtt                              <  1:526027/149‑1 (MQ=255)
                                                                                                                        tacgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGCCATTTATATTCCACAGCCAAACGTATAATATATTCCGGCACTATTTATCATTTCATTAACAACTGAAACCTTATTAAAACATTAGCTAGACCGGGTAATTCAcc                          >  4:609997/2‑148 (MQ=255)
                                                                                                                        aacgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACt                          >  4:82448/1‑149 (MQ=255)
                                                                                                                        aacgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACt                          >  4:814640/1‑149 (MQ=255)
                                                                                                                            aacgCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTc                      >  1:111707/1‑149 (MQ=255)
                                                                                                                            aacgCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTc                      <  3:757752/149‑1 (MQ=255)
                                                                                                                               gCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGaa                   <  4:279140/149‑1 (MQ=255)
                                                                                                                                ccTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAAt                  <  2:357889/149‑1 (MQ=255)
                                                                                                                                 cTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAAtt                 <  2:111707/149‑1 (MQ=255)
                                                                                                                                     ccTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTatat             >  4:694443/1‑149 (MQ=255)
                                                                                                                                       tGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATAttt           >  4:380870/1‑149 (MQ=255)
                                                                                                                                       tGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATAttt           >  1:414397/1‑149 (MQ=255)
                                                                                                                                       tGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATAttt           >  3:414579/1‑149 (MQ=255)
                                                                                                                                        gACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATAtttt          <  1:616103/149‑1 (MQ=255)
                                                                                                                                         gaaTAAAAATGCCGGAGATATTCTCCGGCACTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTc         >  4:318868/3‑149 (MQ=39)
                                                                                                                                         aCATAAAAGTGCCGGAGAATATCTCCTGCATTTTTATTCCACAGCCAAACTCATAACATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTACCCAGTCCGGGTAATTCACTATTCGAATTATATTTTc         <  2:112649/149‑1 (MQ=255)
                                                                                                                                         aCATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTc         >  3:804177/1‑149 (MQ=255)
                                                                                                                                          aaTAAAAATGCCGGAGATATTCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCg        >  1:204124/2‑149 (MQ=255)
                                                                                                                                           aTAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcaca                                                                                                                   <  3:930298/41‑1 (MQ=255)
                                                                                                                                           aTAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcaca                                                                                                                   >  4:930298/1‑41 (MQ=255)
                                                                                                                                           aTAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGc       >  1:362673/1‑149 (MQ=255)
                                                                                                                                           aTAAAAGTGCCGGAGAATATCTACGGCATTTTTATTCCACAGCCAAACTCACAATATATCCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATGAAACATTCGCCAGTCCGGGTAATTCACTATTCGAATTATATTATCGc       >  4:674151/1‑149 (MQ=255)
                                                                                                                                            tAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCAc                           <  2:691920/128‑1 (MQ=255)
                                                                                                                                            tAAAAATGCCGGAGATATTCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGt                                  >  1:797480/1‑121 (MQ=255)
                                                                                                                                            tAAAAATGCCGGAGATATTCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGt                                  <  2:797480/121‑1 (MQ=39)
                                                                                                                                               aaGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCg   >  2:524153/1‑149 (MQ=255)
                                                                                                                                                aGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGa  <  3:312911/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGA  >  NC_000913/3539660‑3539952

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: