Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,512,350:1 +G 6.9% 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.9% 810.2 / 57.1 189coding (65/957 nt)fecDferric citrate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base . (85/91);  new base G (6/7);  total (91/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.31e-01

Only 100 of 200 total aligned reads displayed.

AGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGG  >  NC_000913/4512203‑4512496
                                                                                                                                                    |                                                                                                                                                  
aGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑g                                                                                                                                                   >  2:1463400/1‑149 (MQ=255)
aGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑g                                                                                                                                                   <  4:1267135/149‑1 (MQ=255)
aGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑g                                                                                                                                                   <  2:721772/149‑1 (MQ=255)
 gCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑gg                                                                                                                                                  >  2:1250934/1‑149 (MQ=255)
  cACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGa                                                                                                                                                 <  1:455334/149‑1 (MQ=255)
  cACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGa                                                                                                                                                 <  1:1154843/149‑1 (MQ=255)
   aCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAg                                                                                                                                                <  2:1433580/149‑1 (MQ=255)
   aCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACT‑CCATAT‑gggag                                                                                                                                               <  2:160126/149‑4 (MQ=255)
    cgcccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGt                                                                                                                                               <  2:408964/149‑1 (MQ=255)
     gcccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACT‑CCATAT‑gggagtg                                                                                                                                             <  3:1079962/149‑6 (MQ=255)
       ccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGAt                                                                                                                                            <  1:1073062/149‑1 (MQ=255)
       ccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACT‑CCATATGGGAGTGAt                                                                                                                                            <  2:963879/149‑1 (MQ=255)
          ccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAAc                                                                                                                                         >  2:1550766/1‑149 (MQ=255)
            aCGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAg                                                                                                                                       >  2:1344324/1‑149 (MQ=255)
             cGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAgc                                                                                                                                      <  2:1220149/149‑1 (MQ=255)
             cGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAgc                                                                                                                                      <  3:1334524/149‑1 (MQ=255)
                 gAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGcac                                                                                                                                  >  2:1539062/1‑149 (MQ=255)
                       tGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCg                                                                                                                            >  4:1397026/1‑149 (MQ=255)
                         cACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGcc                                                                                                                          >  1:1034397/1‑149 (MQ=255)
                          aCCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCa                                                                                                                         <  1:839023/149‑1 (MQ=255)
                            cGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGa                                                                                                                                                 >  4:801856/1‑123 (MQ=255)
                            cGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAgg                                                                                                                       >  2:701574/1‑149 (MQ=255)
                             gACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAggg                                                                                                                      >  4:493618/1‑149 (MQ=255)
                             gACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAggg                                                                                                                      <  2:927044/149‑1 (MQ=255)
                                aaaCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCaa                                                                                                                   >  3:1490672/1‑149 (MQ=255)
                                aaaCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCaa                                                                                                                   >  1:1531676/1‑149 (MQ=255)
                                 aaCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGGGTGATAACAGCGCACAGCCCGCCAGGGCAAg                                                                                                                  <  4:571126/149‑1 (MQ=255)
                                  aCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAgg                                                                                                                 <  1:887764/149‑1 (MQ=255)
                                  aCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAgg                                                                                                                 <  4:172830/149‑1 (MQ=255)
                                    aGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGt                                                                                                               <  3:1427350/149‑1 (MQ=255)
                                     gTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAg                                                                                                                                                <  2:1307042/115‑1 (MQ=255)
                                     gTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCTCACAGCCCGCCAGGGCAAGGGtg                                                                                                              <  4:1106287/149‑1 (MQ=255)
                                       gCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtgat                                                                                                            <  3:932088/149‑1 (MQ=255)
                                        ccAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtgatg                                                                                                           >  3:773313/1‑149 (MQ=255)
                                          agcagcAAGCGGGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑g                                                                                                                                                   >  4:1078258/1‑107 (MQ=255)
                                             agcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAAt                                                                                                      <  1:325033/149‑1 (MQ=255)
                                              gcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATa                                                                                                     >  1:769360/1‑149 (MQ=255)
                                              gcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCCGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑CGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATa                                                                                                     >  2:643629/1‑149 (MQ=255)
                                                aaGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAc                                                                                                   <  4:347758/149‑1 (MQ=255)
                                                 aGCGCGGCAGTGGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAcc                                                                                                  >  3:1358561/1‑149 (MQ=255)
                                                              gTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCa                                                                                     >  1:1344583/1‑149 (MQ=255)
                                                               tACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCAt                                                                                    <  3:685840/149‑1 (MQ=255)
                                                                  tCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAgc                                                                                                                                      >  4:1478986/1‑96 (MQ=255)
                                                                  tCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAAtt                                                                                         <  1:1317671/141‑1 (MQ=255)
                                                                  tCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGAGCAATTTTCATCCt                                                                                 >  2:776193/1‑149 (MQ=255)
                                                                   ccATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTc                                                                                >  4:1037719/1‑149 (MQ=255)
                                                                     aTCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGc                                                                              >  1:12635/1‑149 (MQ=255)
                                                                     aTCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGc                                                                              >  1:128537/1‑149 (MQ=255)
                                                                     aTCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGc                                                                              >  2:30959/1‑149 (MQ=255)
                                                                     aTCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGc                                                                              >  1:1163103/1‑149 (MQ=255)
                                                                      tCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcc                                                                             >  3:140066/1‑149 (MQ=255)
                                                                          tACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAAt                                                                                                      <  1:1091223/120‑1 (MQ=255)
                                                                          tACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAAt                                                                                                      >  2:1091223/1‑120 (MQ=255)
                                                                           aCATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCctcctc                                                                        <  2:924638/149‑1 (MQ=255)
                                                                            cATAATAATGCTCGTGTCCGGCCTGCCAGCCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTca                                                                       >  1:637273/1‑149 (MQ=255)
                                                                             ataataATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAAttt                                                                                        >  2:533878/1‑131 (MQ=255)
                                                                             ataataATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAAttt                                                                                        <  1:533878/131‑1 (MQ=255)
                                                                             ataataATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTcac                                                                      <  1:264066/149‑1 (MQ=255)
                                                                             ataataATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTcac                                                                      >  3:222957/1‑149 (MQ=255)
                                                                             ataataATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTcac                                                                      >  2:115835/1‑149 (MQ=255)
                                                                               aataatGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACaa                                                                    <  4:1001717/149‑1 (MQ=255)
                                                                                ataatGCTCGTGTCCGACCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAg                                                                   >  3:466754/1‑149 (MQ=255)
                                                                                    tGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAg                                                               <  2:896860/149‑1 (MQ=255)
                                                                                     gCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGt                                                                                                               <  1:1178129/100‑1 (MQ=255)
                                                                                        cGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACaa                                                           <  4:1376006/149‑1 (MQ=255)
                                                                                             ccGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCaa                                                      <  3:394477/149‑1 (MQ=255)
                                                                                             ccGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCaa                                                      <  1:1083170/149‑1 (MQ=255)
                                                                                             ccGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCaa                                                      <  1:935344/149‑1 (MQ=255)
                                                                                              cGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAg                                                     <  4:898276/149‑1 (MQ=255)
                                                                                              cGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAg                                                     <  2:765019/149‑1 (MQ=255)
                                                                                                 ccTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGc                                                  >  2:1346542/1‑149 (MQ=255)
                                                                                                     ccAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGcc                                              >  3:1151179/1‑149 (MQ=255)
                                                                                                      cAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCa                                             <  2:1163103/149‑1 (MQ=255)
                                                                                                            gTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAc                                                                                                   <  4:858405/89‑1 (MQ=255)
                                                                                                            gTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAc                                                                                                   >  3:858405/1‑89 (MQ=255)
                                                                                                            gTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCa                                         >  1:1589376/1‑147 (MQ=255)
                                                                                                             tCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAgcg                                      >  2:43153/1‑149 (MQ=255)
                                                                                                               agcagcGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAA‑TTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCa                                   >  3:445396/1‑149 (MQ=255)
                                                                                                                     gcgcgcCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACt                              <  1:200044/149‑1 (MQ=255)
                                                                                                                       gcgcCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTgc                            <  3:1064525/149‑1 (MQ=255)
                                                                                                                       gcgcCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTgc                            <  4:1151179/149‑1 (MQ=255)
                                                                                                                         gcCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAAACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTgcgc                          >  1:589901/1‑149 (MQ=255)
                                                                                                                               gcacggGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCgg                    <  3:493618/149‑1 (MQ=255)
                                                                                                                               gcacggGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCgg                    >  3:497950/1‑149 (MQ=255)
                                                                                                                               gcacggGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCgg                    <  4:466754/149‑1 (MQ=255)
                                                                                                                                   ggggATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAg                >  4:518218/1‑149 (MQ=255)
                                                                                                                                    ggggcaaCTCCCATAT‑GGAGTGCTAACAGCGCAGAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGAAATTTTAATCCTCGCCTCCTCACAAGCCAGACAAAGCAATGGCTGCCAATCAGAGCCAGCACTGCGCCCGCGGGCAGa               >  2:912420/7‑149 (MQ=255)
                                                                                                                                     ttaTCACT‑CCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcctcct                                                                         >  4:1416774/3‑90 (MQ=255)
                                                                                                                                     ttaTCACT‑CCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAg                                                     >  1:1423350/3‑110 (MQ=255)
                                                                                                                                     ttaTCACT‑CCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAt              >  4:252171/3‑149 (MQ=255)
                                                                                                                                     ggATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAt              <  2:1112807/149‑1 (MQ=255)
                                                                                                                                            tCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCATACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCAggg       >  1:418851/1‑149 (MQ=255)
                                                                                                                                            tCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCaa                                            >  1:918913/1‑112 (MQ=255)
                                                                                                                                             tccATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCggg       <  3:252171/148‑1 (MQ=255)
                                                                                                                                             cccATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCgggg      >  3:147467/1‑149 (MQ=255)
                                                                                                                                               cATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGaa    <  2:117776/149‑1 (MQ=255)
                                                                                                                                               cATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGaa    <  1:1506066/149‑1 (MQ=255)
                                                                                                                                               cATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGaa    <  3:1110178/149‑1 (MQ=255)
                                                                                                                                               cATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGaa    <  1:1060543/149‑1 (MQ=255)
                                                                                                                                                 tat‑GGAGTGATAACAGCGCACATCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAgg  >  1:967562/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGG  >  NC_000913/4512203‑4512496

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