Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 623,054 A→T 6.4% N252I (AAC→ATC)  entS → enterobactin exporter, iron‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913623,0540AT 6.4% 197.0 / 3.8 68N252I (AAC→ATC) entSenterobactin exporter, iron‑regulated
Reads supporting (aligned to +/- strand):  ref base A (34/29);  new base T (2/2);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

CCGCGTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCC  >  NC_000913/622912‑623197
                                                                                                                                              |                                                                                                                                               
ccGCGTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTgg                                                                                                                                             >  6:214230/1‑147 (MQ=255)
  gcgTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTgg                                                                                                                                             >  7:111347/1‑145 (MQ=255)
   cgTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGTGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGAt                                                                                                                                        <  3:178144/149‑1 (MQ=255)
   cgTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGAt                                                                                                                                        <  8:191694/149‑1 (MQ=255)
    gTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATg                                                                                                                                       <  7:187172/149‑1 (MQ=255)
     tGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGt                                                                                                                                      <  4:190199/149‑1 (MQ=255)
     tGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGt                                                                                                                                      >  6:148210/1‑149 (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                                                                                                                              >  2:76081/1‑148 (MQ=255)
              gTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGca                                                                                                                              >  3:142488/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                                                                                                                               <  6:288333/115‑4 (MQ=255)
                                              tGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACATCTGCCAGTTGTCAGccag                                                                                                                               >  1:208444/1‑112 (MQ=255)
                                              tGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACATCTGCCAGTTGTCAGccag                                                                                                                               >  5:288333/1‑112 (MQ=255)
                                              tGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACATCTGCCAGTTGTCAGccag                                                                                                                               <  2:208444/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                               >  7:15861/1‑149 (MQ=255)
                                                                                                            gcgGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTAtt                               <  4:209427/149‑1 (MQ=255)
                                                                                                            gcgGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTAtt                               <  1:105944/149‑1 (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)
                                                                                                                                              |                                                                                                                                               
CCGCGTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCC  >  NC_000913/622912‑623197

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: