Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,978,708 C→G 7.6% F45L (TTC→TTG)  cmk → cytidylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,978,7080CG 7.6% 201.5 / 7.4 66F45L (TTC→TTG) cmkcytidylate kinase
Reads supporting (aligned to +/- strand):  ref base C (33/28);  new base G (3/2);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

AAGAGGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATCCG  >  NC_002947/1978561‑1978854
                                                                                                                                                    |                                                                                                                                                  
aagaGGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTCCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTc                                                                                                                                                    >  2:135148/1‑148 (MQ=255)
  gagGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAAt                                                                                                                                                 >  3:179418/1‑149 (MQ=255)
     ggCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTTCAGGACTCCGGCGCGCAGTATCGCTTGCTGGCATTCAATga                                                                                                                                               >  8:314624/1‑147 (MQ=255)
      gCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCa                                                                                                                                             <  6:149667/149‑1 (MQ=255)
          gTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCaa                                                                                                                                         <  4:179418/149‑1 (MQ=255)
             gTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCa                                                                                                                                      >  3:39538/1‑149 (MQ=255)
              tGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCAc                                                                                                                                     >  7:167158/1‑149 (MQ=255)
               gAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCTGCCGGGGTGATTTCCCCCGAGCGGGGGTGGAACATGGTGGCATCCGGCGCGCGGTAACGCTTGATGGCATTCAATGCCAGCAACCAc                                                                                                                                     >  1:106639/1‑149 (MQ=255)
                aaTTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACg                                                                                                                                    <  5:47230/148‑1 (MQ=255)
                 atacAATGGCACCGGATATCACCATCGACGGACCAAGCGGCTCCCGCAAGGCCACGGTC‑GCCGGGCTTATTGCCCGCGACCTGGGCTGGAAGCTGCTTGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGc                                                                                                                                  <  6:117167/146‑1 (MQ=255)
                     aaTGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTc                                                                                                                               <  2:7770/148‑1 (MQ=255)
                         gCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCt                                                                                                                          <  4:39538/149‑1 (MQ=255)
                          cACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCAGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCt                                                                                                                          >  1:42324/1‑148 (MQ=255)
                           aCCGGTTATCACCATCGACGGACCTAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGa                                                                                                                        <  7:314624/149‑1 (MQ=255)
                           aCCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCACTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGa                                                                                                                        <  1:224737/149‑1 (MQ=255)
                               gTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCaa                                                                                                                    >  3:112102/1‑149 (MQ=255)
                                   tCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACg                                                                                                                  >  5:352781/1‑147 (MQ=255)
                                      ccATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGCCCTGACCAACGAagag                                                                                                             >  6:53201/1‑149 (MQ=255)
                                      ccATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAagag                                                                                                             >  3:261718/1‑149 (MQ=255)
                                        aTCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgct                                                                                                           <  3:183459/149‑1 (MQ=255)
                                           gACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAAc                                                                                                                                        >  8:141641/1‑117 (MQ=255)
                                           gACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAAc                                                                                                                                        <  7:141641/117‑1 (MQ=255)
                                            aCGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgctgct                                                                                                        >  7:180551/1‑148 (MQ=255)
                                             ctgACCTAGCGGGTCGGGCAAGGGCAGGGTC‑GCCTGGGTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTTCTTGCATTCAATGCCAGCAACCACGGCGTAGACCTTACCAACGAAGAGctgctg                                                                                                       <  2:42324/146‑1 (MQ=255)
                                              ggACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCTACCACGGCGTCGGCCTGACCAACGAAGAGCTGCTGAc                                                                                                     >  6:31578/1‑149 (MQ=255)
                                                  cAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGACCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGa                                                                                                                        >  1:73376/1‑126 (MQ=255)
                                                  cAAGCGGAGCGGGCAAGGGCACGGTC‑GCCTGGCTGATTGCCCGCGAGCTGGTCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGa                                                                                                                        <  2:73376/126‑1 (MQ=255)
                                                    aGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAAc                                                                                                                   <  6:95733/129‑1 (MQ=255)
                                                    aGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAAc                                                                                                                   >  5:95733/1‑129 (MQ=255)
                                                          tCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGTACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTgccgc                                                                                         <  8:167158/149‑1 (MQ=255)
                                                                  gggCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTg                                                                                  >  3:357103/1‑148 (MQ=255)
                                                                   ggCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGa                                                                                >  7:60684/1‑149 (MQ=255)
                                                                      aCGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCGTACCAAGCTTGCCGCCCAGCTGGAtgt                                                                             >  3:196830/1‑149 (MQ=255)
                                                                        ggTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGAc                                                                                                                            >  8:72711/1‑100 (MQ=255)
                                                                        ggTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGAc                                                                                                                            <  7:72711/100‑1 (MQ=255)
                                                                          tC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAg                                                                         >  7:271775/1‑149 (MQ=255)
                                                                          tC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAg                                                                         <  4:101914/149‑1 (MQ=255)
                                                                           c‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGc                                                                                                                                           >  7:223502/1‑82 (MQ=255)
                                                                           c‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGc                                                                                                                                           <  8:223502/82‑1 (MQ=255)
                                                                             gCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGt                                                                        >  3:14031/1‑148 (MQ=255)
                                                                                      aTTGGCCGCGAGCTGGGCTGGATGCTGCTGGACTCCGGCGCTCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAATAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCt                                                                <  5:335063/147‑1 (MQ=255)
                                                                                       ttGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGcc                                                             <  8:60684/149‑1 (MQ=255)
                                                                                             gcgAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGccc                                                       <  4:357103/149‑1 (MQ=255)
                                                                                              cgAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCg                                                      <  6:352781/149‑1 (MQ=255)
                                                                                                 gctgggctggAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTa                                                   <  2:106639/149‑1 (MQ=255)
                                                                                                   tgggctggAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAg                                                 >  3:93714/1‑149 (MQ=255)
                                                                                                            aGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTacaa                                           >  1:158873/1‑146 (MQ=255)
                                                                                                             gctgctGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGa                                       <  1:135148/149‑1 (MQ=255)
                                                                                                             gctgctGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCAAATGGATGTGCAGTTCATCGCTGCCGAGCCCg                                                      >  6:10662/1‑134 (MQ=255)
                                                                                                                       tCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGCCCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCtt                               >  7:33310/1‑147 (MQ=255)
                                                                                                                         cGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAgg                           >  4:112384/1‑149 (MQ=255)
                                                                                                                         cGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAgg                           <  1:10565/149‑1 (MQ=255)
                                                                                                                          ggCGCGCTGTATCACTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAtgtg                                                                            <  5:222063/99‑1 (MQ=255)
                                                                                                                          ggCGCGCTGTATCACTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAtgtg                                                                            >  6:222063/1‑99 (MQ=255)
                                                                                                                             gcgcTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGa                       <  4:93714/149‑1 (MQ=255)
                                                                                                                               gcTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAgg                     <  4:14031/149‑1 (MQ=255)
                                                                                                                               gcTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAg                      >  6:123027/1‑148 (MQ=255)
                                                                                                                                   tATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGt                 <  4:112102/149‑1 (MQ=255)
                                                                                                                                           gCTGGCATTGAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGCGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACg         >  8:584/1‑149 (MQ=255)
                                                                                                                                            cTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACg         <  2:158873/148‑1 (MQ=255)
                                                                                                                                            cTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGt        >  7:221333/1‑149 (MQ=255)
                                                                                                                                             tGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGt        >  7:293175/1‑148 (MQ=255)
                                                                                                                                                 aTTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATcc   >  7:245968/1‑149 (MQ=255)
                                                                                                                                                 aTTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGATGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATcc   >  7:194748/1‑149 (MQ=255)
                                                                                                                                                  ttCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGa                    <  2:330537/131‑1 (MQ=255)
                                                                                                                                                  ttCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGa                    >  1:330537/1‑131 (MQ=255)
                                                                                                                                                  ttCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATCCg  >  6:264378/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAGAGGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTC‑GCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATCCG  >  NC_002947/1978561‑1978854

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