Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,539,817 A→T 6.9% intergenic (+111/‑346) 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,8170AT 6.9% 340.7 / 12.0 101intergenic (+111/‑346)yhgF/feoAputative transcriptional accessory protein/ferrous iron transporter, protein A
Reads supporting (aligned to +/- strand):  ref base A (52/42);  new base T (5/2);  total (57/44)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

Only 100 of 106 total aligned reads displayed.

TGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGAG  >  NC_000913/3539669‑3539964
                                                                                                                                                    |                                                                                                                                                   
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)
                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)
                                        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                                                                                                         >  1:770632/1‑140 (MQ=255)
                                                     tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCg                                                                                                         <  2:770632/140‑1 (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                                                                                           >  3:227025/1‑149 (MQ=255)
                                                          ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcat                                                                                           >  2:221409/1‑149 (MQ=255)
                                                          ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcat                                                                                           >  1:687303/1‑149 (MQ=255)
                                                          ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcat                                                                                           >  4:160883/1‑149 (MQ=255)
                                                          ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcat                                                                                           >  4:757752/1‑149 (MQ=255)
                                                          ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcat                                                                                           >  4:814582/1‑149 (MQ=255)
                                                          ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAGATATTCTCCGGCACTTTTAt                                                                                                                                     <  3:248211/107‑1 (MQ=38)
                                                          ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAGATATTCTCCGGCACTTTTAt                                                                                                                                     >  4:248211/1‑107 (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: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)
                                                                       gAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAGAAGTGCCGGAGAATATCTCCGGCAtt                                                                                                                                          <  3:173703/89‑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:295891/149‑1 (MQ=255)
                                                                           tCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGaaa                                                                          <  1:1040811/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)
                                                                                  gCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaat                                                                   >  2:404276/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                                                          >  1:271082/1‑149 (MQ=255)
                                                                                           gCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTa                                                          >  3:339643/1‑149 (MQ=255)
                                                                                             gTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGc                                                        <  3:20491/149‑1 (MQ=255)
                                                                                             gTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGc                                                        <  1:950859/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                              <  4:364929/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                       >  1:414397/1‑149 (MQ=255)
                                                                                                                              tGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATAttt                       >  3:414579/1‑149 (MQ=255)
                                                                                                                              tGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATAttt                       >  4:380870/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)
                                                                                                                                 aaTAAAAATGCCGGAGATATTCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCg                    >  1:204124/2‑149 (MQ=255)
                                                                                                                                  aTAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcaca                                                                                                                               >  4:930298/1‑41 (MQ=255)
                                                                                                                                  aTAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcaca                                                                                                                               <  3:930298/41‑1 (MQ=255)
                                                                                                                                  aTAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGc                   >  1:362673/1‑149 (MQ=255)
                                                                                                                                  aTAAAAGTGCCGGAGAATATCTACGGCATTTTTATTCCACAGCCAAACTCACAATATATCCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATGAAACATTCGCCAGTCCGGGTAATTCACTATTCGAATTATATTATCGc                   >  4:674151/1‑149 (MQ=255)
                                                                                                                                   tAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCAc                                       >  1:691920/1‑128 (MQ=255)
                                                                                                                                   tAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCAc                                       <  2:691920/128‑1 (MQ=255)
                                                                                                                                   tAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCt                  >  3:581163/1‑149 (MQ=255)
                                                                                                                                   tAAAAATGCCGGAGATATTCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGt                                              <  2:797480/121‑1 (MQ=39)
                                                                                                                                   tAAAAATGCCGGAGATATTCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGt                                              >  1:797480/1‑121 (MQ=255)
                                                                                                                                      aaGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCg               >  2:524153/1‑149 (MQ=255)
                                                                                                                                      aaGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCg               >  2:1083239/1‑149 (MQ=255)
                                                                                                                                       aatGCCGGAGATATTCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATGCACTATTCGAATTATATTTTCGCTACGa              >  2:123659/3‑149 (MQ=255)
                                                                                                                                       aGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGa              <  3:312911/149‑1 (MQ=255)
                                                                                                                                             ggAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGTGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAAc        <  1:1083239/149‑1 (MQ=255)
                                                                                                                                              gagaATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAAcc       >  4:637245/1‑149 (MQ=255)
                                                                                                                                              gagaATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAAcc       >  1:478445/1‑149 (MQ=255)
                                                                                                                                                  atatCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGa   <  2:677354/149‑1 (MQ=255)
                                                                                                                                                   tatCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGAg  <  1:233238/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGAG  >  NC_000913/3539669‑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: