Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,910,826 G→T 8.3% A46S (GCG→TCG)  PP_3451 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,910,8260GT 8.3% 182.8 / 5.7 60A46S (GCG→TCG) PP_3451hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/22);  new base T (2/3);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

TTCTGATGCGTACGCTCATCGCCTTCTGCCTGCTGGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGAGGGCTTGCAGCTGGTGCTGGGGCTGTACATGGACCCG  >  NC_002947/3910686‑3910972
                                                                                                                                            |                                                                                                                                                  
ttCTGATGCGTACGCTCATCGCCTTCTGCCTGCTGGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGccc                                                                                                                                              >  3:321940/1‑147 (MQ=255)
     aTGCGTACGCTCATCGCCTTCTGCCTGCTGGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTg                                                                                                                                        <  6:49861/148‑1 (MQ=255)
        cGTACGCTCATCGCCTTCTGCCTGCTGGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGTCAAGCCACCTTGGGCCGCGGTGGctt                                                                                                                                    <  4:202066/149‑3 (MQ=255)
         gTACGCTCATCGCCTTCTGCCTGCTGGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCTTGGGCCGCGGTGGcttg                                                                                                                                   <  5:39045/149‑4 (MQ=255)
         gTACGCTCATCGCCTTCTGCCTGCTGGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCaa                                                                                                                                    >  2:12588/1‑148 (MQ=255)
            cGCTCATCGCCTTCTGCCTGCTGGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAAc                                                                                                                                 <  4:321940/148‑1 (MQ=255)
              cTCATCGCCTTCTGCCTGCTGGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAAcc                                                                                                                                >  4:366545/1‑147 (MQ=255)
              cTCATCGCCTTCTGCCTGCTGGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAAcc                                                                                                                                <  3:366545/147‑1 (MQ=255)
               tCATCGCCTTCTGCCTGCTGGCCCTGTGCCTACCGGCCACCGCTGCCCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCTGCAACTGTGGCAAGCCACCTCGGCCCAAGGTGGCAAAcccag                                                                                                                              >  4:273484/1‑146 (MQ=255)
                 aTCGCCTTCTGCCTGCTGGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTg                                                                                                                             <  1:283745/147‑1 (MQ=255)
                      cTTCTGCCTGCTGGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTg                                                                                                                       >  3:62210/1‑148 (MQ=255)
                          tgcctgcTGGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTg                                                                                                                       <  4:212064/144‑1 (MQ=255)
                          tgcctgcTGGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTg                                                                                                                       >  3:212064/1‑144 (MQ=255)
                           gcctgcTGGCCCTGTGCCTACCGGCCCCCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCCACGGTGGCAAGCCACCGCGGCCCAAGGTGGCACAACCCTGAGCGGGCAATgg                                                                                                                 >  8:336996/1‑149 (MQ=255)
                            cctgcTGGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATg                                                                                                                  <  1:12588/147‑1 (MQ=255)
                             ctgctgGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCTTGGGCCGCGGTGGc                                                                                                                                      <  3:361275/126‑1 (MQ=255)
                             ctgctgGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCTTGGGCCGCGGTGGc                                                                                                                                      >  4:361275/1‑126 (MQ=255)
                                     ccTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATg                                                                                                       >  4:226634/1‑149 (MQ=255)
                                      cTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATgg                                                                                                      >  4:46239/1‑149 (MQ=255)
                                            cTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAg                                                                                                >  2:11065/1‑149 (MQ=255)
                                             tACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTg                                                                                                                                        >  3:345302/1‑108 (MQ=255)
                                             tACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTg                                                                                                                                        <  4:345302/108‑1 (MQ=255)
                                             tACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGc                                                                                               >  5:94055/1‑149 (MQ=255)
                                               ccGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGa                                                                                             <  2:262958/149‑1 (MQ=255)
                                                     aCCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTcgg                                                                                       <  6:94055/149‑1 (MQ=255)
                                                       cGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTcggcg                                                                                     >  5:72797/1‑149 (MQ=255)
                                                       cGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTcggcg                                                                                     >  8:127127/1‑149 (MQ=255)
                                                          tGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTg                                                                                  >  5:119804/1‑149 (MQ=255)
                                                            aCCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGcc                                                                                >  8:271338/1‑149 (MQ=255)
                                                            aCCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGcc                                                                                >  6:304644/1‑149 (MQ=255)
                                                                    cTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGggctg                                                                        >  8:229828/1‑149 (MQ=255)
                                                                    cTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGCGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGggctg                                                                        >  2:229471/1‑149 (MQ=255)
                                                                         cTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTcggc                                                                                      <  6:244384/130‑1 (MQ=255)
                                                                          tACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCtggctggc                                                                  >  8:90068/1‑149 (MQ=255)
                                                                           aCCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCtggctggc                                                                  >  5:265636/1‑148 (MQ=255)
                                                                             cAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCaa                                                               <  7:90068/149‑1 (MQ=255)
                                                                             cAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCa                                                                >  2:80718/1‑148 (MQ=255)
                                                                               aCCCCTCCACCGCGATGCCAGCGTCACACCCGCCCCGTGGCAGCACCTGTGGCAAGCCACCGCGGCCCAAGGGGGCAAAACCCTGATCGGGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACa                                                             <  3:283845/149‑1 (MQ=255)
                                                                                  ccTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGcc                                                          >  7:113456/1‑149 (MQ=255)
                                                                                   cTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGc                                                           >  4:376940/1‑147 (MQ=255)
                                                                                    tCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTg                                                        >  3:380757/1‑149 (MQ=255)
                                                                                      aaCCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTgc                                                       <  4:62210/148‑1 (MQ=255)
                                                                                         cgcgATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCAc                                                   >  8:356602/1‑149 (MQ=255)
                                                                                          gcgATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCAcc                                                  >  7:325778/1‑149 (MQ=255)
                                                                                              tGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACcgcgc                                              <  5:74583/149‑1 (MQ=255)
                                                                                                     gTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGCGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGa                                       <  1:80718/149‑1 (MQ=255)
                                                                                                             cGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGAGGGCTTGc                               <  7:336996/149‑1 (MQ=255)
                                                                                                                   gTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGAGGGCTTGCAGCTgg                         <  3:226634/149‑1 (MQ=255)
                                                                                                                      gcagcaACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGAGGGCTTGCAGCTGGTGc                      >  6:27476/1‑149 (MQ=255)
                                                                                                                      gcagcaACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGAGGGCTTGCAGCTGGTGc                      <  3:46239/149‑1 (MQ=255)
                                                                                                                         gcaACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGAGGGCTTGCAGCTGGTGCTgg                   >  7:283092/1‑149 (MQ=255)
                                                                                                                         gcaACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGAGGGCTTGCAGCTGGTGCTgg                   <  8:283092/149‑1 (MQ=255)
                                                                                                                          caACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGAGGGCTTGCAGCTGGTGCTggg                  <  7:271338/149‑1 (MQ=255)
                                                                                                                           aaCTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGAGGGCTTGCAGCTGGTGCTgggg                 <  7:356602/149‑1 (MQ=255)
                                                                                                                           aaCTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGAGGGCTTGCAGCTGGTGCTgggg                 >  2:159207/1‑149 (MQ=255)
                                                                                                                            aCTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGAGGGCTTGCAGCTGGTGCTGGGGc                <  7:229828/149‑1 (MQ=255)
                                                                                                                                tgGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGAGGGCTTGCAGCTGGTGCTGGGGCTGTa            >  1:350993/1‑149 (MQ=255)
                                                                                                                                  gCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGAGGGCTTGCAGCTGGTGCTGGGGCTGTACa          >  4:99092/1‑149 (MQ=255)
                                                                                                                                   cAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGAGGGCTTGCAGCTGGTGCTGGGGCTGTACAt         >  5:346424/1‑149 (MQ=255)
                                                                                                                                          ccGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGAGGGCTTGCAGCTGGTGCTGGGGCTGTACATGGACCCg  >  1:159059/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
TTCTGATGCGTACGCTCATCGCCTTCTGCCTGCTGGCCCTGTGCCTACCGGCCACCGCTGACCAGCGCCTGTTCTACCAACCCCTCAACCGCGATGCCAGCGTCACACCCGCCCAGTGGCAGCAACTGTGGCAAGCCACCGCGGCCCAAGGTGGCAAAACCCTGATCGTGCAATGGAGTGCCTATGGCGACAGCGATTTCGGCGGTGCCCAGGGCTGGCTGGCCAACAGCCTGCGCACCGCGCGTGCCGAGGGCTTGCAGCTGGTGCTGGGGCTGTACATGGACCCG  >  NC_002947/3910686‑3910972

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