Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,886,588 A→G 7.2% I419V (ATC→GTC)  rmuC → DNA recombination limiting inversions protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,886,5880AG 7.2% 199.5 / 4.7 69I419V (ATC→GTC) rmuCDNA recombination limiting inversions protein
Reads supporting (aligned to +/- strand):  ref base A (30/34);  new base G (3/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

CAGATCGTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTT  >  NC_002947/1886447‑1886736
                                                                                                                                             |                                                                                                                                                     
cAGATCGTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGa                                                                                                                                                <  7:205129/149‑1 (MQ=255)
 aGATCGTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGAc                                                                                                                                               <  4:143171/149‑1 (MQ=255)
    tCGTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCt                                                                                                                                             <  4:187370/148‑1 (MQ=255)
    tCGTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTg                                                                                                                                           >  7:224507/1‑149 (MQ=255)
     cGTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTg                                                                                                                                           >  8:220705/1‑148 (MQ=255)
      gTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTgg                                                                                                                                          <  4:93263/148‑1 (MQ=255)
      gTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGa                                                                                                                                         >  6:65226/1‑149 (MQ=255)
      gTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGa                                                                                                                                         <  2:110445/149‑1 (MQ=255)
       tGATTGTCAGCCCTACCACGTTGCTGGCCACCTTTCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGAc                                                                                                                                        >  6:148878/1‑149 (MQ=255)
             tCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACt                                                                                                                                   <  6:229355/148‑1 (MQ=255)
             tCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTg                                                                                                                                  >  2:60666/1‑149 (MQ=255)
               aGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTggg                                                                                                                                <  8:224507/149‑1 (MQ=255)
               aGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTggg                                                                                                                                >  6:193022/1‑149 (MQ=255)
               aGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTgg                                                                                                                                 <  5:174845/148‑1 (MQ=255)
                gCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGc                                                                                                                               >  4:760/1‑149 (MQ=255)
                  ccTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAg                                                                                                                             >  2:72975/1‑149 (MQ=255)
                       cACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATc                                                                                                                                                     <  1:56573/121‑1 (MQ=255)
                       cACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATc                                                                                                                                                     >  2:56573/1‑121 (MQ=255)
                        aCGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCtt                                                                                                                       <  4:83722/149‑1 (MQ=255)
                          gTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGc                                                                                                                         >  3:150679/1‑145 (MQ=255)
                          gTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGc                                                                                                                         <  4:150679/145‑1 (MQ=255)
                          gTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTgc                                                                                                                     <  2:280323/149‑1 (MQ=255)
                               tGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGcagcag                                                                                                                >  8:324014/1‑149 (MQ=255)
                               tGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGcagcag                                                                                                                >  2:278960/1‑149 (MQ=255)
                                  ccACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTg                                                                                                             <  5:173798/149‑1 (MQ=255)
                                  ccACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTg                                                                                                             <  1:2244/149‑1 (MQ=255)
                                       ttGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCtt                                                                                                                       <  2:309209/134‑1 (MQ=255)
                                       ttGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCtt                                                                                                                       >  1:309209/1‑134 (MQ=255)
                                            ggTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCtt                                                                                                   <  4:101380/149‑1 (MQ=255)
                                            ggTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCt                                                                                                    >  4:340260/1‑148 (MQ=255)
                                              ttATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTAt                                                                                                 <  7:324014/149‑1 (MQ=255)
                                                aTCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCa                                                                                                                                                   <  1:137514/98‑1 (MQ=255)
                                                aTCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCa                                                                                                                                                   >  2:137514/1‑98 (MQ=255)
                                                    aCAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCa                                                                                           <  5:148878/149‑1 (MQ=255)
                                                         cTGTGGAAGCAGGAGCGTCAAGGGCATAACGCCCGCGAAATTGCCGCTCGGGCCGGCTGGCTCTACGACAATTTCGTGCTGTTCATCCAGGACCTGGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAgcgc                                                                                       <  1:278960/149‑1 (MQ=255)
                                                                          tCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAg                                                                     <  6:209813/149‑1 (MQ=255)
                                                                            aaGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAgg                                                                    <  1:72975/148‑1 (MQ=255)
                                                                                gCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTg                                                               <  7:220705/149‑1 (MQ=255)
                                                                                   gAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGc                                                             >  5:330912/1‑148 (MQ=255)
                                                                                          cgcgAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCTGCa                                                                                           <  4:65448/111‑1 (MQ=255)
                                                                                          cgcgAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCa                                                                                           >  3:65448/1‑111 (MQ=255)
                                                                                           gcgAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGt                                                     >  5:330295/1‑148 (MQ=255)
                                                                                            cgAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTa                                                   <  6:348016/149‑1 (MQ=255)
                                                                                            cgAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTa                                                   <  4:170946/149‑1 (MQ=255)
                                                                                                aTTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTAGTGCTGTTCATCCAGGACCTGGACGGACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCg                                               >  7:311930/1‑149 (MQ=255)
                                                                                                    ccGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGc                                           >  7:148864/1‑149 (MQ=255)
                                                                                                    ccGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGc                                           <  5:65226/149‑1 (MQ=255)
                                                                                                      gAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGa                                         <  3:760/149‑1 (MQ=255)
                                                                                                             ccGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCt                                   >  6:193132/1‑148 (MQ=255)
                                                                                                                 ctggctCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGagcag                                     <  5:299989/142‑1 (MQ=255)
                                                                                                                 ctggctCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGagcag                                     >  6:299989/1‑142 (MQ=255)
                                                                                                                      tctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGctgctg                         >  3:133753/1‑149 (MQ=255)
                                                                                                                           gACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTgg                        <  4:63583/145‑1 (MQ=255)
                                                                                                                           gACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgcgc                    >  3:99292/1‑149 (MQ=255)
                                                                                                                           gACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgcgc                    <  5:193022/149‑1 (MQ=255)
                                                                                                                           gACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTCGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTgg                        >  3:63583/1‑145 (MQ=255)
                                                                                                                                 ttCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGc              <  2:151198/149‑1 (MQ=255)
                                                                                                                                 ttCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGc              <  5:193132/149‑1 (MQ=255)
                                                                                                                                 ttCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGAAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGc              >  3:282610/1‑149 (MQ=255)
                                                                                                                                 ttCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGCCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGc              >  1:151198/1‑149 (MQ=255)
                                                                                                                                      gCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGc          >  1:139831/1‑148 (MQ=255)
                                                                                                                                      ccagTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTg                                                       >  2:59330/4‑103 (MQ=39)
                                                                                                                                      ccagTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTg                                                       <  1:59330/100‑1 (MQ=39)
                                                                                                                                      ccagTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTg                                                       <  1:53867/100‑1 (MQ=39)
                                                                                                                                      ccagTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTg                                                       >  2:53867/4‑103 (MQ=39)
                                                                                                                                        tGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGc          >  1:8378/1‑146 (MQ=255)
                                                                                                                                         gTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCaa        >  5:122059/1‑147 (MQ=255)
                                                                                                                                         gTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGGAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCaa        >  7:44228/1‑147 (MQ=255)
                                                                                                                                            cATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGc    <  6:330295/148‑1 (MQ=255)
                                                                                                                                             aTCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGc          <  3:225613/141‑1 (MQ=255)
                                                                                                                                             aTCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGc          >  4:225613/1‑141 (MQ=255)
                                                                                                                                             aTCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCtt  <  4:99292/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
CAGATCGTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTT  >  NC_002947/1886447‑1886736

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