Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,187,133 G→A 6.0% M58I (ATG→ATA rcnB → periplasmic modulator of Ni and Co efflux

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,187,1330GA 6.0% 217.9 / 2.6 67M58I (ATG→ATArcnBperiplasmic modulator of Ni and Co efflux
Reads supporting (aligned to +/- strand):  ref base G (27/36);  new base A (2/2);  total (29/38)
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.78e-01

GGTGCGCTTTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACG  >  NC_000913/2186987‑2187272
                                                                                                                                                  |                                                                                                                                           
ggTGCGCTTTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGTGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgt                                                                                                                                            <  6:268567/148‑1 (MQ=255)
ggTGCGCTTTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgt                                                                                                                                            <  8:51469/148‑1 (MQ=255)
          tGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGtaccgtaccg                                                                                                                                   <  4:272225/147‑1 (MQ=255)
            cTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATg                                                                                                                                <  2:64108/148‑1 (MQ=255)
                ttACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGTGGACACCGTACCAGAAATATACCTTACCGATTAGTa                                                                                                                            >  6:175994/1‑148 (MQ=255)
                ttACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATATACCGTACCGATGAGTa                                                                                                                            >  6:181795/1‑148 (MQ=255)
                 tACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGt                                                                                                                             >  8:131207/1‑146 (MQ=255)
                 tACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGt                                                                                                                             <  7:131207/146‑1 (MQ=255)
                   ccAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTacaa                                                                                                                         <  4:178492/148‑1 (MQ=255)
                    cAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTacaa                                                                                                                         >  4:310692/1‑147 (MQ=255)
                    cAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTacaa                                                                                                                         >  8:240288/1‑147 (MQ=255)
                         ccgccTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTaa                                                                                                                   <  1:167740/148‑1 (MQ=255)
                           gccTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAg                                                                                                                  >  7:204015/1‑147 (MQ=255)
                           gccTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGc                                                                                                                 <  5:96643/148‑1 (MQ=255)
                            ccTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCa                                                                                                                <  2:398420/148‑1 (MQ=255)
                                 gCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTgg                                                                                                            >  3:19107/1‑147 (MQ=255)
                                  ccGCACCAGCCAACTGGGGTTAGACCAATACCTCTGGAATTTCTAATTATGAGTTAAGTAGTTTCATTGCTGACTTTAATCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCa                                                                                                          <  6:36786/148‑1 (MQ=255)
                                         aGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCg                                                                                                    >  5:61003/1‑147 (MQ=255)
                                         aGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGATAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGa                                                                                                                                  <  7:43155/117‑1 (MQ=255)
                                         aGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGATAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGa                                                                                                                                  >  8:43155/1‑117 (MQ=255)
                                           ccACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTaa                                                                                                 <  5:289216/148‑1 (MQ=255)
                                             aCCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTAAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAAcc                                                                                               <  8:46312/148‑1 (MQ=255)
                                                 cgGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGcc                                                                                           >  6:282386/1‑148 (MQ=255)
                                                 cgGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGcc                                                                                           <  1:75801/148‑1 (MQ=255)
                                                    gTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGTATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCcgc                                                                                        >  2:58048/1‑148 (MQ=255)
                                                     ttCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCcgcg                                                                                       <  7:240288/148‑1 (MQ=255)
                                                      tCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCgcgc                                                                                      >  4:41458/1‑148 (MQ=255)
                                                        gACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGGATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGcc                                                                                     <  2:202723/147‑1 (MQ=255)
                                                        gACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGcc                                                                                     <  6:270003/147‑1 (MQ=255)
                                                         aCCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAACCCAGGGGACACCGTACCAGAAATATACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTg                                                                                   <  5:181795/148‑1 (MQ=255)
                                                         aCCAATACCTCCGTAAGTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATATACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTg                                                                                   <  5:175994/148‑1 (MQ=255)
                                                                ccTCGGGAATTTCTAAGTATGAGGTAAGTAGTTTCATTGCTGACTATAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGAGGCCg                                                                             <  6:82730/147‑1 (MQ=255)
                                                                  tCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGa                                                                          >  7:151866/1‑148 (MQ=255)
                                                                   cGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTg                                                                                                      >  8:49388/1‑119 (MQ=255)
                                                                   cGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTg                                                                                                      <  7:49388/119‑1 (MQ=255)
                                                                   cGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGAc                                                                         >  2:204803/1‑148 (MQ=255)
                                                                     ggAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGTATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGc                                                                       <  1:58048/148‑1 (MQ=255)
                                                                     ggAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGc                                                                       <  6:61003/148‑1 (MQ=255)
                                                                     ggAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGc                                                                       <  4:19107/148‑1 (MQ=255)
                                                                     ggAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGc                                                                       <  8:151866/148‑1 (MQ=255)
                                                                           tCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGa                                                                 >  8:222509/1‑148 (MQ=255)
                                                                            cTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGa                                                                 <  5:282386/147‑1 (MQ=255)
                                                                               aGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGAcc                                                               <  3:290663/146‑1 (MQ=255)
                                                                               aGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGAcc                                                               >  4:290663/1‑146 (MQ=255)
                                                                                gTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCtat                                                            <  3:310692/148‑1 (MQ=255)
                                                                                  aTGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTatat                                                          >  6:90307/1‑148 (MQ=255)
                                                                                     aGTTAATTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATgg                                                        <  1:30214/147‑1 (MQ=255)
                                                                                     aGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATggg                                                       <  3:41458/148‑1 (MQ=255)
                                                                                      gTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGt                                                      >  6:200270/1‑148 (MQ=255)
                                                                                              ttttCATTGCTGACTTTAATCATTTAAAAACAGGTGACACCGTACCAGAAATGTAACGTACCTATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTa                                              <  8:204015/147‑1 (MQ=255)
                                                                                                ttCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGcc                                                                              <  8:236614/114‑1 (MQ=255)
                                                                                                ttCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGcc                                                                              >  7:236614/1‑114 (MQ=255)
                                                                                                           cTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATTTGGGTGGCGCGTACGTGTTGATCAGc                                 >  4:11526/1‑148 (MQ=255)
                                                                                                             ttAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGa                               >  3:114359/1‑148 (MQ=255)
                                                                                                             ttAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGAACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGa                               >  8:278149/1‑148 (MQ=255)
                                                                                                              tAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGCGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCg                                <  3:324095/146‑1 (MQ=255)
                                                                                                                  cATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACAcc                           >  4:101506/1‑147 (MQ=255)
                                                                                                                    tttCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGACTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTg                                                                                   <  7:311823/89‑1 (MQ=255)
                                                                                                                     ttCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACg                       >  6:176918/1‑148 (MQ=255)
                                                                                                                     ttCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACg                       >  1:233391/1‑148 (MQ=255)
                                                                                                                        aaaCCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGt                     >  1:393856/1‑147 (MQ=255)
                                                                                                                          aCCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTatat                                                          <  6:212660/108‑1 (MQ=255)
                                                                                                                          aCCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTatat                                                          >  5:212660/1‑108 (MQ=255)
                                                                                                                             aGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATc               <  2:9490/148‑1 (MQ=255)
                                                                                                                              ggggACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGa                                      >  6:16030/1‑124 (MQ=255)
                                                                                                                              ggggACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGa                                      <  5:16030/124‑1 (MQ=255)
                                                                                                                              ggggACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCa              <  1:204803/148‑1 (MQ=255)
                                                                                                                                 gACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTa           <  5:90307/148‑1 (MQ=255)
                                                                                                                                          ccAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTacg  <  5:200270/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
GGTGCGCTTTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACG  >  NC_000913/2186987‑2187272

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