Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 623,062 A→T 5.6% M255L (ATG→TTG)  entS → enterobactin exporter, iron‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913623,0620AT 5.6% 226.8 / 3.5 72M255L (ATG→TTG) entSenterobactin exporter, iron‑regulated
Reads supporting (aligned to +/- strand):  ref base A (37/31);  new base T (2/2);  total (39/33)
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.54e-01

CGTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCATTGGTCT  >  NC_000913/622915‑623210
                                                                                                                                                   |                                                                                                                                                    
cgTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGTGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGAt                                                                                                                                                     <  3:178144/149‑1 (MQ=255)
cgTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGAt                                                                                                                                                     <  8:191694/149‑1 (MQ=255)
 gTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATg                                                                                                                                                    <  7:187172/149‑1 (MQ=255)
  tGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGt                                                                                                                                                   >  6:148210/1‑149 (MQ=255)
  tGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGt                                                                                                                                                   <  4:190199/149‑1 (MQ=255)
  ggAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTGTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGAtct                                                                                                                                                   >  1:361703/2‑147 (MQ=255)
    aGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGt                                                                                                                                                   >  1:251474/1‑147 (MQ=255)
        tCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCgg                                                                                                                                             >  6:121975/1‑149 (MQ=255)
        tCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCgg                                                                                                                                             <  7:238789/149‑1 (MQ=255)
           gTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGca                                                                                                                                           >  3:142488/1‑148 (MQ=255)
           gTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGca                                                                                                                                           >  2:76081/1‑148 (MQ=255)
             tGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTCTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGca                                                                                                                                           >  3:90718/1‑146 (MQ=255)
                aaTCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGAtt                                                                                                                                     >  2:105944/1‑149 (MQ=255)
                    aTTACTGGCTGGATTTCGTTTTCTTCTTGCGAGCCCGCTGTTGTGCGGGTATGGTCGTCTGGGTGGTTTTTTGACGATGGCGAGCGCGGGGCGGGTAGTGTATGCGGCGCTTGCTGAAAACTGGCATATGTCAGGGGCTCAGATTGGtt                                                                                                                                 <  6:153548/149‑1 (MQ=255)
                     ttACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAGCTGGCAGATGTCAGCGGCACAGATTGGttt                                                                                                                                >  2:4198/1‑149 (MQ=255)
                        cTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGCGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTtct                                                                                                                             >  5:354768/1‑149 (MQ=255)
                        cTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTtc                                                                                                                              <  3:179642/148‑1 (MQ=255)
                           gCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGAtt                                                                                                                                     >  7:314490/1‑138 (MQ=255)
                           gCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGAtt                                                                                                                                     <  8:314490/138‑1 (MQ=255)
                            cAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTAc                                                                                                                         >  7:171710/1‑149 (MQ=255)
                             aGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTAcg                                                                                                                        >  3:209427/1‑149 (MQ=255)
                                aTTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACgcgg                                                                                                                     <  2:169621/149‑1 (MQ=255)
                                aTTTCGGTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATTGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACgcgg                                                                                                                     <  5:281402/149‑1 (MQ=255)
                                  ttCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACgcggc                                                                                                                    >  3:362048/1‑148 (MQ=255)
                                   tCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGGTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGa                                                                                                                  >  1:61954/1‑149 (MQ=255)
                                    cGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGa                                                                                                                  <  2:74530/148‑1 (MQ=255)
                                       tttCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATccc                                                                                                              >  7:97096/1‑149 (MQ=255)
                                           tGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACATCTGCCAGTTGTCAGccag                                                                                                                                            >  5:288333/1‑112 (MQ=255)
                                           tGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACATCTGCCAGTTGTCAGccag                                                                                                                                            <  2:208444/115‑4 (MQ=255)
                                           tGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACATCTGCCAGTTGTCAGccag                                                                                                                                            >  1:208444/1‑112 (MQ=255)
                                           tGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACATCTGCCAGTTGTCAGccag                                                                                                                                            <  6:288333/115‑4 (MQ=255)
                                             cTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCgg                                                                                                        >  2:333201/1‑149 (MQ=255)
                                                 ccAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGcgcg                                                                                                    <  5:148210/149‑1 (MQ=255)
                                                    gCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCt                                                                                                 <  4:142488/149‑1 (MQ=255)
                                                     cccGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTa                                                                                                <  2:361988/149‑1 (MQ=255)
                                                      ccGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTAt                                                                                               >  1:390265/1‑149 (MQ=255)
                                                            gtggGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGt                                                                                           <  6:201144/147‑1 (MQ=255)
                                                              gggCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGt                                                                                       <  5:129362/149‑1 (MQ=255)
                                                               ggCGGGATTGCGCGGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGtt                                                                                      <  6:11583/149‑1 (MQ=255)
                                                                     aTTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAg                                                                                <  5:121975/149‑1 (MQ=255)
                                                                           ctgctgGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTctct                                                                                                                           >  6:201703/1‑100 (MQ=255)
                                                                           ctgctgGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTctct                                                                                                                           <  5:201703/100‑1 (MQ=255)
                                                                            tgctgGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAg                                                                         >  5:166303/1‑149 (MQ=255)
                                                                                 ggtggtTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTgg                                                                     >  2:236133/1‑148 (MQ=255)
                                                                                  gtggtTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGca                                                                   <  5:214230/149‑1 (MQ=255)
                                                                                    ggtTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTAAGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGcaca                                                                 >  1:254843/1‑149 (MQ=255)
                                                                                     gtTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGcaca                                                                 >  3:302694/1‑148 (MQ=255)
                                                                                     gtTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGAAAAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGAGGGAAGCTGGcaca                                                                 >  6:177408/1‑148 (MQ=255)
                                                                                        tATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGt                                                             >  1:89361/1‑149 (MQ=255)
                                                                                                aTGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGAc                                                      >  8:13452/1‑148 (MQ=255)
                                                                                                 tGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCa                                                    <  7:13452/149‑1 (MQ=255)
                                                                                                 tGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCa                                                    <  4:103055/149‑1 (MQ=255)
                                                                                                   gCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAgg                                                  <  2:89361/149‑1 (MQ=255)
                                                                                                    cGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAggg                                                 >  1:129809/1‑149 (MQ=255)
                                                                                                         gcgGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTAtt                                            <  1:105944/149‑1 (MQ=255)
                                                                                                         gcgGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTAtt                                            <  4:209427/149‑1 (MQ=255)
                                                                                                         gcgGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTAtt                                            >  7:15861/1‑149 (MQ=255)
                                                                                                             tGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCTCAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGTCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGAtg                                        <  8:82862/149‑1 (MQ=255)
                                                                                                                gggTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGAtgctg                                     <  4:362048/149‑1 (MQ=255)
                                                                                                                gggTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGAtgctg                                     <  2:172847/149‑1 (MQ=255)
                                                                                                                 ggTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGAtgctgc                                    >  3:337344/1‑149 (MQ=255)
                                                                                                                     cTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCtctc                                >  6:351760/1‑149 (MQ=255)
                                                                                                                              gcgcTGGCTGACAAGTGGCAGAGGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGAGCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAATCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTggg                       <  5:177408/149‑1 (MQ=255)
                                                                                                                                      tgACAACTGGCAGATTTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTcc               >  3:10041/1‑149 (MQ=255)
                                                                                                                                      tgACAACTGGCAGATGTCAGCGGCCCAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTcc               >  2:398296/1‑149 (MQ=255)
                                                                                                                                      tgACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTcc               >  3:313107/1‑149 (MQ=255)
                                                                                                                                            cTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATACCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCa         >  2:253499/1‑149 (MQ=255)
                                                                                                                                             tGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCAt        >  6:201908/1‑149 (MQ=255)
                                                                                                                                               gCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCATTg      <  8:171710/148‑1 (MQ=255)
                                                                                                                                                cAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGAtgc                                       <  5:121921/115‑1 (MQ=255)
                                                                                                                                                cAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGAtgc                                       >  6:121921/1‑115 (MQ=255)
                                                                                                                                                 aGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCATTGGt    >  2:194529/1‑149 (MQ=255)
                                                                                                                                                   aTGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCATTGGTCt  <  2:129809/149‑1 (MQ=255)
                                                                                                                                                   aTGGCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCATTGGTCt  <  1:76081/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCATTGGTCT  >  NC_000913/622915‑623210

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: