Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 905,720 C→A 67.0% L270M (CTG→ATG)  amiD → 1,6‑anhydro‑N‑acetylmuramyl‑L‑alanine amidase, Zn‑dependent; OM lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913905,7200CA67.0% 158.8 / 129.0 133L270M (CTG→ATG) amiD1,6‑anhydro‑N‑acetylmuramyl‑L‑alanine amidase, Zn‑dependent; OM lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (22/22);  new base A (45/44);  total (67/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 134 total aligned reads displayed.

GCAGCGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCGGCTGAT  >  NC_000913/905623‑905819
                                                                                                 |                                                                                                   
gcagcGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTg                                                                                                   >  2:1383/1‑100 (MQ=255)
gcagcGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGAGAAGCGGATGAAGAAACTCAGGCGATTGACGAAGCATTGATg                                                                                                   >  1:2869924/1‑100 (MQ=255)
 cagcGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATgg                                                                                                  >  2:2554207/1‑100 (MQ=255)
  agcGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGga                                                                                                 <  2:2403060/100‑1 (MQ=255)
   gcggcgCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACTTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGgag                                                                                                <  2:2983164/100‑1 (MQ=255)
   gcggcgCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGTAGCATTGATgggg                                                                                                <  2:2300299/100‑3 (MQ=255)
   gcggcgCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGgag                                                                                                <  2:666227/100‑1 (MQ=255)
    cggcgCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGgaga                                                                                               <  2:2540327/100‑1 (MQ=255)
         cgTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATAc                                                                                          >  1:1724965/1‑100 (MQ=255)
         cgTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATAc                                                                                          >  2:326759/1‑100 (MQ=255)
         cgTGATTATGGCATTCCAGATGCATTTCCGCCCGACGATATATAACGGCGAAGCGGAAGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGCAATAc                                                                                          >  1:2385980/1‑100 (MQ=255)
             aTTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGc                                                                                      <  2:207644/100‑1 (MQ=255)
             aTTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGc                                                                                      >  1:2060304/1‑100 (MQ=255)
             aTTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGc                                                                                      <  1:2270477/100‑1 (MQ=255)
              ttATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCa                                                                                     <  1:1993327/100‑1 (MQ=255)
                aTGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAgg                                                                                   <  2:212138/100‑1 (MQ=255)
                aTGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAgg                                                                                   <  1:2878208/100‑1 (MQ=255)
                   gCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGTATTGCTGGAGAAATACGGGCAGGAtt                                                                                >  2:1184353/1‑100 (MQ=255)
                   gCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGAtt                                                                                <  1:3033776/100‑1 (MQ=255)
                   gCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGAtt                                                                                >  2:2911364/1‑100 (MQ=255)
                      ttCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAgc                                                                             <  2:1832035/100‑1 (MQ=255)
                        ccAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAgcgc                                                                           <  2:34069/100‑1 (MQ=255)
                         cAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGGATTGATGGAGAAATACGGGCAGGATTAGcgcg                                                                          >  1:3351125/1‑100 (MQ=255)
                         cAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGcgcg                                                                          >  1:3357495/1‑100 (MQ=255)
                          aaaTGCATTTCCGACCGACGTTATATAACGGAGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCgg                                                                         >  2:2282581/3‑100 (MQ=255)
                          aGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCgg                                                                         <  1:1710438/100‑1 (MQ=255)
                             tGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAg                                                                      <  2:596490/100‑1 (MQ=255)
                             tGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAg                                                                      >  2:3032984/1‑100 (MQ=255)
                             tGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAg                                                                      >  1:572990/1‑100 (MQ=255)
                                 tttCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAACCTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGtttt                                                                  >  1:1640384/1‑100 (MQ=255)
                                 tttCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGtttt                                                                  >  1:3036860/1‑100 (MQ=255)
                                 tttCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGtttt                                                                  <  2:882681/100‑1 (MQ=255)
                                  ttCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTc                                                                 >  2:946652/1‑100 (MQ=255)
                                  ttCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTc                                                                 <  1:1029937/100‑1 (MQ=255)
                                    ccGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCg                                                               >  1:1067286/1‑100 (MQ=255)
                                    ccGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCg                                                               <  1:499269/100‑1 (MQ=255)
                                     cGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGt                                                              >  1:2697722/1‑100 (MQ=255)
                                      gCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTg                                                             <  1:960787/100‑1 (MQ=255)
                                      gCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTg                                                             >  2:1519730/1‑100 (MQ=255)
                                      gCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTg                                                             <  2:1806051/100‑1 (MQ=255)
                                       cccGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTgg                                                            <  2:651725/100‑1 (MQ=255)
                                       cccGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTgg                                                            >  1:302898/1‑100 (MQ=255)
                                        ccGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGt                                                           >  2:2102058/1‑100 (MQ=255)
                                        ccGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGt                                                           <  2:336720/100‑1 (MQ=255)
                                         cGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTc                                                          >  1:3082285/1‑100 (MQ=255)
                                         cGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTc                                                          >  1:966496/1‑100 (MQ=255)
                                         cGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTc                                                          <  1:3070891/100‑1 (MQ=255)
                                          gACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTcg                                                         >  2:3159460/1‑100 (MQ=255)
                                          gACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTcg                                                         <  1:3159433/100‑1 (MQ=255)
                                           aCGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTcgc                                                        <  2:3259183/100‑1 (MQ=255)
                                            cGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTcgcg                                                       >  2:1121053/1‑100 (MQ=255)
                                             gTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGt                                                      <  1:2448559/100‑1 (MQ=255)
                                              ttATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTa                                                     <  1:3259781/100‑1 (MQ=255)
                                                 tataACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGCGGAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGc                                                   <  1:99010/99‑1 (MQ=255)
                                                 tataACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGcc                                                  >  1:909864/1‑100 (MQ=255)
                                                  ataACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCa                                                 >  1:328516/1‑100 (MQ=255)
                                                   taACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAg                                                >  1:96215/1‑100 (MQ=255)
                                                   taACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAg                                                <  2:2412733/100‑1 (MQ=255)
                                                    aaCGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGGCAgg                                               >  1:3305058/1‑100 (MQ=255)
                                                        gCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAggcggc                                           >  2:2158658/1‑100 (MQ=255)
                                                        gCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAggcggc                                           <  1:964044/100‑1 (MQ=255)
                                                         cGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCa                                          >  1:1241180/1‑100 (MQ=255)
                                                          gAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAg                                         >  1:2349643/1‑100 (MQ=255)
                                                          gAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAg                                         <  1:3089984/100‑1 (MQ=255)
                                                           aaGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGt                                        >  2:3007303/1‑100 (MQ=255)
                                                            aGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGtt                                       <  2:328517/100‑1 (MQ=255)
                                                             gCGGATGCAGAAACTCAGGCGATTGACGAAGCATTGATGGAGAAATACGGGCAGGATGAGCGCGGCAGTTTTCCGTGGGCGCGTAGCCAGGCGGCAGttt                                      <  2:2787159/100‑1 (MQ=255)
                                                                 aTGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTtctc                                  <  2:2054604/100‑1 (MQ=255)
                                                                 aTGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTtctc                                  >  1:3204064/1‑100 (MQ=255)
                                                                 aTGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTtctc                                  <  2:1257410/100‑1 (MQ=255)
                                                                  tGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCg                                 >  1:103728/1‑100 (MQ=255)
                                                                    cAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGAt                               <  1:1591453/100‑1 (MQ=255)
                                                                    cAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGAt                               >  1:644249/1‑100 (MQ=255)
                                                                     aGAAACTCAGGCGATTGCCGAAGCATTGATAGAGAAACACGGGCAGGTTTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATa                              <  1:2289251/100‑1 (MQ=255)
                                                                         aCTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCtt                          <  2:2352815/100‑1 (MQ=255)
                                                                         aCTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCtt                          <  1:1223014/100‑1 (MQ=255)
                                                                         aCTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCtt                          <  1:441334/100‑1 (MQ=255)
                                                                          ctgaGGCGATTGCCGAAGCATTGCTGGAGCACCACGGGCAGGAGTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCtt                          <  1:2023421/96‑1 (MQ=255)
                                                                          cTCAGGCGACTGCCGAGGCAGTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGGGGTCGCGTAGCGAGGCGGCAGTTTTCTCg                                 <  1:2802075/92‑1 (MQ=255)
                                                                           tCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCa                        <  2:1067289/100‑1 (MQ=255)
                                                                           tCAGGCGATTACCGCAGCCTTGATGGAGAAATACGGGCCGGATTAGCGCGGAAGTTTTCCGTGGTCGCGTAGCCAGGGGGTAGTTTTGTCGCTACCTTc                         >  1:143638/1‑99 (MQ=255)
                                                                            cAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCAt                       >  1:1677897/1‑100 (MQ=255)
                                                                               gCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCa                    >  2:2225119/1‑100 (MQ=255)
                                                                               gCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAACCAGACGGCAGTTTTCTCGATACCTTCATCCa                    >  2:40397/1‑100 (MQ=255)
                                                                                 gATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAgg                  <  2:3186323/100‑1 (MQ=255)
                                                                                   ttGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAggg                 <  1:1715255/99‑1 (MQ=255)
                                                                                    tGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtg               >  2:2565400/1‑100 (MQ=255)
                                                                                     gCCGAAGCCTTGATGGAGAAATACGGGCAGGCTTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGAGGCAGTTTTCTCGATACCTTCATCCAggg                 >  1:631285/1‑97 (MQ=255)
                                                                                     gCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtga              >  1:51093/1‑100 (MQ=255)
                                                                                      ccGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgat             >  2:3000559/1‑100 (MQ=255)
                                                                                       cGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgatg            <  1:2757040/100‑1 (MQ=255)
                                                                                        gAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgatga           >  1:3319715/1‑100 (MQ=255)
                                                                                        gAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgatga           >  1:3344176/1‑100 (MQ=255)
                                                                                        gAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgatga           >  1:2233625/1‑100 (MQ=255)
                                                                                        gAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgatga           >  1:277549/1‑100 (MQ=255)
                                                                                            cATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCgg       <  2:1263737/100‑1 (MQ=255)
                                                                                             aTTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCGGc      <  2:1204655/100‑1 (MQ=255)
                                                                                              ttGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCGGCt     >  2:20648/1‑100 (MQ=255)
                                                                                               tGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCGGCTg    <  2:1320917/100‑1 (MQ=255)
                                                                                                 cTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCGGCTGAt  >  2:2651452/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GCAGCGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCGGCTGAT  >  NC_000913/905623‑905819

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