Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,511,689 G→T 5.1% A242A (GCC→GCA fecD ← ferric citrate ABC transporter membrane subunit FecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,511,6890GT 5.1% 339.7 / 2.4 101A242A (GCC→GCAfecDferric citrate ABC transporter membrane subunit FecD
Reads supporting (aligned to +/- strand):  ref base G (50/45);  new base T (3/2);  total (54/47)
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.00e+00

Only 100 of 105 total aligned reads displayed.

TGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAA  >  NC_000913/4511543‑4511838
                                                                                                                                                  |                                                                                                                                                     
tGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAggcg                                                                                                                                                    <  1:585900/150‑1 (MQ=255)
  aataatTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCGggcggc                                                                                                                                                  >  1:399573/1‑150 (MQ=255)
  aataatTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAggcggc                                                                                                                                                  <  1:897059/150‑1 (MQ=255)
      atTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACg                                                                                                                                              <  2:690798/150‑1 (MQ=255)
        tctcGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCCTGCGGcccgatt                                                                                                                                            <  2:892843/150‑6 (MQ=255)
        tctcGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGcc                                                                                                                                            <  2:1518187/150‑1 (MQ=255)
         ctcGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCCCCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCg                                                                                                                                           >  1:533836/1‑150 (MQ=255)
         ctcGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCg                                                                                                                                           >  2:663387/1‑150 (MQ=255)
         ctcGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCg                                                                                                                                           >  2:657806/1‑150 (MQ=255)
          tcGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCCTGCGGcccgattag                                                                                                                                          <  1:166327/150‑8 (MQ=255)
           cGCCAGCATATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGGCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGt                                                                                                                                         <  2:1586888/150‑1 (MQ=255)
           cGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGt                                                                                                                                         <  1:840439/150‑1 (MQ=255)
           cGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGt                                                                                                                                         <  1:531360/150‑1 (MQ=255)
              cagcagATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGa                                                                                                                                      <  1:1283355/150‑1 (MQ=255)
              cagcagATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGGCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGa                                                                                                                                      <  1:1628037/150‑1 (MQ=255)
                gcagATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATg                                                                                                                                    >  2:1001049/1‑150 (MQ=255)
                 cagATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGt                                                                                                                                   >  1:992143/1‑150 (MQ=255)
                 cagATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGt                                                                                                                                   >  1:4742/1‑150 (MQ=255)
                       ggCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGc                                                                                                                             >  2:712215/1‑150 (MQ=255)
                       ggCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGc                                                                                                                             >  1:368832/1‑150 (MQ=255)
                          aaCCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGAc                                                                                                                          <  1:233620/150‑1 (MQ=255)
                             caccaACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGc                                                                                                                       <  2:1526007/150‑1 (MQ=255)
                              accaACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAAGCGGCCACGCCGGTAGATGTCATGGCGACAGCt                                                                                                                      >  1:1154535/1‑150 (MQ=255)
                               ccaACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTa                                                                                                                     >  1:940007/1‑150 (MQ=255)
                                caACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAg                                                                                                                    >  2:1218560/1‑150 (MQ=255)
                                 aaCAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGt                                                                                                                   <  1:713598/150‑1 (MQ=255)
                                 aaCAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGt                                                                                                                   <  2:940007/150‑1 (MQ=255)
                                  aCAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTa                                                                                                                  <  1:80713/150‑1 (MQ=255)
                                    aGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAAc                                                                                                                >  1:1007896/1‑150 (MQ=255)
                                     gCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACa                                                                                                               <  1:1210715/150‑1 (MQ=255)
                                         cGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGc                                                                                                           >  2:110235/1‑150 (MQ=255)
                                         cGCACCTGTCAGGGCAGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGc                                                                                                           >  2:1206464/1‑150 (MQ=255)
                                          gCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGcc                                                                                                          <  1:954839/150‑1 (MQ=255)
                                           cACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGccc                                                                                                         >  2:254610/1‑150 (MQ=255)
                                           cACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGccc                                                                                                         >  1:1187420/1‑150 (MQ=255)
                                              cTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGa                                                                                                      >  2:294422/1‑150 (MQ=255)
                                                  cAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCg                                                                                                  <  1:84972/150‑1 (MQ=255)
                                                    gggCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCggg                                                                                                >  1:1568187/1‑150 (MQ=255)
                                                    gggCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCggg                                                                                                <  1:468136/150‑1 (MQ=255)
                                                      gCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTa                                                                                              <  2:371862/150‑1 (MQ=255)
                                                          aaaCAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGCAGCGGGTATggg                                                                                          >  1:737513/1‑150 (MQ=255)
                                                          aaaCAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATggg                                                                                          <  2:792874/150‑1 (MQ=255)
                                                           aaCAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATgggg                                                                                         <  2:368832/150‑1 (MQ=255)
                                                           aaCAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATgggg                                                                                         <  1:1415131/150‑1 (MQ=255)
                                                              aGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCAc                                                                                      >  1:1406304/1‑150 (MQ=255)
                                                                 cagcagTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCgc                                                                                   >  1:881170/1‑149 (MQ=255)
                                                                 cagcagTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGa                                                                                   >  2:956115/1‑150 (MQ=255)
                                                                 cagcagTCTGCGGAGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGCCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCgc                                                                                   >  2:1363247/1‑149 (MQ=255)
                                                                  agcagTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGac                                                                                  <  1:294422/150‑1 (MQ=255)
                                                                   gcagTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGaca                                                                                 <  1:712215/150‑1 (MQ=255)
                                                                          gCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAg                                                                          >  1:490666/1‑150 (MQ=255)
                                                                          gCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGGAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAAGCGGGTATGGGGCACCGACCCACCGAg                                                                          >  1:496197/1‑150 (MQ=255)
                                                                              tGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCgtg                                                                      <  1:572479/150‑1 (MQ=255)
                                                                                 cGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCgtggtg                                                                   >  2:1596181/1‑150 (MQ=255)
                                                                                 cGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCgtggt                                                                    >  2:406951/1‑149 (MQ=255)
                                                                                 cGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCgtggt                                                                    <  1:406951/149‑1 (MQ=255)
                                                                                  gTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGtggtgg                                                                  <  2:4742/150‑1 (MQ=255)
                                                                                  gTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGtggtgg                                                                  <  1:874237/150‑1 (MQ=255)
                                                                                    ccACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGgcg                                                                >  1:943226/1‑150 (MQ=255)
                                                                                     cACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGgcgc                                                               >  2:191873/1‑150 (MQ=255)
                                                                                       ccGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGgcgcgc                                                             >  2:666162/1‑150 (MQ=255)
                                                                                              tGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGc                                                      >  1:39898/1‑150 (MQ=255)
                                                                                               gCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGcc                                                     <  1:663387/150‑1 (MQ=255)
                                                                                               gCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGcc                                                     >  2:544826/1‑150 (MQ=255)
                                                                                                  aCGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCAcc                                                                                     >  2:194893/1‑115 (MQ=255)
                                                                                                  aCGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCAcc                                                                                     <  1:194893/115‑1 (MQ=255)
                                                                                                  aCGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAg                                                  >  2:1116174/1‑150 (MQ=255)
                                                                                                  aCGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAg                                                  <  2:1568187/150‑1 (MQ=255)
                                                                                                  aCGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCAGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAg                                                  <  1:1206464/150‑1 (MQ=255)
                                                                                                   cGCATCATATGCGGCACCACGAGAACAATAAATCTAATAGGGCCGCATGCGGCAACGCATGTAGATTTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACCCACCGAGCGTGGTGGCGCGCGCATCGCCGCGt                                                 >  2:1251809/1‑150 (MQ=255)
                                                                                                    gCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTg                                                <  2:1406304/150‑1 (MQ=255)
                                                                                                    gCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTg                                                >  1:1277867/1‑150 (MQ=255)
                                                                                                     catcatATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGc                                               <  1:1218560/150‑1 (MQ=255)
                                                                                                        catATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAag                                            >  1:402641/1‑150 (MQ=255)
                                                                                                            tGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAaggagg                                        >  1:817502/1‑150 (MQ=255)
                                                                                                               ggCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCg                                     <  2:191226/150‑1 (MQ=255)
                                                                                                                    cacGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATc                                <  1:1596181/150‑1 (MQ=255)
                                                                                                                       gagaCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATcgcg                             >  2:259526/1‑150 (MQ=255)
                                                                                                                        agaCCAATAAAGGTAATCGGGCCGCAGGCGGACACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCAAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCgg                            >  2:80073/1‑150 (MQ=255)
                                                                                                                        agaCCAATAAAGCTAATCGGGCCGCGGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCgg                            <  2:399573/150‑1 (MQ=255)
                                                                                                                        agaCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCgg                            <  2:238240/150‑1 (MQ=255)
                                                                                                                        agaCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCgg                            <  1:294144/150‑1 (MQ=255)
                                                                                                                        agaCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCgg                            <  2:679133/150‑1 (MQ=255)
                                                                                                                            cAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCaaa                        >  1:1621467/1‑150 (MQ=255)
                                                                                                                                 aaGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACCCCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTc                   >  1:347412/1‑150 (MQ=255)
                                                                                                                                  aGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTACCAACGCCCAGACTCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTcc                  >  2:442651/1‑149 (MQ=255)
                                                                                                                                  aGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCa                  <  1:635393/150‑1 (MQ=255)
                                                                                                                                  aGCTAATCGGGCCGCAGGCGGCCACGACGGTAGATGTCATGGCGACAGCTAGTAACAAAGCACAGAATCTGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCa                  >  2:1309518/1‑150 (MQ=255)
                                                                                                                                    ccggcgtggCCGCCTGAGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTTACAAAGGGCAGAATAGGGTATGGGGCGCCGACACACCGAGCGGGGTGGCGCGCGCATAGCAGAGTGCAAGGAGGTGGAGATAGCGGCAAAAACTCAgg                >  2:417543/8‑150 (MQ=255)
                                                                                                                                     tAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGc               <  1:1322336/150‑1 (MQ=255)
                                                                                                                                        cgtggCCGCCTGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAAGCGGGTCTGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTcc            >  2:164957/4‑149 (MQ=255)
                                                                                                                                         cGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAg           >  2:196389/1‑150 (MQ=255)
                                                                                                                                          gggCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGc          >  1:632340/1‑150 (MQ=255)
                                                                                                                                          gggCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGc          <  1:259526/150‑1 (MQ=255)
                                                                                                                                           ggCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCg         >  2:834705/1‑150 (MQ=255)
                                                                                                                                            gCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCgg        >  1:1357169/1‑150 (MQ=255)
                                                                                                                                            gCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCgg        >  1:1006266/1‑150 (MQ=255)
                                                                                                                                             ccGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAAc                     >  1:998236/1‑136 (MQ=255)
                                                                                                                                              cGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCa      <  1:1001049/150‑1 (MQ=255)
                                                                                                                                                  ggcggcCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGaaa  >  1:797492/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGCATCGCCGAGTGCAAGGAGGTCGAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAA  >  NC_000913/4511543‑4511838

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: