Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,524,934 G→C 5.5% R122P (CGC→CCC) ‡ PP_4013 → MutT/nudix family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,524,9340GC 5.5% 237.3 / 3.1 73R122P (CGC→CCC) ‡PP_4013MutT/nudix family protein
Reads supporting (aligned to +/- strand):  ref base G (29/40);  new base C (2/2);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

ACCGCCCCCAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGACCAGCCCA  >  NC_002947/4524786‑4525081
                                                                                                                                                     |                                                                                                                                                   
acCGCCCCCAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGgcg                                                                                                                                                     >  3:25438/1‑149 (MQ=255)
       ccAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAg                                                                                                                                              <  2:190799/149‑1 (MQ=255)
        cAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGc                                                                                                                                             <  1:119391/149‑1 (MQ=255)
         aGCTACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCt                                                                                                                                            <  4:180249/149‑1 (MQ=255)
         aGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGCGCCAGCGAgtg                                                                                                                                            <  1:165438/149‑3 (MQ=255)
         aGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGCGCCAGCGAgtg                                                                                                                                            <  5:177929/149‑3 (MQ=255)
          gCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTg                                                                                                                                           <  2:265446/149‑1 (MQ=255)
          gCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTg                                                                                                                                           <  1:102311/149‑1 (MQ=255)
                gCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCcg                                                                                                                                     >  2:64007/1‑149 (MQ=255)
                 cGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATCTCGTCCGCGCCGTGTGGCTGACCCGCGACGAAGTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGcgccgc                                                                                                                                    <  5:260618/149‑1 (MQ=255)
                     aCCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCTGTGCCGTGTGGCTGACGCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCt                                                                                                                                 <  5:274459/148‑1 (MQ=255)
                     aCCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCt                                                                                                                                 <  3:84681/148‑1 (MQ=255)
                      ccTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGc                                                                                                                               <  8:24167/149‑1 (MQ=255)
                       cTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGcc                                                                                                                              >  4:346651/1‑149 (MQ=255)
                         aCCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCt                                                                                                                                 >  1:380178/1‑144 (MQ=255)
                         aCCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCt                                                                                                                                 <  2:380178/144‑1 (MQ=255)
                         aCCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTg                                                                                                                            >  3:172340/1‑149 (MQ=255)
                           cAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGa                                                                                                                          <  1:362651/149‑1 (MQ=255)
                                    tGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCgcgc                                                                                                                                  >  2:163656/1‑132 (MQ=255)
                                    tGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCgcgc                                                                                                                                  <  1:163656/132‑1 (MQ=255)
                                     gCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCt                                                                                                                 >  8:144812/1‑148 (MQ=255)
                                       ttCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCa                                                                                                               >  6:245180/1‑148 (MQ=255)
                                         cgccgcCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCaa                                                                                                              <  4:25438/147‑1 (MQ=255)
                                           acgcccgccCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCTACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGc                                                                                                          <  8:184077/148‑1 (MQ=255)
                                            cgcccgccCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGcc                                                                                                         <  7:123570/149‑1 (MQ=255)
                                             gcccgcccCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCAGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGCCTACCTCAAAGGccc                                                                                                        >  8:260520/1‑149 (MQ=255)
                                             gcccgcccCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGACTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGccc                                                                                                        >  3:75358/1‑149 (MQ=255)
                                              cccgcccCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGGGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCg                                                                                                       <  3:15853/149‑1 (MQ=255)
                                               ccgcccCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGc                                                                                                      >  7:320395/1‑149 (MQ=255)
                                                   cccGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTTGATGACTACCTCAAAGGCCCGCTGca                                                                                                  <  5:45271/149‑1 (MQ=255)
                                                   cccGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTAc                                                                                                                   >  3:276590/1‑132 (MQ=255)
                                                   cccGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTAc                                                                                                                   <  4:276590/132‑1 (MQ=255)
                                                     cGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGTCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGca                                                                                                  >  3:123062/1‑147 (MQ=255)
                                                       agcgcCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGc                                                                                              <  6:172411/148‑1 (MQ=255)
                                                         cgcCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGTCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGcc                                                                                             <  4:123062/148‑1 (MQ=255)
                                                         cgcCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGcc                                                                                             <  8:41394/148‑1 (MQ=255)
                                                          gcCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTc                                                                                           <  6:190218/149‑1 (MQ=255)
                                                             aTCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGAc                                                                                        <  2:382908/149‑1 (MQ=255)
                                                                 cGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTg                                                                                     >  8:181469/1‑148 (MQ=255)
                                                                      aCCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTAcgc                                                                               <  7:144812/149‑1 (MQ=255)
                                                                          ggCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTg                                                                                     >  5:270642/1‑139 (MQ=255)
                                                                          ggCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTg                                                                                     <  6:270642/139‑1 (MQ=255)
                                                                               tGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGAc                                                                            <  4:146613/143‑1 (MQ=255)
                                                                               tGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGAc                                                                            >  3:146613/1‑143 (MQ=255)
                                                                                  aCAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACgcggc                                                                   >  3:331988/1‑149 (MQ=255)
                                                                                    aGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGgcgc                                                                 >  1:183300/1‑149 (MQ=255)
                                                                                      cGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTCGCCGACCCCCCTCGCTGGCGCAGCCAGCTCGTGCCGCGCTGCCAGGATGACTACCTCAAAGGCCCGCTGCACATCCTCGACCTGCTACGCGACTTCCGCGGagccg                                                               >  8:115125/1‑149 (MQ=255)
                                                                                       gATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGcgccgc                                                              <  5:245180/149‑1 (MQ=255)
                                                                                           tCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCggg                                                           <  5:326075/148‑1 (MQ=255)
                                                                                             gTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGa                                                                             <  1:103700/128‑1 (MQ=255)
                                                                                             gTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGa                                                                             >  2:103700/1‑128 (MQ=255)
                                                                                             gTCCGCGCCGTGTGACTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGttt                                                        <  4:75358/149‑1 (MQ=255)
                                                                                               ccGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGa                                                      <  4:172340/149‑1 (MQ=255)
                                                                                                cgcgCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGAt                                                     <  3:346651/149‑1 (MQ=255)
                                                                                                         tgGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGc                                            <  4:166527/149‑1 (MQ=255)
                                                                                                          ggCTGACCCGCGACGAACCGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCg                                           >  3:259840/1‑149 (MQ=255)
                                                                                                           gCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCGTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCg                                           >  1:69339/1‑149 (MQ=255)
                                                                                                             tGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGttt                                        >  3:350074/1‑149 (MQ=255)
                                                                                                               aCCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGttttt                                      <  7:181469/149‑1 (MQ=255)
                                                                                                                      tcgaACTGCGGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCACGCTGCAGAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCtt                               <  3:272453/148‑1 (MQ=255)
                                                                                                                        gaACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATTACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGa                             <  8:335370/149‑1 (MQ=255)
                                                                                                                        gaACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCttta                             <  4:331988/149‑3 (MQ=255)
                                                                                                                        gaACTGCTGGCCGACCCCACTCGCTCGCGCAGCGAGCTGGTGACGCGCTGCCTGGAAGACTACCTCAAAGGCACGCTCCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGAAAGAATCGGCGATTTTCCCCCttta                             >  2:232881/1‑147 (MQ=255)
                                                                                                                           ctgcagGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTCGATGACTACCTCAAAGGCCCGCAGCACAGCCTCGACCTGCTACGCGACTCACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATa                           >  2:388958/1‑148 (MQ=255)
                                                                                                                                 gCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCg                    >  2:145212/1‑149 (MQ=255)
                                                                                                                                 gCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCg                    >  6:95191/1‑149 (MQ=255)
                                                                                                                                   cGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGt                  <  7:315693/149‑1 (MQ=255)
                                                                                                                                      accagcTCGCTGCGCCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTAcc               >  8:286359/6‑149 (MQ=255)
                                                                                                                                       cccACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTAcc               >  5:359066/1‑148 (MQ=255)
                                                                                                                                       ccagcTCGCTGCGCCAGCGAGTTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTAcc               >  7:332263/5‑148 (MQ=255)
                                                                                                                                         cACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCAt            >  2:229822/1‑149 (MQ=255)
                                                                                                                                         cACTCGCTGGCGCAGCGAGATGGTGCAGCGCTGCCTGGATGACTACCTCAAAGGCCAGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCCGGTTTGATAGAATCGGAGTTTTTCCCCCTTGATACACACAGGTACCCAt            >  4:236392/1‑149 (MQ=255)
                                                                                                                                           cTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGa          <  4:259840/149‑1 (MQ=255)
                                                                                                                                              gCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGAc         <  3:236392/147‑1 (MQ=255)
                                                                                                                                                   gcgcAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGAccagcc    >  1:73154/1‑148 (MQ=255)
                                                                                                                                                   gcgcAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGAccagcc    <  1:64007/148‑1 (MQ=255)
                                                                                                                                                    cgcAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGAccagccca  <  1:229822/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACCGCCCCCAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGACCAGCCCA  >  NC_002947/4524786‑4525081

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: