Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 956,584 A→T 7.8% intergenic (+17/‑97) alaC → / → PP_0818 aminotransferase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947956,5840AT 7.8% 182.8 / 5.1 64intergenic (+17/‑97)alaC/PP_0818aminotransferase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/30);  new base T (3/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

GGGATCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTCGAC  >  NC_002947/956439‑956731
                                                                                                                                                 |                                                                                                                                                   
gggATCGGCTTTGGTGACTATGGGGACGACCATGGGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAagc                                                                                                                                                  <  2:290916/149‑1 (MQ=255)
  gATCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAagcag                                                                                                                                                <  4:1130/149‑1 (MQ=255)
    tCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGcc                                                                                                                                              >  7:108145/1‑149 (MQ=255)
       gCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGcccc                                                                                                                                            >  4:206872/1‑148 (MQ=255)
         tttGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAAc                                                                                                                                         >  7:119251/1‑149 (MQ=255)
           ttgTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCTGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCTTCGCAGCCCCatt                                                                                                                                         <  5:161684/145‑3 (MQ=255)
             gTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGcc                                                                                                                                     >  6:211040/1‑149 (MQ=255)
             gTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGc                                                                                                                                      <  6:372158/148‑1 (MQ=255)
             gTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGACAACCCGCACCGTATTCGCCAGGCCAGCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGc                                                                                                                                      <  8:72583/148‑1 (MQ=255)
              tGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCct                                                                                                                                    >  4:319729/1‑149 (MQ=255)
                 cTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCtt                                                                                                                                 <  3:206872/149‑1 (MQ=255)
                       ggACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGt                                                                                                                             >  4:103876/1‑147 (MQ=255)
                       ggACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGa                                                                                                                           <  1:297464/149‑1 (MQ=255)
                             ccATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTc                                                                                                                     <  5:211040/149‑1 (MQ=255)
                             ccATGTGCGCTTTGCGTTGATCGAGAACCAGGACCATATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTc                                                                                                                     >  2:182378/1‑149 (MQ=255)
                               aTGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCg                                                                                                                   >  5:161268/1‑149 (MQ=255)
                                tgtgCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTccac                                                                                                                  <  1:99370/149‑3 (MQ=255)
                                    cgcTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCgcctgc                                                                                                              <  3:319729/149‑1 (MQ=255)
                                          gCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCTAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCt                                                                                                         <  2:215442/148‑1 (MQ=255)
                                          gCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTg                                                                                                        >  1:170156/1‑149 (MQ=255)
                                           cGTTGACCGACAACCAGCACCGTTTTCGCCAGGCCATCCGCGGTAGCCGGCAAATGTTCCGGGCTGTTGGTTGGGCGCGCAAGTACGCCAAATGGGGCGGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTgg                                                                                                       <  6:6187/149‑1 (MQ=255)
                                            gTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGt                                                                                                      >  7:189201/1‑149 (MQ=255)
                                            gTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGt                                                                                                      >  2:119569/1‑149 (MQ=255)
                                             ttGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGtt                                                                                                     <  8:189201/149‑1 (MQ=255)
                                                    agaACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAACGGGGCTGCTTCGCAGCCCCattgg                                                                                                                                       <  3:169775/108‑5 (MQ=255)
                                                    agaACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAACGGGGCTGCTTCGCAGCCCCattgg                                                                                                                                       >  4:169775/1‑104 (MQ=255)
                                                         cAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAagcagc                                                                                                                                               >  3:188283/1‑95 (MQ=255)
                                                         cAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAagcagc                                                                                                                                               <  4:188283/95‑1 (MQ=255)
                                                               cGTCTTCGCCAGGCCAGCGGCGGGATCCGGCAAATGTGCCGGGCTGAGGGTTTGGCTCGCAAGTAAGCCCAATGGGGCGGCGAACCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTc                                                                                   <  4:36155/149‑1 (MQ=255)
                                                               cGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTAc                                                                                                   >  3:190968/1‑133 (MQ=255)
                                                               cGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTc                                                                                   >  4:71051/1‑149 (MQ=255)
                                                                 tATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAg                                                                                 <  5:337792/149‑1 (MQ=255)
                                                                      gacAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTAGTAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACt                                                                             >  2:251384/3‑148 (MQ=255)
                                                                      gCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCAcg                                                                             >  3:22365/1‑147 (MQ=255)
                                                                      gCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTc                                                                            <  6:161268/149‑1 (MQ=255)
                                                                         aGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCAtt                                                                         <  7:121744/149‑1 (MQ=255)
                                                                          ggCCATCCGCGGGATCCGGCAAATGTTCCGGGCTTATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCAttt                                                                        <  3:71051/149‑1 (MQ=255)
                                                                                      gATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCaaa                                                               >  2:295115/1‑146 (MQ=255)
                                                                                          cGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATcc                                                        >  7:167495/1‑149 (MQ=255)
                                                                                           ggCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATccc                                                       <  8:167495/149‑1 (MQ=255)
                                                                                           ggCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATccc                                                       >  4:155329/1‑149 (MQ=255)
                                                                                           ggCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGACTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATccc                                                       >  2:223676/1‑149 (MQ=255)
                                                                                                    tCCGGGCTGTTGTTTTGGCGCGAAAGTAAGCCAATTGGGGCTGCAAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAgt                                               <  5:57897/148‑1 (MQ=255)
                                                                                                         gCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTTCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTc                                         <  1:81634/149‑1 (MQ=255)
                                                                                                           tGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCAc                                       >  7:126552/1‑149 (MQ=255)
                                                                                                             aTGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCg                                     >  1:216398/1‑149 (MQ=255)
                                                                                                                gctaggCTCGCAAGATAGCCCATTGGGCCTGCGACGCAGCCCCAACTGCCTCTTGGTTGAGGTTTCCGCCTGCCTTCTGGTTACCCGTACGCATACCGTCAGCACTCATTTCTCACCAAAGAGATCCCGCATGAGTGTCTTCACCGCCt                                  <  8:103376/145‑1 (MQ=255)
                                                                                                                 tttGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTa                                 >  5:318705/1‑149 (MQ=255)
                                                                                                                        cGCACGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGCGGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGc                          <  8:168218/149‑1 (MQ=255)
                                                                                                                         gttaGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCg                         <  1:251384/146‑1 (MQ=255)
                                                                                                                         gCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCg                         <  1:295115/149‑1 (MQ=255)
                                                                                                                             gTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGa                                                 <  5:345227/121‑1 (MQ=255)
                                                                                                                             gTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGa                                                 >  6:345227/1‑121 (MQ=255)
                                                                                                                              tAAGCACAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACt                    <  8:108145/149‑1 (MQ=255)
                                                                                                                                  cccAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCGTCACCGCCTACTTCTGCGGCACTGGCt                >  7:136164/1‑149 (MQ=255)
                                                                                                                                      ttGGGGCTGCTTCGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCGTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCa            >  1:50294/2‑149 (MQ=255)
                                                                                                                                         gggCTGCTTCGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCg         >  2:54225/1‑149 (MQ=255)
                                                                                                                                           gCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGc                 <  1:209688/139‑1 (MQ=255)
                                                                                                                                           gCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGc                 >  2:209688/1‑139 (MQ=255)
                                                                                                                                             tGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTc     >  8:165644/1‑149 (MQ=255)
                                                                                                                                               ggAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTcga   <  8:119251/148‑1 (MQ=255)
                                                                                                                                               cGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTcga   <  1:182378/149‑1 (MQ=255)
                                                                                                                                               cGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTcga   >  8:272311/1‑149 (MQ=255)
                                                                                                                                                gAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTcgac  >  1:313796/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GGGATCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTCGAC  >  NC_002947/956439‑956731

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