Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,524,933 C→G 5.5% R122G (CGC→GGC) ‡ 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,9330CG 5.5% 240.1 / 3.4 73R122G (CGC→GGC) ‡PP_4013MutT/nudix family protein
Reads supporting (aligned to +/- strand):  ref base C (29/40);  new base G (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.64e-01

CACCGCCCCCAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGACCAGCCCA  >  NC_002947/4524785‑4525081
                                                                                                                                                     |                                                                                                                                                    
cacCGCCCCCAGCTACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGgc                                                                                                                                                      >  3:180249/1‑149 (MQ=255)
cacCGCCCCCAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGCCCCCACTCGCTGgc                                                                                                                                                      >  8:315693/1‑149 (MQ=255)
 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                                                                                                                                           <  1:102311/149‑1 (MQ=255)
           gCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTg                                                                                                                                           <  2:265446/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                    >  6:95191/1‑149 (MQ=255)
                                                                                                                                  gCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCg                    >  2:145212/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)
                                                                                                                                                     |                                                                                                                                                    
CACCGCCCCCAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGACCAGCCCA  >  NC_002947/4524785‑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: