Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,385,573 C→A 6.8% G295V (GGC→GTC)  oprD ← basic amino acid specific porin OprD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,385,5730CA 6.8% 196.8 / 4.1 67G295V (GGC→GTC) oprDbasic amino acid specific porin OprD
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base A (2/3);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.97e-01

AGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCATAGATGTTGAA  >  NC_002947/1385426‑1385718
                                                                                                                                                   |                                                                                                                                                 
aGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAgtcgccgtcgcc                                                                                                                                                  >  7:95693/1‑149 (MQ=255)
  gCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCgg                                                                                                                                                >  8:137666/1‑149 (MQ=255)
   ccTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGt                                                                                                                                               <  1:240460/149‑1 (MQ=255)
    cTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGt                                                                                                                                               >  8:50136/1‑148 (MQ=255)
       gAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCTTCGATACCATAACCATAGTCGCCGTCGCCGGTGa                                                                                                                                             <  1:241344/147‑1 (MQ=255)
       gAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACt                                                                                                                                           >  5:308858/1‑149 (MQ=255)
        aaTGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCATCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCACCGGCGACGGCGACta                                                                                                                                          <  6:9055/149‑2 (MQ=255)
         aTGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCACCGGCGACGGCGACtat                                                                                                                                         <  1:214591/149‑3 (MQ=255)
          tGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCACCGGCGACGGCGACtatg                                                                                                                                        <  7:239314/149‑4 (MQ=255)
            ccGTCTTCTGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGCATTTCTCGTCTTCGGCGTTGACGTCGGAGCGGGCCACCGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATCGTCGCCGTCGCCGGTGACTTTCt                                                                                                                                       <  8:73458/148‑1 (MQ=255)
             cGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTgg                                                                                                                                     >  8:237442/1‑149 (MQ=255)
                       gAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAggg                                                                                                                           >  3:238474/1‑149 (MQ=255)
                                gTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTg                                                                                                                                      <  8:256330/129‑1 (MQ=255)
                                gTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTg                                                                                                                                      >  7:256330/1‑129 (MQ=255)
                                gTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGTGTGAAGgtgt                                                                                                                  <  1:146231/149‑1 (MQ=255)
                                 tCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGtgtg                                                                                                                 <  2:70945/149‑1 (MQ=255)
                                 tCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCAGAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGgtgt                                                                                                                  >  5:279994/1‑148 (MQ=255)
                                      gCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAgcgc                                                                                                            >  7:210415/1‑149 (MQ=255)
                                          gCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGAt                                                                                                        <  4:238474/149‑1 (MQ=255)
                                          gCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGAt                                                                                                        <  1:154682/149‑1 (MQ=255)
                                           ccTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGAtg                                                                                                       <  3:172006/149‑1 (MQ=255)
                                              gCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTTGTAGGCCAGGGTGAAGTTGTGAGCGCCGAtgttg                                                                                                    >  8:345681/1‑149 (MQ=255)
                                               cctggATTTCTGGTCTTGGGCGTTGTAGTCGGAGGGGTCCACTGAGTTGGACAGGATGATGGTGCCGCCACCGTCTATACCATAACCATAGTCGACGTCGCCGGTGACTTTCTGGTAGGACAGGGTGATGGTGTGAGCGCCGAtgttgt                                                                                                   <  2:288501/146‑1 (MQ=255)
                                                cAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGAtgttg                                                                                                    >  5:19701/1‑147 (MQ=255)
                                                cAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGAtgttg                                                                                                    >  6:147838/1‑147 (MQ=255)
                                                cAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTa                                                                                                  <  1:306007/149‑1 (MQ=255)
                                                cAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTa                                                                                                  <  4:22239/149‑1 (MQ=255)
                                                    aTTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGgcg                                                                                              >  2:91066/1‑149 (MQ=255)
                                                         tcGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCAGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGa                                                                                         >  4:161557/1‑149 (MQ=255)
                                                           gTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGa                                                                                         >  1:164091/1‑147 (MQ=255)
                                                                cGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCt                                                                                   <  6:305251/148‑1 (MQ=255)
                                                                cGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCt                                                                                   >  1:379381/1‑148 (MQ=255)
                                                                  gCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATCCCATAACCATAGTCGCCGTCGCCGGTGACTATCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGaa                                                                                >  3:304740/1‑149 (MQ=255)
                                                                     ttGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAAcgc                                                                             >  3:18308/1‑149 (MQ=255)
                                                                       gAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAAcgcgc                                                                           >  8:242976/1‑149 (MQ=255)
                                                                       gAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAAcgcg                                                                            <  7:50136/148‑1 (MQ=255)
                                                                       gAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCAACCGTCGATACCATAACCATAGTAGCCGTAGCCGCTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTTTAGGCGCACGACAGGCTGAAcgcgc                                                                           >  4:261779/1‑149 (MQ=255)
                                                                         aGTCTGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCgg                                                                         <  6:308858/149‑1 (MQ=255)
                                                                            cGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGtt                                                                       <  7:137666/148‑1 (MQ=255)
                                                                             ggAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGccc                                                                                           >  1:317075/1‑127 (MQ=255)
                                                                             ggAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGccc                                                                                           <  2:317075/127‑1 (MQ=255)
                                                                                      cACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACTATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGt                                                            >  1:204075/1‑149 (MQ=255)
                                                                                      cACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAAcgcg                                                                            <  1:115036/133‑1 (MQ=255)
                                                                                      cACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAAcgcg                                                                            >  2:115036/1‑133 (MQ=255)
                                                                                      cACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGt                                                            >  7:154049/1‑149 (MQ=255)
                                                                                            cgtGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGcc                                                        <  5:18836/145‑1 (MQ=255)
                                                                                             ttGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCAt                                                      <  3:322161/148‑1 (MQ=255)
                                                                                             ttGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATc                                                     <  5:147838/149‑1 (MQ=255)
                                                                                                       aGATGGTGCCGCCACCGTCGATACCAAAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCgg                                           <  6:279994/149‑1 (MQ=255)
                                                                                                           ggTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACt                                       >  4:51794/1‑149 (MQ=255)
                                                                                                               ccgccACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCgg                                    <  3:243206/148‑1 (MQ=255)
                                                                                                                cgccACTGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTCAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCgg                                    >  6:227186/1‑147 (MQ=255)
                                                                                                                cgccACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCa                                  <  8:210415/149‑1 (MQ=255)
                                                                                                                cgccACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCa                                  <  1:46463/149‑1 (MQ=255)
                                                                                                                 gccACCGTAGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTTTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCag                                 <  7:345681/149‑1 (MQ=255)
                                                                                                                    aCCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGAt                                                                                                        <  2:321546/75‑1 (MQ=255)
                                                                                                                    aCCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGAt                                                                                                        >  1:321546/1‑75 (MQ=255)
                                                                                                                          gATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCg                         <  3:51794/148‑1 (MQ=255)
                                                                                                                          gATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGc                        <  4:18308/149‑1 (MQ=255)
                                                                                                                              ccataaccataGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCtt                    <  3:161557/149‑1 (MQ=255)
                                                                                                                               cataaccataGTCGCCGTCGACGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGAGAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTg                   >  8:166872/1‑149 (MQ=255)
                                                                                                                                  aaccataGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCtt                    >  1:19808/1‑145 (MQ=255)
                                                                                                                                  aaccataGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCtt                    <  2:19808/145‑1 (MQ=255)
                                                                                                                                   cagaaaGTCACCGGCGACGGCGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGgtgt               >  2:105519/6‑149 (MQ=255)
                                                                                                                                    ccataGGCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTc              <  5:64015/149‑1 (MQ=255)
                                                                                                                                     cataGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTg                   <  4:198270/143‑1 (MQ=255)
                                                                                                                                     cataGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTg                   >  3:198270/1‑143 (MQ=255)
                                                                                                                                      ataGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCa             >  1:228976/1‑148 (MQ=255)
                                                                                                                                          tcgccggcgccgGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCTGATATGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCATAgtt        <  7:242976/149‑3 (MQ=255)
                                                                                                                                              cgtcgccgGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCATAGATGtt     >  5:247274/1‑148 (MQ=255)
                                                                                                                                              cgtcgccgGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTAGTCCAGCTTGTCGCCATCGAAGGCGCGGTACACGGCAGAGCCTTCGCTCTTGGTGTCATAGAAGtt     >  5:247271/1‑148 (MQ=255)
                                                                                                                                                tcgccgGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCATAGATGTTGaa  <  8:95693/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCATAGATGTTGAA  >  NC_002947/1385426‑1385718

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