Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,512,350:1 +G 6.2% coding (65/957 nt) fecD ← ferric citrate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,512,3501.G 6.2% 353.2 / 18.4 81coding (65/957 nt)fecDferric citrate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base . (44/32);  new base G (3/2);  total (47/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

AGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAG  >  NC_000913/4512203‑4512495
                                                                                                                                                    |                                                                                                                                                 
aGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATg                                                                                                                                                  <  1:602473/149‑1 (MQ=255)
aGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATg                                                                                                                                                  >  1:306723/1‑149 (MQ=255)
  cACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGa                                                                                                                                                <  2:315953/149‑1 (MQ=255)
  cACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGa                                                                                                                                                <  4:29512/149‑1 (MQ=255)
   aCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAg                                                                                                                                               <  3:126664/149‑1 (MQ=255)
    taaccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGt                                                                                                                                              >  2:77402/4‑149 (MQ=255)
    cgcccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGt                                                                                                                                              <  4:373049/149‑1 (MQ=255)
        cgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCATATGGGAGTGAtc                                                                                                                                          <  3:276536/149‑2 (MQ=255)
                          aCCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCa                                                                                                                        >  3:497492/1‑149 (MQ=255)
                          aCCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCa                                                                                                                        <  2:608541/149‑1 (MQ=255)
                             gACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAggg                                                                                                                     >  4:239340/1‑149 (MQ=255)
                             gACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAggg                                                                                                                     >  4:515255/1‑149 (MQ=255)
                             gACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAggg                                                                                                                     >  3:614545/1‑149 (MQ=255)
                                aaaCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCaa                                                                                                                  >  1:260061/1‑149 (MQ=255)
                                aaaCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCaa                                                                                                                  >  2:332072/1‑149 (MQ=255)
                                         cagcagCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtgatga                                                                                                         >  4:287726/1‑149 (MQ=255)
                                          agcagcAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGaa                                                                                                        >  4:4145/1‑149 (MQ=255)
                                            cagcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGaaaa                                                                                                      <  4:463358/149‑1 (MQ=255)
                                            cagcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGaaaa                                                                                                      <  4:175041/149‑1 (MQ=255)
                                            cagcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGaaaa                                                                                                      <  2:177975/149‑1 (MQ=255)
                                                aaGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAc                                                                                                  <  4:214258/149‑1 (MQ=255)
                                                        cAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCaa                                                                                          >  2:12091/1‑149 (MQ=255)
                                                          gTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAAtt                                                                                        <  3:75104/149‑1 (MQ=255)
                                                              gTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCa                                                                                    >  3:556605/1‑149 (MQ=255)
                                                                  tCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCt                                                                                >  3:619369/1‑149 (MQ=255)
                                                                   ccATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTc                                                                               >  1:89153/1‑149 (MQ=255)
                                                                    cATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCg                                                                              >  1:247550/1‑149 (MQ=255)
                                                                     aTCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGc                                                                             >  2:70847/1‑149 (MQ=255)
                                                                     aTCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGc                                                                             >  2:299231/1‑149 (MQ=255)
                                                                      tCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcc                                                                            >  4:326972/1‑149 (MQ=255)
                                                                      tCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcc                                                                            >  1:316310/1‑149 (MQ=255)
                                                                         aTACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcctcc                                                                         >  3:442928/1‑149 (MQ=255)
                                                                          tACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcctcct                                                                        <  1:332072/149‑1 (MQ=255)
                                                                             ataataATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTcac                                                                     >  2:336716/1‑149 (MQ=255)
                                                                             ataataATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTcac                                                                     >  2:131031/1‑149 (MQ=255)
                                                                              taataaTGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTcaca                                                                    >  1:154803/1‑149 (MQ=255)
                                                                               aataatGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACaa                                                                   >  4:67064/1‑149 (MQ=255)
                                                                               aataatGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACaa                                                                   <  4:614545/149‑1 (MQ=255)
                                                                               aataatGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACaa                                                                   >  3:271439/1‑149 (MQ=255)
                                                                                   atGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCa                                                               <  3:143225/149‑1 (MQ=255)
                                                                                   atGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCa                                                               >  2:141703/1‑149 (MQ=255)
                                                                                   atGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCa                                                               <  1:299231/149‑1 (MQ=255)
                                                                                    tGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAg                                                              >  2:64846/1‑149 (MQ=255)
                                                                                     gCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATc                                                                                  <  2:170849/128‑1 (MQ=255)
                                                                                     gCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATc                                                                                  >  1:170849/1‑128 (MQ=255)
                                                                                      cTCGTGTCCGGCCTCCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGAc                                                            >  1:468317/1‑149 (MQ=255)
                                                                                       tCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATaa                                                                                                                                         <  3:210309/71‑1 (MQ=255)
                                                                                       tCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATaa                                                                                                                                         >  4:210309/1‑71 (MQ=255)
                                                                                       tCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACa                                                           >  2:154892/1‑149 (MQ=255)
                                                                                       tCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACa                                                           >  4:516063/1‑149 (MQ=255)
                                                                                         gtgtCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAATCCAGACaaa                                                         <  1:536007/149‑1 (MQ=255)
                                                                                             ccGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCaa                                                     <  2:89153/149‑1 (MQ=255)
                                                                                             ccGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCaa                                                     <  4:289843/149‑1 (MQ=255)
                                                                                              cGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAg                                                    <  3:183831/149‑1 (MQ=255)
                                                                                                           ggTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCa                                                               >  4:211132/1‑125 (MQ=255)
                                                                                                           ggTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCa                                                               <  3:211132/125‑1 (MQ=255)
                                                                                                                gcagcGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCa                                  >  4:85190/1‑149 (MQ=255)
                                                                                                                         gcCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTgcgc                         <  3:172049/149‑1 (MQ=255)
                                                                                                                            acggcacggGCATTACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTg                      >  2:25484/1‑149 (MQ=255)
                                                                                                                              ggcacggGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCg                    <  4:427624/149‑1 (MQ=255)
                                                                                                                                  cggGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCa                <  4:271439/149‑1 (MQ=255)
                                                                                                                                   ggggATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAg                                                    <  3:271573/112‑1 (MQ=255)
                                                                                                                                   ggggATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAg                                                    <  3:137663/112‑1 (MQ=255)
                                                                                                                                   ggggATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAg                                                    >  4:271573/1‑112 (MQ=255)
                                                                                                                                   ggggATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAg                                                    >  4:137663/1‑112 (MQ=255)
                                                                                                                                   ggggATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAg               >  1:67648/1‑149 (MQ=255)
                                                                                                                                   ggggATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAg               >  4:413931/1‑149 (MQ=255)
                                                                                                                                     ttaTCACTCCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAt             >  1:159177/3‑149 (MQ=255)
                                                                                                                                     ggATCACGCCCATATGGAGTGATAACAGCGCACAGCCCGCCTGGGCAAGGGTGATGAAAATAACCAGCGCAATGTTCATCCTCGCCGCCTCACTAGCCAGACGAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAgct             <  4:607908/149‑3 (MQ=255)
                                                                                                                                       aTCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAtct           <  4:497492/149‑1 (MQ=255)
                                                                                                                                       aTCACTCCAATATGGAGTGATAACAGCGCACAGCCCTCCAGGGCAAGGGTTATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGAGCTGCCAATCAGCGCCAGCACTACGCCTGCGGGCAGAtct           >  2:213271/1‑149 (MQ=255)
                                                                                                                                        tCACTCCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGa                                                             >  2:436825/1‑98 (MQ=255)
                                                                                                                                        tCACTCCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGa                                                             <  1:436825/98‑1 (MQ=255)
                                                                                                                                         cACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTcc         <  1:118887/149‑1 (MQ=255)
                                                                                                                                          cactccATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTcc         >  1:276020/5‑149 (MQ=255)
                                                                                                                                          aCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCg        >  4:561357/1‑149 (MQ=255)
                                                                                                                                           cTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCgg       <  4:219204/149‑1 (MQ=255)
                                                                                                                                            tCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCggg      >  2:180835/1‑149 (MQ=255)
                                                                                                                                            tCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCggg      >  2:240122/1‑149 (MQ=255)
                                                                                                                                            tCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCggg      >  3:84016/1‑149 (MQ=255)
                                                                                                                                              ccATATGTAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGa    >  2:335948/1‑149 (MQ=255)
                                                                                                                                               cATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGaa   <  2:306723/149‑1 (MQ=255)
                                                                                                                                               cATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGaa   <  3:85190/149‑1 (MQ=255)
                                                                                                                                                atatGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAg  <  1:12091/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAG  >  NC_000913/4512203‑4512495

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: