Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 274,082 A→C 5.9% I163S (ATC→AGC)  PP_0221 ← methionine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947274,0820AC 5.9% 216.6 / 2.4 67I163S (ATC→AGC) PP_0221methionine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (36/27);  new base C (2/2);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ACAGCAACGCCGAGGTAGCATCGAAGGTGTTGGCCAGGCGGATGTAGGCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTG  >  NC_002947/273939‑274228
                                                                                                                                               |                                                                                                                                                  
aCAGCAACGCCGAGGTAGCATCGAAGGTGTTGGCCAGGCGGATGTAGGCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGt                                                                                                                                                 >  2:149915/1‑147 (MQ=255)
 cAGCAACGCCGAGGTAGCATCGAAGGTGTTGGCCAGGCGGATGTAGGCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAg                                                                                                                                              <  4:17993/149‑1 (MQ=255)
    cAACGCCGAGGTAGCATCGAAGGTGTTGGCCAGGCGGATGTAGGCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCt                                                                                                                                           >  2:60981/1‑149 (MQ=255)
       cGCCGAGGTAGCATCGAAGGTGTTGGCCAGGCGGATGTAGGCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCgg                                                                                                                                        <  1:149915/149‑1 (MQ=255)
            aGGTAGCATCGAAGGTGTTGGCCAGGCGGATGTAGGCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGcc                                                                                                                                   >  2:34643/1‑149 (MQ=255)
              gTAGCATCGAAGGTGTTGGCCAGGCGGATGTAGGCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCtt                                                                                                                                 >  6:236798/1‑149 (MQ=255)
              gTAGCATCGAAGGTGTTGGCCAGGCGGATGTAGGCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCtt                                                                                                                                 <  7:130541/149‑1 (MQ=255)
                  cATCGAAGGTGTTGGCCAGGCGGATGTAGGCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAATGATCTTCAGCTTTTTAGGGTTGGCGAGGCTGACATCCTCgc                                                                                                                                        >  8:2615/1‑137 (MQ=255)
                  cATCGAAGGTGTTGGCCAGGCGGATGTAGGCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGCTGACATCCTCgc                                                                                                                                        >  8:6308/1‑137 (MQ=255)
                  cATCGAAGGTGTTGGCCAGGCGGATGTAGGCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGCTGACATCCTCgc                                                                                                                                        <  7:6308/138‑2 (MQ=255)
                  cATCGAAGGTGTTGGCCAGGCGGATGTAGGCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGCTGACATCCTCgc                                                                                                                                        <  7:2615/138‑2 (MQ=255)
                         ggTGTTGGCCAGGCGGATGTAGGCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGc                                                                                                                      <  3:98432/149‑1 (MQ=255)
                                gCCAGGCGGATGTAGGCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCt                                                                                                                <  4:246760/148‑1 (MQ=255)
                                  ctggCGGATGGAGGCCGGGTAACCGTGCACGAGGTCGGCGTCGTGGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCGCGGTGGCCTTGTAGGCCACGCCTGGCTTGa                                                                                                             <  2:132798/147‑1 (MQ=255)
                                           gTAGGCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCATGTAGCCCACGCCGGGCTTGAGGGTGATc                                                                                                     >  2:252194/1‑148 (MQ=255)
                                              ggCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGAGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAAcc                                                                                                 >  2:124595/1‑149 (MQ=255)
                                               gCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAAccc                                                                                                >  6:222204/1‑149 (MQ=255)
                                               gCCGGGT‑ACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCg                                                                                               >  6:82028/1‑149 (MQ=255)
                                                        cccTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGc                                                                                       >  3:165350/1‑149 (MQ=255)
                                                        cccTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGc                                                                                       >  3:168761/1‑149 (MQ=255)
                                                         ccTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTgg                                                                                                                                     >  2:242616/1‑102 (MQ=255)
                                                         ccTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTgg                                                                                                                                     <  1:242616/102‑1 (MQ=255)
                                                                   gTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCAGGCCTTGGCcagcagcggcc                                                                            >  1:18006/1‑148 (MQ=255)
                                                                     cGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAgg                                                                          >  7:62546/1‑149 (MQ=255)
                                                                      ggCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAg                                                                                                                             >  3:53373/1‑97 (MQ=255)
                                                                      ggCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAg                                                                                                                             <  4:53373/97‑1 (MQ=255)
                                                                          tcgtcgtAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCTGCTGCTGCTGGCCaan                                                                     <  4:168761/149‑3 (MQ=255)
                                                                          tcgtcgtAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCTGCTGCTGCTGGCCaag                                                                     <  4:165350/149‑3 (MQ=255)
                                                                                tAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAgg                                                                                                           <  2:110755/105‑1 (MQ=255)
                                                                                tAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAgg                                                                                                           >  1:110755/1‑105 (MQ=255)
                                                                                    cgcgTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTg                                                           >  2:173489/1‑149 (MQ=255)
                                                                                     gcgTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCcagcagcagc                                                                             <  7:253066/130‑1 (MQ=255)
                                                                                     gcgTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCcagcagcagc                                                                             >  8:253066/1‑130 (MQ=255)
                                                                                        tACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCg                                                       >  8:242886/1‑149 (MQ=255)
                                                                                           cAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAgg                                                                                                           >  4:210524/1‑94 (MQ=255)
                                                                                           cAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAgg                                                                                                           <  3:210524/94‑1 (MQ=255)
                                                                                           cAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGt                                                    >  4:55004/1‑149 (MQ=255)
                                                                                                 cACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTgg                                              <  6:241416/149‑1 (MQ=255)
                                                                                                  aCCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGc                                             >  3:144009/1‑149 (MQ=255)
                                                                                                  aCCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGc                                             >  6:31273/1‑149 (MQ=255)
                                                                                                   ccGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGc                                             >  5:72856/1‑148 (MQ=255)
                                                                                                      ccTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATg                                         >  1:131030/1‑149 (MQ=255)
                                                                                                       cTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATgg                                        <  1:173489/149‑1 (MQ=255)
                                                                                                        tCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGACCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGc                                       >  2:224267/1‑149 (MQ=255)
                                                                                                          gACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGc                                       <  2:131030/147‑1 (MQ=255)
                                                                                                          gACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCa                                     <  3:55004/149‑1 (MQ=255)
                                                                                                             cTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCgg                                                      <  8:106384/129‑1 (MQ=255)
                                                                                                              tGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCt                                 >  3:262503/1‑149 (MQ=255)
                                                                                                               gAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTg                                >  1:138240/1‑149 (MQ=255)
                                                                                                               gAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTg                                >  5:122107/1‑149 (MQ=255)
                                                                                                                      cTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCtt                                                                                          >  3:37429/1‑84 (MQ=255)
                                                                                                                      cTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCtt                                                                                          <  4:37429/84‑1 (MQ=255)
                                                                                                                        tCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGACGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCtt                         >  7:173817/1‑147 (MQ=255)
                                                                                                                         cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCa                       <  5:236798/148‑1 (MQ=255)
                                                                                                                         cAGCTTTTTAGAGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAg                      >  6:52548/1‑149 (MQ=255)
                                                                                                                          aGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAgg                     <  4:144009/149‑1 (MQ=255)
                                                                                                                          aGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAgg                     <  8:62546/149‑1 (MQ=255)
                                                                                                                                        ggCGAGGATGTCAGCCTCGGTGGCCATGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAg        <  1:252194/148‑1 (MQ=255)
                                                                                                                                         gCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCt      >  8:142507/1‑149 (MQ=255)
                                                                                                                                           gAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAgg                     <  5:174451/132‑1 (MQ=255)
                                                                                                                                           gAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAgg                     >  6:174451/1‑132 (MQ=255)
                                                                                                                                           gAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCttt    >  7:145049/1‑149 (MQ=255)
                                                                                                                                            aGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCtttt   <  8:173817/149‑1 (MQ=255)
                                                                                                                                            aGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCtttt   <  6:122107/149‑1 (MQ=255)
                                                                                                                                             ggATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTg  >  2:247546/1‑149 (MQ=255)
                                                                                                                                              gATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCtttt   <  5:196958/147‑1 (MQ=255)
                                                                                                                                              gATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCtttt   >  6:196958/1‑147 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
ACAGCAACGCCGAGGTAGCATCGAAGGTGTTGGCCAGGCGGATGTAGGCCGGGTAACCCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTG  >  NC_002947/273939‑274228

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