Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,882,291:1 +G 6.2% coding (906/1416 nt) oprN → multidrug RND transporter outer membrane protein OprN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,882,2911.G 6.2% 261.1 / 10.5 64coding (906/1416 nt)oprNmultidrug RND transporter outer membrane protein OprN
Reads supporting (aligned to +/- strand):  ref base . (27/33);  new base G (2/2);  total (29/35)
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.03e-01

GACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCT  >  NC_002947/3882145‑3882431
                                                                                                                                                   |                                                                                                                                            
gACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCATCATAGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGTGGCCACGCCAACGTTGGCggtg                                                                                                                                          <  1:1142916/150‑4 (MQ=255)
 aCGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACAGCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTgg                                                                                                                                          <  2:555606/150‑1 (MQ=255)
  cgcgCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGc                                                                                                                                         <  2:390956/150‑1 (MQ=255)
            gTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGAc                                                                                                                               <  1:1192650/150‑1 (MQ=255)
            gTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGAc                                                                                                                               >  1:331219/1‑150 (MQ=255)
                acatGTCGCCCAAGGCCCTGCCAGCCATCGCCAGGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGt                                                                                                                           >  2:55274/4‑150 (MQ=255)
                      cGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCcgc                                                                                                                     <  2:1403869/150‑1 (MQ=255)
                         ccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCcgcgcg                                                                                                                  >  1:1269158/1‑150 (MQ=255)
                          cAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGt                                                                                                                 >  1:634338/1‑150 (MQ=255)
                            aGGCCCTGCCGGCCATCGCAAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCa                                                                                                               >  2:203118/1‑150 (MQ=255)
                               cccTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGcc                                                                                                            >  2:1253121/1‑150 (MQ=255)
                                ccTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCCAAGAGCCGAGCGCCAGCTGGCGGCGGCCTTGGCCAACATTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCt                                                                                                           >  1:674021/1‑150 (MQ=255)
                                   gCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAg                                                                                                        <  1:1050010/150‑1 (MQ=255)
                                    ccGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGt                                                                                                       <  2:690614/150‑1 (MQ=255)
                                    ccGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGt                                                                                                       <  2:1234449/150‑1 (MQ=255)
                                    ccGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGAGGAACTGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGt                                                                                                       <  1:1085616/150‑1 (MQ=255)
                                       gCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACGCCAACGTTGGCGGTGgccgccgccagctg                                                                                                                             >  2:916096/1‑120 (MQ=38)
                                       gCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACGCCAACGTTGGCGGTGgccgccgccagctg                                                                                                                             <  1:916096/125‑6 (MQ=38)
                                        ccATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCt                                                                                                   <  1:437211/150‑1 (MQ=255)
                                          aTCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTc                                                                                                 >  1:1240881/1‑150 (MQ=255)
                                          aTCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTc                                                                                                 >  2:1222342/1‑150 (MQ=255)
                                           tCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGc                                                                                                                                         >  1:1118371/1‑109 (MQ=255)
                                           tCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGc                                                                                                                                         <  2:1118371/109‑1 (MQ=255)
                                             gCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCAGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACATGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTc                                                                                              >  2:303402/1‑150 (MQ=255)
                                              ccAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCg                                                                                             >  2:602561/1‑150 (MQ=255)
                                               cAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTg                                                                                                                                           >  1:1079353/1‑103 (MQ=255)
                                               cAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTg                                                                                                                                           <  2:1079353/103‑1 (MQ=255)
                                               cAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCgg                                                                                            <  1:892830/150‑1 (MQ=255)
                                                aaGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGc                                                                                           <  1:120757/150‑1 (MQ=255)
                                                aaGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGc                                                                                           <  2:1269158/150‑1 (MQ=255)
                                                  ggCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCtt                                                                                         <  1:567703/150‑1 (MQ=255)
                                                   gcgcTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCttt                                                                                        >  1:631953/1‑150 (MQ=255)
                                                         ccGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTAccgcc                                                                                  <  1:983944/150‑1 (MQ=255)
                                                           ggtcggtGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACcgccgc                                                                                >  1:753013/1‑150 (MQ=255)
                                                                 tGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCgg                                                                          <  1:514714/150‑1 (MQ=255)
                                                                             gctgctGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGt                                                                                                                                            >  1:83717/1‑72 (MQ=255)
                                                                             gctgctGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGt                                                                                                                                            <  2:83717/72‑1 (MQ=255)
                                                                              ctgctgCGCCGCCGCCCGGACCTCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGTATCGTCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGc                                                             <  1:602561/150‑1 (MQ=255)
                                                                               tgctgcGTCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGct                                                            >  2:315332/1‑150 (MQ=255)
                                                                                   gcgcCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTc                                                        <  2:1099496/150‑1 (MQ=255)
                                                                                   gcgcCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTc                                                        <  1:1322403/150‑1 (MQ=255)
                                                                                      ccgccgccCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTcggc                                                     <  2:1240881/150‑1 (MQ=255)
                                                                                       cgccgccCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTcggcg                                                    <  2:1345567/150‑1 (MQ=255)
                                                                                         ccgccCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCcg                                                                                                                      >  2:208516/1‑82 (MQ=255)
                                                                                         ccgccCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCcg                                                                                                                      <  1:208516/82‑1 (MQ=255)
                                                                                         ccgccCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCggcggc                                                  >  1:350497/1‑150 (MQ=255)
                                                                                                   aTCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTgg                                        >  1:805018/1‑150 (MQ=255)
                                                                                                    tcagtaGCGCCGAGCGACTGCTAGCTGCGGCGACGGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTggg                                       >  2:1156531/6‑150 (MQ=255)
                                                                                                      cGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCa                                     <  1:1222342/150‑1 (MQ=255)
                                                                                                       gCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTACCCGGCTTTTCCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCAc                                    >  2:1261661/1‑150 (MQ=255)
                                                                                                        cAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTc                                                                       <  1:54678/114‑1 (MQ=255)
                                                                                                        cAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTc                                                                       >  2:54678/1‑114 (MQ=255)
                                                                                                            gcCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGc                               <  1:734716/150‑1 (MQ=255)
                                                                                                            gcCGAGCGCCAGCTCGCGGCGGCCACCGCCAATGTTGGCGTGGCCACCGCCGCCCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACTGCCGCCCGCGGCTCGCAGATCCGCCCTTCGGAGGCCAACGCCGGGGCACTGGGc                               >  2:124559/1‑150 (MQ=255)
                                                                                                                   gcCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGcc                        >  2:530436/1‑149 (MQ=255)
                                                                                                                   gcCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCa                        <  2:631953/150‑1 (MQ=255)
                                                                                                                           gcggcggcCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTgg                <  2:753013/150‑1 (MQ=255)
                                                                                                                            cggcggcCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTggg               <  2:331219/150‑1 (MQ=255)
                                                                                                                                   ccacGCCAACGTTGGCGGTGGCCGCCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGggcg             <  1:458997/143‑1 (MQ=255)
                                                                                                                                   ccacGCCAACGTTGGCGGTGGCCGCCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGggcg             >  2:458997/4‑146 (MQ=255)
                                                                                                                                    cacGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTc       >  2:1244455/3‑150 (MQ=255)
                                                                                                                                      cGCCAACGTTGGCGTGGCCACCGCCGCCCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGa     <  2:232706/150‑1 (MQ=255)
                                                                                                                                         cAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCt  <  2:585367/150‑1 (MQ=255)
                                                                                                                                         cAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCt  >  1:505898/1‑150 (MQ=255)
                                                                                                                                         cAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCt  <  2:805018/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
GACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCT  >  NC_002947/3882145‑3882431

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