Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 621,506 C→A 5.7% T228T (ACG→ACT)  fepD ← ferrienterobactin ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913621,5060CA 5.7% 270.7 / 5.9 87T228T (ACG→ACT) fepDferrienterobactin ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (39/43);  new base A (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 = 3.60e-01

TGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACT  >  NC_000913/621357‑621654
                                                                                                                                                     |                                                                                                                                                    
tgatCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGcccc                                                                                                                                                      >  1:109235/1‑150 (MQ=255)
   tCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTc                                                                                                                                                   <  2:652644/150‑1 (MQ=255)
    cGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTcg                                                                                                                                                  >  1:102995/1‑150 (MQ=255)
    cGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTcg                                                                                                                                                  <  1:142865/150‑1 (MQ=255)
    cGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTcg                                                                                                                                                  <  2:653153/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:902634/1‑150 (MQ=255)
                    aCGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAg                                                                                                                                  <  2:94122/150‑1 (MQ=255)
                    aCGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAg                                                                                                                                  <  1:682546/150‑1 (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: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                                                                  <  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)
                                                                                                                                                    ccGTCGCGGTGGCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGAGCAGTTTCTGAGGCGTTCAGCACCGGGATCATCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACt  <  2:546943/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACT  >  NC_000913/621357‑621654

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: ‑