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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913623,0580GC 5.9% 231.4 / 5.7 68W253C (TGG→TGCentSenterobactin exporter, iron‑regulated
Reads supporting (aligned to +/- strand):  ref base G (35/29);  new base C (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 = 7.67e-01

CCGCGTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCATTG  >  NC_000913/622912‑623206
                                                                                                                                                  |                                                                                                                                                    
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                                                                                                                                               >  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                                                                                                                                         <  7:238789/149‑1 (MQ=255)
           tCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCgg                                                                                                                                         >  6:121975/1‑149 (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                                                                                                                                 <  8:314490/138‑1 (MQ=255)
                              gCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGAtt                                                                                                                                 >  7:314490/1‑138 (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                                                                                                                                        >  1:208444/1‑112 (MQ=255)
                                              tGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACATCTGCCAGTTGTCAGccag                                                                                                                                        <  2:208444/115‑4 (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                                                                                                                       <  5:201703/100‑1 (MQ=255)
                                                                              ctgctgGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTctct                                                                                                                       >  6:201703/1‑100 (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                                                <  4:103055/149‑1 (MQ=255)
                                                                                                    tGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCa                                                <  7:13452/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)
                                                                                                                                               cTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATACCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCa     >  2:253499/1‑149 (MQ=255)
                                                                                                                                                tGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCAt    >  6:201908/1‑149 (MQ=255)
                                                                                                                                                  gCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCATTg  <  8:171710/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCGCGTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTCGCCATTG  >  NC_000913/622912‑623206

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: