Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 622,455 G→C 5.9% W52C (TGG→TGC entS → enterobactin exporter, iron‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913622,4550GC 5.9% 365.8 / 12.5 119W52C (TGG→TGCentSenterobactin exporter, iron‑regulated
Reads supporting (aligned to +/- strand):  ref base G (52/60);  new base C (3/4);  total (55/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-01

Only 100 of 122 total aligned reads displayed.

AAACAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGCCTTA  >  NC_000913/622306‑622603
                                                                                                                                                     |                                                                                                                                                    
aaaCAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTgg                                                                                                                                                      >  2:1791566/1‑150 (MQ=255)
  aCAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCa                                                                                                                                                    <  1:710371/150‑1 (MQ=255)
        cTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTggg                                                                                                                                              >  1:485540/1‑150 (MQ=255)
         tGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTgggg                                                                                                                                             <  1:1022394/150‑1 (MQ=255)
           gctgctTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTggaatg                                                                                                                                           <  1:515563/150‑5 (MQ=255)
           gctgctTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTggaatg                                                                                                                                           <  1:416443/150‑5 (MQ=255)
           gctgctTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCACTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTggaatg                                                                                                                                           <  2:1593569/150‑5 (MQ=255)
           gctgctTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCACTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTggaatg                                                                                                                                           <  2:473102/150‑5 (MQ=255)
           gctgctTAACCTCAGCCTGGTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGAAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCt                                                                                                                                           <  1:93581/150‑1 (MQ=255)
            ctgctTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCtt                                                                                                                                          <  2:89717/150‑1 (MQ=255)
                   aCCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGGTTTCGGTGa                                                                                                                                   >  1:1704550/1‑150 (MQ=255)
                   aCCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGa                                                                                                                                   >  1:589959/1‑150 (MQ=255)
                   aCCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGa                                                                                                                                   >  1:108539/1‑150 (MQ=255)
                    ccTCAGCCTGTTGAAACCGCACCCGGCGTTTCGCGCACTAATCATCGCTCGTTTCATCTCAATTGAGGCTCTGGGTTTGCTCGGCGTCGCGGTGCCTGTGCAGATCCAGATGATGAAACATTCAAGCTGGCAGGTGGGGCTTTCGGTGAc                                                                                                                                  >  1:1388869/1‑150 (MQ=255)
                    ccTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCACTTGTGGCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGAc                                                                                                                                  <  2:1782792/150‑1 (MQ=255)
                    ccTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGAc                                                                                                                                  <  2:277009/150‑1 (MQ=255)
                    ccTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGAc                                                                                                                                  <  2:1556015/150‑1 (MQ=255)
                     cTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACg                                                                                                                                 <  2:1045743/150‑1 (MQ=255)
                       cAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACTTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCt                                                                                                                               >  2:519395/1‑150 (MQ=255)
                        aGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTAATCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTg                                                                                                                              <  2:1250652/150‑1 (MQ=255)
                         gCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGa                                                                                                                             >  2:1441223/1‑150 (MQ=255)
                         gCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGa                                                                                                                             >  2:1001306/1‑150 (MQ=255)
                           cTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCAGGCAGGGGGGGCTTTCGGTGCCGCTGAcc                                                                                                                           >  2:228071/1‑150 (MQ=255)
                            tGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACcg                                                                                                                          >  1:922599/1‑150 (MQ=255)
                               tGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACcggcg                                                                                                                       <  2:590406/150‑1 (MQ=255)
                               tGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACcggcg                                                                                                                       >  2:421135/1‑150 (MQ=255)
                                tgaaaCGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATACAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACcggcgg                                                                                                                      >  2:1102636/3‑150 (MQ=255)
                                 aaaaCGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGGCTCTGGGTTGGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCggcggc                                                                                                                     <  1:1268445/150‑1 (MQ=255)
                                    aCGGACCCGGCGTTGCGCGCAGTAGTCCGCGCGCGCTGCATCTCAAGTGTGTCGCGGGGTTTGCTCGGCGGCGCGGGGCCGGTGCAGATCCAGATGAGGACACATTCCCCCTGGCAGGTGGGGCTTTCGGTGACTCTGACCGGCGGcgcg                                                                                                                  <  1:266778/150‑1 (MQ=255)
                                     cGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGa                                                                                                                 <  2:626907/150‑1 (MQ=255)
                                       cccccGGCGTGACGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGGCGCGGGGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtg                                                                                                               <  1:989190/148‑1 (MQ=255)
                                       cACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtg                                                                                                               <  2:1342356/150‑1 (MQ=255)
                                       cACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtg                                                                                                               >  1:1242188/1‑150 (MQ=255)
                                       cACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTCTCATCTCAATTGTCTCTCTGGGTGTGCTCGGCGGCGCGGTGCCGGGGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtg                                                                                                               <  1:118909/150‑1 (MQ=255)
                                               gTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTgg                                                                                                       <  2:730388/150‑1 (MQ=255)
                                               gTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTgg                                                                                                       <  2:1690844/150‑1 (MQ=255)
                                               gTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGGCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTgg                                                                                                       <  2:360072/150‑1 (MQ=255)
                                                    gcgcAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCAGACCGGCGGCGCGATGTTTGTTGGccggg                                                                                                  >  1:1151937/1‑147 (MQ=255)
                                                       cAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGTCGGCGCGATGTTTGTTGGCCTGATgg                                                                                               >  1:1384988/1‑150 (MQ=255)
                                                           aTTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGGGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTcgg                                                                                           <  2:1434438/150‑1 (MQ=255)
                                                            ttCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTcggc                                                                                          >  2:1015123/1‑150 (MQ=255)
                                                             tCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTcggcg                                                                                         <  1:793993/150‑1 (MQ=255)
                                                              ccTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTcggcgg                                                                                        <  1:523435/150‑1 (MQ=255)
                                                              ccTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTcggcgg                                                                                        <  1:1144855/150‑1 (MQ=255)
                                                               ctcgctcgTTTCATCTCAAATGAGTCTCTGGGATTGCTCGGCGTCGCGGGGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGGGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGgt                                                                                       <  1:80929/150‑1 (MQ=255)
                                                                tcgctcgTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGTGCGATGTTTGTTGGCCTGATGGTCGGCGgtg                                                                                      <  1:1791566/150‑1 (MQ=255)
                                                                tcgcgcgTTTCATCTCAATTGGGTCTCTGGGTTGGCTCGGCGTCGCGGTGCCGGGGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGgtg                                                                                      <  1:48782/150‑1 (MQ=255)
                                                                     cgTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTg                                                                                 >  1:1535502/1‑150 (MQ=255)
                                                                        ttCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCg                                                                              >  1:1709206/1‑150 (MQ=255)
                                                                        tgcATCTCCCTTGTGCCGCTGGGTTGGCTCGGCGGCGCGGCGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGGGTGCTGGCg                                                                              <  1:111964/148‑1 (MQ=255)
                                                                          cATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGTGGCTTTCGGTGCCGCTGACCGGCTGCGCGATGTATGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGa                                                                            >  2:96362/1‑150 (MQ=255)
                                                                            tctcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATc                                                                          <  2:798117/150‑1 (MQ=255)
                                                                            tctcAAATCTGTCTCTGGGATTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATc                                                                          <  2:1098310/150‑1 (MQ=255)
                                                                             ctcAATTGTGTCGCTGGGTGTGCGCGGCGGCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGGTGGCCTGATGGGCGGCGGTGTGCTGGCGGATCg                                                                         <  2:741005/150‑1 (MQ=255)
                                                                               cAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCt                                                                       >  2:659235/1‑150 (MQ=255)
                                                                                    gtgtCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAg                                                                  >  2:597131/1‑150 (MQ=255)
                                                                                     ggtCTCTGGGTTTGCGCGGCGGCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAgc                                                                 <  2:342040/149‑1 (MQ=255)
                                                                                        ctctGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCa                                                              >  1:1334549/1‑150 (MQ=255)
                                                                                             ggTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaaaaaa                                                         >  1:445882/1‑150 (MQ=255)
                                                                                              gTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAg                                                        <  2:1014062/150‑1 (MQ=255)
                                                                                              gTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAg                                                        >  2:769053/1‑150 (MQ=255)
                                                                                              gTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTCGCGGATCGCTATGAGCGCAAAAAAg                                                        <  1:894034/150‑1 (MQ=255)
                                                                                                ttGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTg                                                      <  2:108539/150‑1 (MQ=255)
                                                                                                  gCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGAt                                                    >  1:256308/1‑150 (MQ=255)
                                                                                                    tCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGAttt                                                  >  1:1010314/1‑150 (MQ=255)
                                                                                                    gcGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGAttt                                                  <  1:771330/149‑1 (MQ=255)
                                                                                                      ggCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTg                                                >  1:1719268/1‑150 (MQ=255)
                                                                                                       gcggcgCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGc                                               <  1:676458/150‑1 (MQ=255)
                                                                                                           cgcgGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGgc                                           >  2:1340235/1‑150 (MQ=255)
                                                                                                              ggtgccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGgcgcg                                        >  2:920361/1‑150 (MQ=255)
                                                                                                              ggtgccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGgcgcg                                        >  2:872880/1‑150 (MQ=255)
                                                                                                               gtgccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGgcgcgc                                       >  1:152328/1‑150 (MQ=255)
                                                                                                                  ccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGGGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGc                                    <  1:1703632/150‑1 (MQ=255)
                                                                                                                     gtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCAcc                                 >  2:509171/1‑150 (MQ=255)
                                                                                                                        cAGTTCCAGATGATGACACATTCCACCTGGGAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGAGAATATGCTGGCGCGCGGCACCtgt                              <  2:1518160/150‑1 (MQ=255)
                                                                                                                         aGATCCAGATGATGACACATTCCCCCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCtgtg                             <  2:1472355/150‑1 (MQ=255)
                                                                                                                          gATCCAGATGATGCCACATTCCTCCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAGGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTgg                            <  1:41996/150‑1 (MQ=255)
                                                                                                                          gATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGGCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTgg                            <  2:1691019/150‑1 (MQ=255)
                                                                                                                          gATCCAGATGAGGACACATTCCACCTGGGAGGGGGGGCTTTCGGTGCCGCTGACCGGCGGCGCGAGGTTTGTTGGCCTGATGGTCGGCGGTGTGGTGGCGGATCGCTTTGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTgg                            <  1:96362/150‑1 (MQ=255)
                                                                                                                             ccAGATGATGACACATTACACCTGGCAGGTGGCGCTTTCGGTGACGCTGACCGGCGGCGAGATGTTTGTAGGCCTGATGTTCGGCGGTGTGCAGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTCCTGACGCTCGGCACCTGTGGCAt                         >  2:113121/1‑150 (MQ=255)
                                                                                                                                 atgatgACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGc                     <  1:155561/150‑1 (MQ=255)
                                                                                                                                    atgaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTCACCGGCGGCGCGATGTTAGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATATTGCTGGCGCGCGGCACCTGTGGCATTGGCTTc                  >  1:498576/1‑150 (MQ=255)
                                                                                                                                    aggccacaTTCCACCTGGCAGCTGGGGCCTGCGGTGACGAGGACCGGCGGCGCGATGTCGGTTGGCCTGATGGGCGCCGGTGTGATGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCGGGCGCGCGGCACCGGTGGCATTGGCTTc                  <  2:289632/146‑1 (MQ=255)
                                                                                                                                     tgaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCTCGGCACCTGTGGTATTGGCTTCa                 >  1:26900/1‑150 (MQ=255)
                                                                                                                                       acacATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCAtt               <  2:589959/150‑1 (MQ=255)
                                                                                                                                        cacaTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTg              <  1:1277374/150‑1 (MQ=255)
                                                                                                                                           gccccACCTGCCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGAc           >  1:663515/4‑150 (MQ=255)
                                                                                                                                           aTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTgggc           <  1:1102636/150‑3 (MQ=255)
                                                                                                                                           aTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGAc           >  2:435081/1‑150 (MQ=255)
                                                                                                                                            ttCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACt          <  1:645757/150‑1 (MQ=255)
                                                                                                                                            ccccACCTGCCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACt          >  1:994975/3‑150 (MQ=255)
                                                                                                                                            ccccACCTGCCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGGTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAACAAGTGATTTGGCTGGCGCGCGGCCCCTGTGGCAGGGGCTGCATTGGAcg          >  2:153196/3‑149 (MQ=255)
                                                                                                                                              ccACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACtgt        >  1:1170561/1‑150 (MQ=255)
                                                                                                                                               cACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACtgtg       >  1:427175/1‑150 (MQ=255)
                                                                                                                                               cACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACtgtg       >  2:204213/1‑150 (MQ=255)
                                                                                                                                               cACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACtgtg       >  1:1043292/1‑150 (MQ=255)
                                                                                                                                               cACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTGGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATAGGCTTCATTGGACtgtg       <  2:342503/150‑1 (MQ=255)
                                                                                                                                                ccgtGGCAGGTGGGGCTTTCGGTGACGCGTTCGTGCGGCTCGATGGTTGTGTGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTTGCTTCATTGGACTGTGc      <  1:774649/147‑1 (MQ=255)
                                                                                                                                                   tGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGCCtt   <  2:1380258/150‑1 (MQ=255)
                                                                                                                                                    ggCAGGTGGGGCTTTGGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGCCTTa  <  1:228071/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAACAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGCCTTA  >  NC_000913/622306‑622603

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: