Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 621,499 C→G 6.3% G231R (GGC→CGC)  fepD ← ferrienterobactin ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913621,4990CG 6.3% 298.2 / 7.7 91G231R (GGC→CGC) fepDferrienterobactin ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (44/41);  new base G (2/4);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

CCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAA  >  NC_000913/621350‑621647
                                                                                                                                                     |                                                                                                                                                    
ccAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGcc                                                                                                                                                      <  2:400540/150‑1 (MQ=255)
ccAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACAACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGcc                                                                                                                                                      <  2:87661/150‑1 (MQ=255)
  aGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCa                                                                                                                                                    >  2:1674191/1‑150 (MQ=255)
  aGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCa                                                                                                                                                    >  1:653153/1‑150 (MQ=255)
    cGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAgc                                                                                                                                                  >  2:953119/1‑150 (MQ=255)
     gatgatCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAgcg                                                                                                                                                 >  1:1444547/1‑150 (MQ=255)
      atgatCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAgcgc                                                                                                                                                >  1:937967/1‑150 (MQ=255)
      atgatCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAgcgc                                                                                                                                                >  1:931220/1‑150 (MQ=255)
       tgatCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGcccc                                                                                                                                               >  1:109235/1‑150 (MQ=255)
          tCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTc                                                                                                                                            <  2:652644/150‑1 (MQ=255)
           cGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTcg                                                                                                                                           <  2:653153/150‑1 (MQ=255)
           cGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTcg                                                                                                                                           >  1:102995/1‑150 (MQ=255)
           cGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTcg                                                                                                                                           <  1:142865/150‑1 (MQ=255)
             gCACCCACCAGCCACCGCGCCATTGGCGGCATCATCAGGCCAATAAATGAAATCGGTCCAACTATTGCCGTCGCACTACCACACCGCACGGGAATCGCCAGCAGACCACTCACAGGTGTGCGCGCCACGCGACTGCCCAGCGCCGTcgcg                                                                                                                                         <  2:662046/150‑1 (MQ=255)
                cccaccAGCCAACGCGCCATATGCGGCATCATCGGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGgtg                                                                                                                                      >  1:283461/1‑150 (MQ=255)
                cccaccAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACGGCGCTGGGCAGTcgcgtgg                                                                                                                                      <  1:1454924/150‑4 (MQ=255)
                 ccaccaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGTCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGgtgt                                                                                                                                     >  2:1448914/1‑150 (MQ=255)
                  caccaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCTGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACGGCGCTGGGCAGTcgcgtggcg                                                                                                                                    <  1:560800/150‑6 (MQ=255)
                  caccaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACGGCGCTGGGCAGTcgcgtggcg                                                                                                                                    <  1:1365561/150‑6 (MQ=255)
                     caGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCt                                                                                                                                 >  1:546943/1‑150 (MQ=255)
                       gCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCACTCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGc                                                                                                                               >  1:235165/1‑150 (MQ=255)
                       gCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGc                                                                                                                               <  2:811084/150‑1 (MQ=255)
                        ccAACGCGCCATATGCGGCATCATCAGGCCAATACAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGcc                                                                                                                              <  1:560625/150‑1 (MQ=255)
                        ccAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGcc                                                                                                                              >  2:1373938/1‑150 (MQ=255)
                         cAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCg                                                                                                                             <  1:320988/150‑1 (MQ=255)
                           aCGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTTTTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAg                                                                                                                           <  1:1159191/150‑1 (MQ=255)
                           aCGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAg                                                                                                                           <  2:94122/150‑1 (MQ=255)
                           aCGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAg                                                                                                                           <  1:682546/150‑1 (MQ=255)
                           aCGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAg                                                                                                                           >  2:902634/1‑150 (MQ=255)
                                     tGCGGCATCGTCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTg                                                                                                                 >  2:1011232/1‑150 (MQ=255)
                                     tGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTg                                                                                                                 >  2:1599173/1‑150 (MQ=255)
                                        ggCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTc                                                                                                              <  1:703590/150‑1 (MQ=255)
                                        ggCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTc                                                                                                              >  2:1679265/1‑150 (MQ=255)
                                        ggCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTc                                                                                                              <  1:1627087/150‑1 (MQ=255)
                                          catcatCAGGCCAATAAAGGCATTCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAg                                                                                                            <  1:464643/150‑1 (MQ=255)
                                                  ggCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGATAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGAc                                                                                                    >  1:1655311/1‑150 (MQ=255)
                                                     cAATAAAGGCAATCGGGCCAAGTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCa                                                                                                 <  1:848902/150‑1 (MQ=255)
                                                     cAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCa                                                                                                 >  2:1344484/1‑150 (MQ=255)
                                                     cAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCa                                                                                                 >  2:34187/1‑150 (MQ=255)
                                                       aTAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGt                                                                                               <  1:311067/150‑1 (MQ=255)
                                                        tAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAACCGGTGCAGCGCGCGACTCAGTa                                                                                              >  1:1591402/1‑150 (MQ=255)
                                                        tAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTa                                                                                              <  1:1674191/150‑1 (MQ=255)
                                                         aaaGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCCAACTGTTCAGCGCGCGACTCAGTaa                                                                                             >  2:1574518/1‑150 (MQ=255)
                                                         aaaGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTaa                                                                                             >  2:1580938/1‑150 (MQ=255)
                                                         aaaGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGGCGAGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTaa                                                                                             <  2:118151/150‑1 (MQ=255)
                                                          acggCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAAt                                                                                            >  1:37766/3‑150 (MQ=255)
                                                                 tCGGGCCAACTATTGCCGTCGCACTCCCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAg                                                                                     <  2:1142081/150‑1 (MQ=255)
                                                                  cGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGt                                                                                    >  2:154176/1‑150 (MQ=255)
                                                                  cGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGt                                                                                    <  2:597826/150‑1 (MQ=255)
                                                                   gggCCAACTATTGCCGTCGCACTACCACAAAGCACGGTACTCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCTCTGCCGAGGATCAAACAGTTCAGCGCGCGACTCAGTAATAGCGCAGtt                                                                                   >  2:1102144/1‑150 (MQ=255)
                                                                    ggCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTg                                                                                  <  1:953119/150‑1 (MQ=255)
                                                                        aaCTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTcc                                                                              <  2:931220/150‑1 (MQ=255)
                                                                        aaCTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGGCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTcc                                                                              <  2:937967/150‑1 (MQ=255)
                                                                              tGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGTCTCAGTAATAGCGCAGTTGCTCCGGCGAt                                                                        >  2:314039/1‑150 (MQ=255)
                                                                              tGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGAt                                                                        <  2:1386199/150‑1 (MQ=255)
                                                                               gCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGCTc                                                                       >  1:761539/1‑150 (MQ=255)
                                                                               gCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATc                                                                       <  2:1011464/150‑1 (MQ=255)
                                                                               gCCGTCGCACAACCACAACGCACGGTAATCGCCAGCAGACCAATCAACTGGGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATc                                                                       <  1:792861/150‑1 (MQ=255)
                                                                                ccGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCa                                                                      >  2:383295/1‑150 (MQ=255)
                                                                                 cGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCac                                                                     >  2:801838/1‑149 (MQ=255)
                                                                                        cTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTGCAGCGCGCGACTCAGGAATAGCGCAGTTGCTCCGGCGATCAGCCCCggg                                                              >  2:1517996/1‑150 (MQ=255)
                                                                                           ccACAACGCACGCTACTCGCCACCAGACCACTCACCTCTGAGCGCGCCACGCGACTGACCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAACCTGTTCAGCGCGCGACTCAGTAATAGCGCAGCTGCTCCGGCGATCAGCACCGGGATc                                                           <  1:277068/150‑1 (MQ=255)
                                                                                           ccACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATc                                                           >  2:381626/1‑150 (MQ=255)
                                                                                           ccACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATc                                                           <  1:931423/150‑1 (MQ=255)
                                                                                             acaAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCGGCACCGGGATccg                                                         >  2:946445/1‑148 (MQ=255)
                                                                                             acaAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAg                                                         >  1:792224/1‑150 (MQ=255)
                                                                                              caAAGCACGGTAATCGCCAGCAGCCCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGc                                                        >  1:231305/1‑150 (MQ=255)
                                                                                                 tgctcGGTATTCGCCAGCAGACCAATCACCTGTGGGCGCGCCACGCTACTGCGCAGCGCCGTCGCGGTGTCGCTGACGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATTGCGCAGTTGCTCCGGGGATCAGCACCGGGTTCAGgacc                                                     <  2:474053/146‑1 (MQ=255)
                                                                                                   cACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGcaccac                                                   >  2:1094718/1‑150 (MQ=255)
                                                                                                         aaTCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTaa                                             <  1:387494/150‑1 (MQ=255)
                                                                                                           tCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATATCGCAGTTGCTCCGGCGATCAGCACCGGGATCATCACCACTTTTAAgg                                           <  2:102995/150‑1 (MQ=255)
                                                                                                           tCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGGAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAgg                                           >  1:1336661/1‑150 (MQ=255)
                                                                                                            cGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGt                                          >  2:1410293/1‑150 (MQ=255)
                                                                                                             gCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGAGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTa                                         <  2:283461/150‑1 (MQ=255)
                                                                                                                agcagACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGt                                      <  1:1373938/150‑1 (MQ=255)
                                                                                                                 gcagACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTa                                     <  1:810543/150‑1 (MQ=255)
                                                                                                                      ccAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTg                                >  1:1508447/1‑150 (MQ=255)
                                                                                                                        aaTCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCg                              <  1:1599173/150‑1 (MQ=255)
                                                                                                                        aaTCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCg                              <  1:1679265/150‑1 (MQ=255)
                                                                                                                         aTCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGa                             <  2:1655311/150‑1 (MQ=255)
                                                                                                                         aTCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGa                             <  1:1344484/150‑1 (MQ=255)
                                                                                                                             aCTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAatat                         >  1:436204/1‑150 (MQ=255)
                                                                                                                                gtgtgCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCTCGACTCAGTATTATTTTTTTTTATGATTCGTTGAGCACCGTGATCTTCACCACTTTTAAGGTATGTAGATTGCGAATATCCa                      <  2:552421/150‑1 (MQ=255)
                                                                                                                                gtgtgCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCa                      >  1:1198570/1‑150 (MQ=255)
                                                                                                                                  gtgCGCGCCACGCGACTGCCCAGCGCCGTCGCGGGGTCGCTGCCGAGGCGCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGAATGCGAATATCCAGc                    <  1:1448914/150‑1 (MQ=255)
                                                                                                                                       cgcCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACGTTTAAGGTATGTAGATTGCGAATATCCAGCGAAcc               >  1:1113009/1‑150 (MQ=255)
                                                                                                                                        gcCACGCGACTGCCCAGCGCCGACGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCg              <  2:235165/150‑1 (MQ=255)
                                                                                                                                           aCGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCt           >  1:1196642/1‑150 (MQ=255)
                                                                                                                                            cgcgACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCtt          >  2:408453/1‑150 (MQ=255)
                                                                                                                                            cgcgACGGCGCTGGGCAGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCtt          >  1:32932/1‑150 (MQ=255)
                                                                                                                                            cgcgACGGCGCTGGGCAGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCtt          >  2:1385628/1‑150 (MQ=255)
                                                                                                                                                 cTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAACCTGTTCTGCGCGAGACTTATTATTTTTTTTTTGAATGATTCGTTGATCACCGTGAGCTTCACCACTTTGAAGGTAGGTAGATTGCGAATAGCCAGCGAACCGGCTGGCCAg     <  2:1289671/150‑1 (MQ=255)
                                                                                                                                                   gCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGaa   <  1:427341/150‑1 (MQ=255)
                                                                                                                                                    cccAGCGCCGTCGCGGTGGCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGaaa  <  1:34187/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAA  >  NC_000913/621350‑621647

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