Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 159,001 C→G 6.0% G463A (GGC→GCC) ‡ PP_0149 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947159,0010CG 6.0% 210.5 / 4.3 67G463A (GGC→GCC) ‡PP_0149hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/29);  new base G (2/2);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

GCGCAGGTCAGCCTTGGCGAACCAGCTACGGCGCTGCTCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTGCGCAGCAGCTTCAAC  >  NC_002947/158853‑159147
                                                                                                                                                    |                                                                                                                                                  
gcgcAGGTCAGCCTTGGCGAACCAGCTACGGCGCTGCTCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGc                                                                                                                                                    >  6:310451/1‑149 (MQ=255)
   cAGGGCAGCCTTGGCGAACCAGCTACGGCGCTGCTCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTTGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCa                                                                                                                                                 <  1:294523/149‑1 (MQ=255)
     ggTCAGCCTTGGCGAACCAGCTACGGCGCTGCTCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAgg                                                                                                                                               <  2:134743/149‑1 (MQ=255)
            cTTGGCGAACCAGCTACGGCGCTGCTCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACg                                                                                                                                        >  6:186310/1‑149 (MQ=255)
                gCGAACCAGCTACGGCGCTGCTCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTc                                                                                                                                    <  5:144502/149‑1 (MQ=255)
                 cGAACCAGCTACGGCGCTGCTCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGTTACGCCTGGGCTACCGCTTcg                                                                                                                                   <  4:399324/149‑1 (MQ=255)
                             ggCGCTGCTCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATc                                                                                                                        >  3:113531/1‑148 (MQ=255)
                                gctgctCGGCATTGAGGTCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTGAGCACCTCGAACATGGCCTGGTACCGGTCTCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGc                                                                                                                    <  5:285208/149‑1 (MQ=255)
                                  tgctCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGTTACGCCTGGGCTACCGCTTcga                                                                                                                                  >  6:69347/1‑132 (MQ=255)
                                  tgctCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGTTACGCCTGGGCTACCGCTTcga                                                                                                                                  <  5:69347/133‑2 (MQ=255)
                                   gctCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTAACGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCAc                                                                                                                  >  3:240617/1‑148 (MQ=255)
                                       ggCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTg                                                                                                              <  3:268835/148‑1 (MQ=255)
                                            tGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAg                                                                                                                                                <  6:110165/109‑1 (MQ=255)
                                            tGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAg                                                                                                                                                >  5:110165/1‑109 (MQ=255)
                                                ggCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGc                                                                                                      >  1:118532/1‑147 (MQ=255)
                                                ggCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCa                                                                                                    <  8:8304/149‑1 (MQ=255)
                                                 gccagccaGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCt                                                                                                               >  8:369444/1‑137 (MQ=255)
                                                 gccagccaGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCt                                                                                                               <  7:369444/137‑1 (MQ=255)
                                                      ccagTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTcccc                                                                                              >  5:63925/1‑149 (MQ=255)
                                                        agTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTg                                                                                            <  6:63925/149‑1 (MQ=255)
                                                          tCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGcc                                                                                          <  5:310451/149‑1 (MQ=255)
                                                            ttCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCt                                                                                        >  4:109439/1‑149 (MQ=255)
                                                               aaCCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCAGCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGa                                                                                     <  4:369247/149‑1 (MQ=255)
                                                               aaCCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGa                                                                                     <  3:109439/149‑1 (MQ=255)
                                                                  cAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCa                                                                                  >  3:24187/1‑149 (MQ=255)
                                                                   agagCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCAc                                                                                 <  4:113531/149‑1 (MQ=255)
                                                                     agCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACgg                                                                               >  5:402625/1‑149 (MQ=255)
                                                                       cttgcttgATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGc                                                                              >  6:313701/1‑148 (MQ=255)
                                                                            ttgATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTg                                                                                                       <  1:55562/118‑1 (MQ=255)
                                                                            ttgATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTg                                                                                                       >  2:55562/1‑118 (MQ=255)
                                                                            ttgATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAAACGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCa                                                                         >  1:116675/1‑148 (MQ=255)
                                                                                ttGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGTAAACCg                                                                    >  6:146712/1‑149 (MQ=255)
                                                                                   ggaaCTGCCTTCTCCCCCCCTCGCACATGGCCGGGTACCGGTCCCCGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTg                                                                 <  3:26047/147‑1 (MQ=255)
                                                                                   gTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTg                                                                 >  5:154628/1‑149 (MQ=255)
                                                                                        tGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCat                                                            >  7:100006/1‑149 (MQ=255)
                                                                                          cTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCatga                                                          <  8:57030/149‑1 (MQ=255)
                                                                                          cTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCatga                                                          >  6:390155/1‑149 (MQ=255)
                                                                                          cTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCatg                                                           <  6:402625/148‑1 (MQ=255)
                                                                                              ctcAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGcc                                                                             <  1:21323/126‑1 (MQ=255)
                                                                                              ctcAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGcc                                                                             >  2:21323/1‑126 (MQ=255)
                                                                                              ctcAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCt                                                      >  5:411446/1‑149 (MQ=255)
                                                                                                  gCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTc                                                   >  6:92395/1‑148 (MQ=255)
                                                                                                      cTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGt                                              >  1:156265/1‑149 (MQ=255)
                                                                                                          aaCATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGa                                                                                     >  8:161209/1‑106 (MQ=255)
                                                                                                          aaCATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGa                                                                                     <  7:161209/106‑1 (MQ=255)
                                                                                                          aaCATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGtt                                          <  2:306710/149‑1 (MQ=255)
                                                                                                             aTGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAg                                       >  2:289802/1‑149 (MQ=255)
                                                                                                             aTGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAg                                       <  8:100006/149‑1 (MQ=255)
                                                                                                                  cTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTg                                  >  7:165348/1‑149 (MQ=255)
                                                                                                                     gTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTgg                                 <  5:186310/147‑1 (MQ=255)
                                                                                                                       aCCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAg                                                >  8:328177/1‑130 (MQ=255)
                                                                                                                       aCCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAg                                                <  7:328177/130‑1 (MQ=255)
                                                                                                                        ccGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTg                             <  2:156265/148‑1 (MQ=255)
                                                                                                                           gTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGt                                           >  6:2553/1‑131 (MQ=255)
                                                                                                                              gcgGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGc                      <  6:154628/149‑1 (MQ=255)
                                                                                                                                 gTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGc                                                                                                                    >  3:90244/1‑52 (MQ=255)
                                                                                                                                 gTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGc                                                                                                                    <  4:90244/52‑1 (MQ=255)
                                                                                                                                  tAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCt                  >  5:57731/1‑149 (MQ=255)
                                                                                                                                   aaGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTg                 >  1:373355/1‑149 (MQ=255)
                                                                                                                                     gAATCGAAGCGGTAGACCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATTATGATGTCCAGGTCGTTGAGGAGTTGATTGGCGGTTCCGCt                  >  6:286808/1‑146 (MQ=255)
                                                                                                                                      aaTCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTgcg               >  3:312269/1‑148 (MQ=255)
                                                                                                                                      aaTCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTgcg               >  2:66090/1‑148 (MQ=255)
                                                                                                                                       cgcGAAGCGTTACGCCTGGGCTACCGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTgcg               >  4:37869/3‑147 (MQ=255)
                                                                                                                                        tCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTGCgcag            >  4:258914/1‑149 (MQ=255)
                                                                                                                                         cGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTGCgcag            >  8:47666/1‑148 (MQ=255)
                                                                                                                                              cGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTGCGCAGCAGCtt      <  5:313701/149‑1 (MQ=255)
                                                                                                                                                 tAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGTTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGAAACCCGTTTTGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTGCGCAGCAGCTTaaa   <  5:286808/149‑1 (MQ=255)
                                                                                                                                                   gCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTGCGCAGCAGCTTCAAc  <  7:47666/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCGCAGGTCAGCCTTGGCGAACCAGCTACGGCGCTGCTCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTGCGCAGCAGCTTCAAC  >  NC_002947/158853‑159147

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