Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,241,025 C→T 12.1% V240M (GTG→ATG)  PP_3718 ← aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,241,0250CT12.1% 190.6 / 14.2 74V240M (GTG→ATG) PP_3718aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (35/30);  new base T (4/5);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TCTGGGGTGACGCCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTG  >  NC_002947/4240877‑4241173
                                                                                                                                                     |                                                                                                                                                    
tCTGGGGTGACGCCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGTGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCAcc                                                                                                                                                     >  2:299675/1‑150 (MQ=255)
tCTGGGGTGACGCCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCAcc                                                                                                                                                     <  1:499659/150‑1 (MQ=255)
tCTGGGGTGACGCCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCGAGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCAcc                                                                                                                                                     >  1:1859508/1‑150 (MQ=255)
    gggTGACGCCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACcccgc                                                                                                                                                 >  2:875953/1‑150 (MQ=255)
     ggTGACGCCAAACAGTTCCTCGCTGGCGAAGCAGCTGCCCAGGCGCGGCCGGCATGGCGC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGCTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACcccgcc                                                                                                                                                <  2:989331/150‑1 (MQ=255)
          cGCCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCGC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGATGGAAGTAGCCCTCGGGCGGGGTGATCATTCCGCCCGAgg                                                                                                                                           <  2:741698/150‑2 (MQ=255)
             cAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGATGCGAC‑CGCTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCa                                                                                                                                        >  2:1709091/1‑150 (MQ=255)
              aaaCAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCAG‑CACTTAGTCGGCGACGAACAGGGGGACGGAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCAt                                                                                                                                       >  1:233182/1‑150 (MQ=255)
                       cTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGCACAAAG‑CACGTCCTTCCTCGGCTACGAACAGGATGTCGTCGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGAtcggctc                                                                                                                              <  1:1867355/150‑1 (MQ=255)
                        tCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGAtcggctcg                                                                                                                             >  2:1627651/1‑150 (MQ=255)
                            gggCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCga                                                                                                                         <  1:269497/149‑1 (MQ=255)
                             ggCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCGATCAtgggctcgggcgg                                                                                                                        <  1:1101042/150‑5 (MQ=255)
                                gaagaaGGTGCCCAGGCGCCCGAACGAGGTCAC‑CGCTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCgaagaa                                                                                                                     >  1:980217/1‑150 (MQ=255)
                                    aagGTGCCCAGGCGCCCGAACGAGGTCAC‑CGCTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATgcc                                                                                                                 >  2:1120194/1‑150 (MQ=255)
                                    aagGTGCCCAGGCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATgcc                                                                                                                 <  2:1731305/150‑1 (MQ=255)
                                    aagGTGCCCAGGCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATgcc                                                                                                                 <  1:103499/150‑1 (MQ=255)
                                     agGTGCCCAGGCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATgccg                                                                                                                <  1:1594734/150‑1 (MQ=255)
                                      ggggCCCAGGCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATgccgc                                                                                                               <  1:26673/147‑1 (MQ=255)
                                      ggTGCCCAGGCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGACTGATCACCCCGCCCGCGCCCATGATCGGCTCGGCGAAGAATgccgc                                                                                                               <  1:337661/150‑1 (MQ=255)
                                        tGCCCAGGCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCAGCGAAGAAGggcgccg                                                                                                             >  1:122422/1‑149 (MQ=255)
                                         gCCCAGGCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCAc                                                                                                            <  2:5366/150‑1 (MQ=255)
                                           ccAGGCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGt                                                                                                          >  2:372731/1‑150 (MQ=255)
                                              ggCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGt                                                                                                       >  1:775806/1‑150 (MQ=255)
                                               gcgcCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTc                                                                                                      <  1:271411/150‑1 (MQ=255)
                                               gcgcCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTc                                                                                                      >  2:23266/1‑150 (MQ=255)
                                                   ccGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCa                                                                                                  >  2:1807945/1‑150 (MQ=255)
                                                       aCGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCa                                                                                                  >  1:1191681/1‑146 (MQ=255)
                                                       aCGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCa                                                                                                  <  2:1191681/146‑1 (MQ=255)
                                                        cGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAg                                                                                             >  1:1582825/1‑150 (MQ=255)
                                                         gAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGGGTGATCATTCCGCCCGAgg                                                                                                                                           >  2:743753/1‑102 (MQ=38)
                                                         gAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGGGTGATCATTCCGCCCGAgg                                                                                                                                           <  1:743753/103‑2 (MQ=38)
                                                             tCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGGATGCCGCCACGTTGTCGGCACCCAGCgaga                                                                                        >  1:1436434/1‑150 (MQ=255)
                                                                    cTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCt                                                                                  >  2:1443648/1‑150 (MQ=255)
                                                                     ttCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCtt                                                                                 <  1:1595251/150‑1 (MQ=255)
                                                                        gtcgGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCTACATGCGCTGGAAGTAGCCCTCAGGCGGAGCGATCCCCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTc                                                                              >  1:1076569/1‑150 (MQ=255)
                                                                        gtcgGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTc                                                                              <  1:1443648/150‑1 (MQ=255)
                                                                              gacgaACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGGGTGATCATTCCGCCCGAggg                                                                                                                                          >  2:1560096/1‑82 (MQ=25)
                                                                              gacgaACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGGGTGATCATTCCGCCCGAggg                                                                                                                                          <  1:1560096/84‑3 (MQ=25)
                                                                                 gaACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTc                                                                     >  1:1597760/1‑150 (MQ=255)
                                                                                 gaACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCAGCTCGGCGAAGAATGTAGCCTCTTTGTCGGCACCTACCGATAGTTTCTTGTCCTCGAACTc                                                                     >  2:116468/1‑150 (MQ=255)
                                                                                  aaCAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCg                                                                    <  2:1440956/150‑1 (MQ=255)
                                                                                   aCAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCgg                                                                   <  1:875953/150‑1 (MQ=255)
                                                                                      ggATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCATCACACGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGcca                                                                <  1:116468/150‑1 (MQ=255)
                                                                                       gATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGccac                                                               >  1:1163128/1‑150 (MQ=255)
                                                                                          gTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAg                                                            >  1:767752/1‑150 (MQ=255)
                                                                                             gTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGGGTGATCATTCCGCCCGAgg                                                                                                                                           >  2:1602261/1‑67 (MQ=21)
                                                                                             gTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGGGTGATCATTCCGCCCGAgg                                                                                                                                           <  1:1602261/68‑2 (MQ=21)
                                                                                             gTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGt                                                                                                       >  1:567006/1‑104 (MQ=255)
                                                                                             gTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGt                                                                                                       <  2:567006/104‑1 (MQ=255)
                                                                                             gTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGaa                                                         >  2:1262513/1‑150 (MQ=255)
                                                                                                  tCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGa                                                    <  1:372731/150‑1 (MQ=255)
                                                                                                   cTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGAACGGTTCAGCCAGGAATGCTGATCTCTTGTCGGTTTCTATCTAGATGGTCTGGTCCTCTAACTCGGCCACCAGGAAGTCGAg                                                   <  2:1162596/150‑1 (MQ=255)
                                                                                                       cacaGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAAc                                               <  2:775806/150‑1 (MQ=255)
                                                                                                         caGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCCCGCTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTc                                             <  2:259682/150‑1 (MQ=255)
                                                                                                           gCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCg                                           >  2:999348/1‑150 (MQ=255)
                                                                                                             tGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGCTCGGCTCGGCGAAGGATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCCGCCACCAGGAAGTCGAGGAACtgacac                                         >  2:399027/1‑145 (MQ=255)
                                                                                                               ccACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTc                                       >  2:604537/1‑150 (MQ=255)
                                                                                                                 acaTGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAg                                                   <  1:756598/136‑1 (MQ=255)
                                                                                                                 acaTGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAg                                                   >  2:756598/1‑136 (MQ=255)
                                                                                                                    tGCGCTGGAAGTAGCCCTCGGGCGGAATGATCGCCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCa                                  <  1:950468/150‑1 (MQ=255)
                                                                                                                    tGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCCAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCCCCCAGCGCGAGGCTCTTGTCCTCGAACTCCGCCCCCAGGAAGTCGAGGAACTCCGCCTCGTccc                                  >  1:675710/1‑149 (MQ=255)
                                                                                                                       gcTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATgt                               >  2:1079679/1‑150 (MQ=255)
                                                                                                                               tAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGgc                       <  1:299675/150‑1 (MQ=255)
                                                                                                                                 gggcTCGGGCGGGGTGATCATTCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCg                                                                    <  2:1346517/100‑1 (MQ=38)
                                                                                                                                 gggcTCGGGCGGGGTGATCATTCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCg                                                                    >  1:1346517/4‑103 (MQ=37)
                                                                                                                                 gCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACATTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGGCGAGGAACTCCGCCCCGTCCATGTCGTCTGgcgg                     >  2:436719/1‑149 (MQ=255)
                                                                                                                                   ccTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGgcgcgc                   >  1:1567935/1‑150 (MQ=255)
                                                                                                                                     tCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCgg                 <  1:999348/150‑1 (MQ=255)
                                                                                                                                        ggCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTAc              <  2:1597760/150‑1 (MQ=255)
                                                                                                                                           ggAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCAAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACggg           <  1:1120194/150‑1 (MQ=255)
                                                                                                                                           ggAATGATCACCCCGCCCGAGCCCATGATAGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCTATCGAGAGGATATTGTCTTCGAACTAGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCGCTGTCGTCTTGCGCGCGGTACggg           >  1:317655/1‑150 (MQ=255)
                                                                                                                                               tGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggttgg       <  1:23266/150‑1 (MQ=255)
                                                                                                                                                 aTCACCCCGCCCGAGCCCATGATGGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggttgggg     >  2:946118/1‑150 (MQ=255)
                                                                                                                                                 aTCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggttgggg     >  2:1116876/1‑150 (MQ=255)
                                                                                                                                                 aTCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggttgggg     <  2:1436434/150‑1 (MQ=255)
                                                                                                                                                  tCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCCGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggttggggt    >  2:336685/1‑150 (MQ=255)
                                                                                                                                                    aCCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggttgggggtg  >  1:284011/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCTGGGGTGACGCCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCAC‑CACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTG  >  NC_002947/4240877‑4241173

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: ‑