Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,517,557 T→A 6.6% E54V (GAA→GTA) ‡ fecR ← anti‑sigma transmembrane signal transducer for ferric citrate transport; periplasmic FecA‑bound ferric citrate sensor and cytoplasmic FecI ECF sigma factor activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,517,5570TA 6.6% 355.0 / 10.8 121E54V (GAA→GTA) ‡fecRanti‑sigma transmembrane signal transducer for ferric citrate transport; periplasmic FecA‑bound ferric citrate sensor and cytoplasmic FecI ECF sigma factor activator
Reads supporting (aligned to +/- strand):  ref base T (61/50);  new base A (3/5);  total (65/56)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

Only 100 of 122 total aligned reads displayed.

CGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGG  >  NC_000913/4517408‑4517701
                                                                                                                                                     |                                                                                                                                                
cGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGtttt                                                                                                                                                  <  2:464417/150‑1 (MQ=255)
 ggTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTc                                                                                                                                                 <  2:735470/150‑1 (MQ=255)
   tttCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCCCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCTCACCGCCAAGCTGGTTGCGCAGGTTTGCaa                                                                                                                                               >  1:950585/1‑150 (MQ=255)
   tttCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCaa                                                                                                                                               >  1:1217854/1‑150 (MQ=255)
     tCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAAGTTTTCAAcc                                                                                                                                             <  2:1279341/150‑1 (MQ=255)
      ccGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCCAccg                                                                                                                                            >  1:756048/1‑149 (MQ=255)
         ccTGCCAGAGCTGCCAGCCCCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGCCGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACctgct                                                                                                                                         <  1:1301168/149‑1 (MQ=255)
          cTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACctgctg                                                                                                                                        <  2:611579/150‑1 (MQ=255)
           tGCCAGAGCTGCCAGCCTCCGCCCGCGCCGCGCAACAGCAGAAATCCATTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTGAAAACctgcgca                                                                                                                                       <  2:1051245/150‑4 (MQ=255)
           tGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCtgctgc                                                                                                                                       >  1:240384/1‑150 (MQ=255)
           tGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTGAAAACctgcgca                                                                                                                                       <  1:1147229/150‑4 (MQ=255)
           tGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTGAAAACctgcgca                                                                                                                                       <  1:1153067/150‑4 (MQ=255)
            gCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTTcc                                                                                                                                      >  2:1533972/1‑150 (MQ=255)
            gCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGcc                                                                                                                                      <  2:428874/150‑1 (MQ=255)
             ccAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCa                                                                                                                                     >  1:526563/1‑150 (MQ=255)
                    tGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCAc                                                                                                                              >  1:432618/1‑150 (MQ=255)
                       cAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTgg                                                                                                                           >  1:1594939/1‑150 (MQ=255)
                         gCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGtt                                                                                                                         <  2:646804/150‑1 (MQ=255)
                          ccTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGttt                                                                                                                        <  1:602750/150‑2 (MQ=255)
                           cTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGc                                                                                                                                  <  2:681586/139‑1 (MQ=255)
                           cTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGc                                                                                                                                  >  1:681586/1‑139 (MQ=255)
                             ccGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcc                                                                                                                     <  1:1629392/150‑1 (MQ=255)
                             ccGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcc                                                                                                                     <  2:565936/150‑1 (MQ=255)
                                 cAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACCCCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCTCACTGGTTatcctgat                                                                                                                 >  2:318491/1‑150 (MQ=255)
                                 cAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctgat                                                                                                                 <  1:550582/150‑1 (MQ=255)
                                 cAGCGCCGAGCAACAGCAGTATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctgatc                                                                                                                >  2:1521499/1‑150 (MQ=255)
                                  aGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctgatc                                                                                                                <  2:1130267/150‑1 (MQ=255)
                                  aGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctgatc                                                                                                                >  2:1225965/1‑150 (MQ=255)
                                   gcgcCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctgatcc                                                                                                               <  1:1677868/150‑1 (MQ=255)
                                      ccGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGGTATCCTGATCCTGt                                                                                                            >  1:501225/1‑150 (MQ=255)
                                        gAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTc                                                                                                          <  2:1262075/150‑1 (MQ=255)
                                        gAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTc                                                                                                          >  2:542229/1‑150 (MQ=255)
                                          gCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTTAt                                                                                                        <  1:1601798/150‑1 (MQ=255)
                                             aCAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATAcc                                                                                                     <  2:944006/150‑1 (MQ=255)
                                             aCAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATAcc                                                                                                     <  2:231266/150‑1 (MQ=255)
                                               agcagTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCAc                                                                                                   >  2:1439141/1‑150 (MQ=255)
                                                gcagTATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGCTCCTGTTCATACCACTg                                                                                                 >  2:930391/1‑150 (MQ=255)
                                                 cagTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTg                                                                                                 <  2:1611824/150‑1 (MQ=255)
                                                   gTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGtt                                                                                               >  1:1496319/1‑150 (MQ=255)
                                                     aaTCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCCCTGGGTATCCTGATCCTGTTCATACCACTGTTGc                                                                                             >  1:113443/1‑150 (MQ=255)
                                                      aTCCTTTCATCACGTGACGGCTGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCTGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGcc                                                                                            <  1:1265733/150‑1 (MQ=255)
                                                         cTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCCCTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAgc                                                                                         >  2:1175401/1‑150 (MQ=255)
                                                         cTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAgc                                                                                         <  1:1225965/150‑1 (MQ=255)
                                                         cTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAgc                                                                                         <  2:432618/150‑1 (MQ=255)
                                                             catcaCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAgcgcgc                                                                                     >  2:702368/1‑150 (MQ=255)
                                                             catcaCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAgcgcgc                                                                                     >  1:1119268/1‑150 (MQ=255)
                                                               tcaCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCtt                                                                                   >  2:608560/1‑150 (MQ=255)
                                                                caCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGCTCCTGTTCATACCACTGTTGCCAGCGCGCTTc                                                                                  >  1:104521/1‑150 (MQ=255)
                                                                 aCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTcc                                                                                 <  1:552812/150‑1 (MQ=255)
                                                                 aCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTcc                                                                                 >  1:1175807/1‑150 (MQ=255)
                                                                           gggTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGc                                                                       <  2:1319419/150‑1 (MQ=255)
                                                                                aGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGAcg                                                                  >  1:949507/1‑150 (MQ=255)
                                                                                 ggCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCACCCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCCGCGCGCTTCCTGTTGTGGGCTGAcgc                                                                 >  1:1485520/1‑150 (MQ=255)
                                                                                 ggCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGAcgc                                                                 >  1:1231045/1‑150 (MQ=255)
                                                                                   cGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACgcgc                                                               <  1:1521499/150‑1 (MQ=255)
                                                                                    gggTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCGGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCGGACGCGCt                                                              <  1:1316244/150‑1 (MQ=255)
                                                                                       tATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGAGCTGACGCGCTcgc                                                           <  2:840206/150‑1 (MQ=255)
                                                                                          cGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGGTCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTcgccgc                                                        >  1:458614/1‑150 (MQ=255)
                                                                                             gCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCAAACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTa                                                     >  2:520469/1‑150 (MQ=255)
                                                                                              cAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAg                                                    >  2:1271587/1‑150 (MQ=255)
                                                                                                   ccccGCTGGCAACGTCGCCAGGCACACCGCCAAGCGGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCAGATCCTGTTCATACCACTGTTGCCCGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACgg                                               <  2:171113/150‑1 (MQ=255)
                                                                                                   cccGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACgg                                               >  2:946818/1‑150 (MQ=255)
                                                                                                      ggCTGGCAACGGCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTGCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATCCCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCGTAGCACGGCAt                                            <  2:79812/150‑1 (MQ=255)
                                                                                                      gccTGGCGACGTTGCCAGGCACACCTCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCAt                                            >  2:292915/3‑150 (MQ=255)
                                                                                                       gCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTcc                                                                                 >  2:213944/1‑112 (MQ=255)
                                                                                                       gCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTcc                                                                                 <  1:213944/112‑1 (MQ=255)
                                                                                                       gCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCACTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATa                                           >  1:330965/1‑150 (MQ=255)
                                                                                                        cTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCAGCTTAGCACGGCATAc                                          >  2:1218483/1‑150 (MQ=255)
                                                                                                         tGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATAcc                                         <  2:526563/150‑1 (MQ=255)
                                                                                                          ggCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCa                                        >  1:1143098/1‑150 (MQ=255)
                                                                                                            cAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGt                                      >  2:430400/1‑150 (MQ=255)
                                                                                                             aaCGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGATCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTg                                     <  2:161001/150‑1 (MQ=255)
                                                                                                              aCGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTgg                                    <  1:608560/150‑1 (MQ=255)
                                                                                                              aCGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTgg                                    <  1:927839/150‑1 (MQ=255)
                                                                                                                gTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGa                                  >  2:1708552/1‑150 (MQ=255)
                                                                                                                   gCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGCAAGc                               >  2:322411/1‑150 (MQ=255)
                                                                                                                     cAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTg                             <  1:57515/150‑1 (MQ=255)
                                                                                                                     cAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTg                             >  1:1391140/1‑150 (MQ=255)
                                                                                                                       ggCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGaa                           >  1:1124560/1‑150 (MQ=255)
                                                                                                                       ggCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGaa                           >  2:525854/1‑150 (MQ=255)
                                                                                                                       ggCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGaa                           >  1:408577/1‑150 (MQ=255)
                                                                                                                        gCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAAc                          <  1:1533972/150‑1 (MQ=255)
                                                                                                                        gCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAAc                          <  1:570028/150‑1 (MQ=255)
                                                                                                                           cacCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCa                       <  1:430400/150‑1 (MQ=255)
                                                                                                                           cacCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCa                       >  1:515899/1‑150 (MQ=255)
                                                                                                                            acCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAg                      >  1:1416817/1‑150 (MQ=255)
                                                                                                                            acCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCAAGTCCCCCTGGTTATCCTGATCCTGTTCATACCACTTTTGCCAGCGCGCTTCCGGTTCTGGGCTGACGCGCTATCCGCTTAGCACGGCATACCAGTTGGAAGCTGAACGCAg                      >  1:1034064/1‑150 (MQ=255)
                                                                                                                                ccAAGCTGGTTGCGCATGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCAAGGCATACCAGTGGGAAGCTGAACGCAGCGcc                  <  2:1213963/150‑1 (MQ=255)
                                                                                                                                    gCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGAc              <  1:1271587/150‑1 (MQ=255)
                                                                                                                                     cTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGCcg             >  1:1547252/1‑150 (MQ=255)
                                                                                                                                      tGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGAcgg            >  1:474715/1‑150 (MQ=255)
                                                                                                                                       ggTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGAcggc           >  1:140545/1‑150 (MQ=255)
                                                                                                                                       ggTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGAcggc           >  2:761732/1‑150 (MQ=255)
                                                                                                                                       ctggcagcAGGTTGAAAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGAcggc           >  1:1324107/7‑150 (MQ=255)
                                                                                                                                          tGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCggg        <  2:761946/150‑1 (MQ=255)
                                                                                                                                          tGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCggg        >  1:975653/1‑150 (MQ=255)
                                                                                                                                           gcgcAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGa       >  2:1562940/1‑150 (MQ=255)
                                                                                                                                            cgcAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGaa      >  1:4915/1‑150 (MQ=255)
                                                                                                                                                ggTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCgg  <  2:1594939/150‑1 (MQ=255)
                                                                                                                                                ggTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCgg  >  1:1145168/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGG  >  NC_000913/4517408‑4517701

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