Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,705,741 C→G 7.1% G21A (GGC→GCC) ‡ PP_5008 ← poly granule‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,705,7410CG 7.1% 195.4 / 4.3 63G21A (GGC→GCC) ‡PP_5008poly granule‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (26/32);  new base G (2/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

GCCGCATCGAGCTCTTTGTCGATGCGCTTCTTGCCGCGCTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTGAGCATTAAGGTAGGGAAGAAAATTAGAA  >  NC_002947/5705594‑5705885
                                                                                                                                                   |                                                                                                                                                
gCCGCATCGAGCTCTTTGTCGATGCGCTTCTTGCCGCGCTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGcc                                                                                                                                                 >  6:146155/1‑149 (MQ=255)
   gCATCGAGCTCTTTGTCGATGCGCTTCTTGCCGCGCTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCcgcg                                                                                                                                              <  6:201870/149‑1 (MQ=255)
   gCATCGAGCTCTTTGTCGATGCGCTTCTTGCCGCGCTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCcgcg                                                                                                                                              <  3:257737/149‑1 (MQ=255)
             ctTTGTCGATGCGCTTCTTGCCGCGCTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGcc                                                                                                                                     >  3:176296/1‑148 (MQ=255)
                 gTCGATGCGCTTCTTGCCGCGCTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGc                                                                                                                                 <  5:146155/148‑1 (MQ=255)
                 gTCGATGCGCTTCTTGCCGCGCTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGc                                                                                                                                 <  2:190895/148‑1 (MQ=255)
                     aTGCGCTTCTTGCCGCGCTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTc                                                                                                                            <  7:327753/149‑1 (MQ=255)
                        cgcTTCTTGCCGCGCTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCa                                                                                                                          <  2:71637/148‑1 (MQ=255)
                        cgcTTCTTGCCGCGCTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAg                                                                                                                         >  8:96777/1‑149 (MQ=255)
                         gcTTCTTGCCGCGCTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGGCTATGCGCGCaagatc                                                                                                                                       >  5:50360/1‑129 (MQ=255)
                         gcTTCTTGCCGCGCTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGGCTATGCGCGCaagatc                                                                                                                                       <  6:50360/134‑6 (MQ=255)
                           ttcttGCCGCGCTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGGCTATGCGCGCaagat                                                                                                                                        >  6:223195/1‑127 (MQ=255)
                           ttcttGCCGCGCTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGGCTATGCGCGCaagat                                                                                                                                        <  5:223195/131‑5 (MQ=255)
                             cttGCCGCGCTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCtt                                                                                                                    >  3:165510/1‑149 (MQ=255)
                                   gcgcTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCt                                                                                                               >  7:266332/1‑148 (MQ=255)
                                   gcgcTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCt                                                                                                               <  8:266332/148‑1 (MQ=255)
                                   gcgcTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCCACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCtt                                                                                                              >  4:110857/1‑149 (MQ=255)
                                     gcTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGTCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCtttt                                                                                                            <  7:324036/149‑1 (MQ=255)
                                     gcTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCtttt                                                                                                            <  4:217504/149‑1 (MQ=255)
                                         ctcGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTttct                                                                                                        >  6:289328/1‑149 (MQ=255)
                                            gACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTttctt                                                                                                       <  4:165510/147‑1 (MQ=255)
                                            gACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTttctt                                                                                                       <  4:176296/147‑1 (MQ=255)
                                                  ttCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCa                                                                                               >  3:279904/1‑149 (MQ=255)
                                                       ccGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTg                                                                                          <  5:186878/149‑1 (MQ=255)
                                                          cccTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGcc                                                                                       <  5:289328/149‑1 (MQ=255)
                                                          cccTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGcc                                                                                       <  6:42440/149‑1 (MQ=255)
                                                                    ctctTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTa                                                                             <  7:96777/149‑1 (MQ=255)
                                                                     tctTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTAc                                                                            >  7:312985/1‑149 (MQ=255)
                                                                       tttGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACct                                                                          >  2:32057/1‑149 (MQ=255)
                                                                            agtagtCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGGAATCACTTTGGCCATGCTGTGTACctcatg                                                                      <  1:32057/148‑1 (MQ=255)
                                                                                gtCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGc                                                                                        <  5:22761/126‑1 (MQ=255)
                                                                                gtCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTATGGCCATGCTGTCTACctcatgcacat                                                                 >  2:320636/1‑149 (MQ=255)
                                                                                 tCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAgg                                                                                                                        <  3:114633/93‑1 (MQ=255)
                                                                                 tCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAgg                                                                                                                        >  4:114633/1‑93 (MQ=255)
                                                                                    gcgcCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAgt                                                             >  2:114497/1‑149 (MQ=255)
                                                                                       ccTTCCTGACCGACGCGGGCGTAGGGGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAgtggt                                                          <  1:320636/149‑1 (MQ=255)
                                                                                         ttCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAgtggtg                                                         >  1:275295/1‑148 (MQ=255)
                                                                                          tCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGtggtgg                                                        >  2:214457/1‑148 (MQ=255)
                                                                                            cTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGaa                                                      >  1:23029/1‑148 (MQ=255)
                                                                                             tGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATGACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGaaaa                                                    >  3:172917/1‑149 (MQ=255)
                                                                                                       gggCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATg                                           >  2:88426/1‑148 (MQ=255)
                                                                                                            tAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCt                                     <  1:214457/149‑1 (MQ=255)
                                                                                                             aGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTg                                    >  1:166213/1‑149 (MQ=255)
                                                                                                             aGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTg                                    >  7:174654/1‑149 (MQ=255)
                                                                                                                  ccGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCtt                               <  8:312985/149‑1 (MQ=255)
                                                                                                                   cGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTg                              <  4:279904/149‑1 (MQ=255)
                                                                                                                   cGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTg                              <  2:275295/149‑1 (MQ=255)
                                                                                                                   cGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATAACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTg                              <  4:172917/149‑1 (MQ=255)
                                                                                                                    gATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTTGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGATGGGCttct                             >  8:54520/1‑147 (MQ=255)
                                                                                                                         cGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTGAGCAt                         >  8:98878/1‑148 (MQ=255)
                                                                                                                                 cagATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTGAGCATTAAGGTAgg                <  2:166213/149‑1 (MQ=255)
                                                                                                                                 cagATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTGAGCATTAAGGTAgg                <  7:98878/149‑1 (MQ=255)
                                                                                                                                   gATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTGAGc                           <  7:69698/136‑1 (MQ=255)
                                                                                                                                   gATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTGAGc                           >  8:69698/1‑136 (MQ=255)
                                                                                                                                     tCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTGAGCATTAAGGTAggg               >  8:35541/1‑146 (MQ=255)
                                                                                                                                       ttGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGtggtgg                                                        >  8:110609/1‑103 (MQ=255)
                                                                                                                                           gcgcATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTGAGCATTAAGGTAGGGAAGAAAAt       >  5:183097/1‑148 (MQ=255)
                                                                                                                                            cgcATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTGAGCATTAAGGTAGGGAAGAAAATTa     <  7:35541/149‑1 (MQ=255)
                                                                                                                                            cgcATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTGAGCATTAAGGTAGGGAAGAAAATTa     <  5:239716/149‑1 (MQ=255)
                                                                                                                                             gcATAGCCGCGCAGCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTGAGCATTAAGGTAGGGAAGAAAATTAg    <  3:110857/149‑1 (MQ=255)
                                                                                                                                              cATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTGAGCATTAAGGTAGGGAAGAAAATTa     <  5:231945/147‑1 (MQ=255)
                                                                                                                                              cATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAAACACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTGAGCATTAAGGTAGGGAAGAAAATTa     >  6:231945/1‑147 (MQ=255)
                                                                                                                                               aTAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTGAGCATTAAGGTAGGGAAGAAAATTAGaa  <  8:174654/149‑1 (MQ=255)
                                                                                                                                               aTAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTGAGCATTAAGGTAGGGAAGAAAATTAGaa  <  1:310046/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCCGCATCGAGCTCTTTGTCGATGCGCTTCTTGCCGCGCTTCTCGACACCTTCACCCGCCCTGACCAGCTCTTTGAAGTAGTCAGCGCCTTCCTGACCGACGCGGGCGTAGGCGCCGATACCGGCCAGCCAGATCTTGCGCGCATAGCCGCGCACCTCGCCAAGCGTGTCCAGGGCTTCGTCTTTTTTCTTCGCAATCACTTTGGCCATGCTGTGTACCTCATGCTCATGAGTGGTGGAAAACCGCCCATGGAGGCTGGGCTTGAGCATTAAGGTAGGGAAGAAAATTAGAA  >  NC_002947/5705594‑5705885

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