Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 621,503 C→G 5.7% A229A (GCG→GCC fepD ← ferrienterobactin ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913621,5030CG 5.7% 293.1 / 7.0 87A229A (GCG→GCCfepDferrienterobactin ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (39/43);  new base G (2/3);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

CGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACG  >  NC_000913/621354‑621649
                                                                                                                                                     |                                                                                                                                                  
cGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAgc                                                                                                                                                    >  2:953119/1‑150 (MQ=255)
 gatgatCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAgcg                                                                                                                                                   >  1:1444547/1‑150 (MQ=255)
  atgatCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAgcgc                                                                                                                                                  >  1:931220/1‑150 (MQ=255)
  atgatCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAgcgc                                                                                                                                                  >  1:937967/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                                                                                                                             <  1:682546/150‑1 (MQ=255)
                       aCGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAg                                                                                                                             <  2:94122/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:1627087/150‑1 (MQ=255)
                                    ggCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTc                                                                                                                <  1:703590/150‑1 (MQ=255)
                                    ggCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTc                                                                                                                >  2:1679265/1‑150 (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:597826/150‑1 (MQ=255)
                                                              cGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGt                                                                                      >  2:154176/1‑150 (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                                                             <  1:931423/150‑1 (MQ=255)
                                                                                       ccACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATc                                                             >  2:381626/1‑150 (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:1679265/150‑1 (MQ=255)
                                                                                                                    aaTCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCg                                <  1:1599173/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            >  2:1385628/1‑150 (MQ=255)
                                                                                                                                        cgcgACGGCGCTGGGCAGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCtt            >  1:32932/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)
                                                                                                                                                  cAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGCTCAGCACCACTGTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACg  >  1:1523107/1‑150 (MQ=255)
                                                                                                                                                  cAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACg  >  1:1415726/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACG  >  NC_000913/621354‑621649

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: