Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,353,367 C→T 5.9% E250E (GAG→GAA)  pnp ← polyribonucleotide nucleotidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,353,3670CT 5.9% 258.8 / 3.7 84E250E (GAG→GAA) pnppolyribonucleotide nucleotidyltransferase
Reads supporting (aligned to +/- strand):  ref base C (44/35);  new base T (2/3);  total (46/38)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.90e-01

CCGAAGATCTCTTTTACTTCGGAAGCCGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAGTACGGCACCCAGCATCT  >  NC_002947/5353221‑5353515
                                                                                                                                                  |                                                                                                                                                    
ccGAAGATCTCTTTTACTTCGGAAGCCGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTc                                                                                                                                                    >  5:311385/1‑149 (MQ=255)
   aaGATCTCTTTTACTTCGGAAGCCGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAGTTCGgc                                                                                                                                                 <  1:53893/149‑1 (MQ=255)
   aaGATCTCTTTTACTTCGGAAGCCGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGgc                                                                                                                                                 >  5:141406/1‑149 (MQ=255)
           ttttACTTCGGAAGCCGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAGTTCGgcgaag                                                                                                                                             <  6:65173/145‑4 (MQ=255)
           ttttACTTCGGAAGCCGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGa                                                                                                                                            <  5:72628/146‑1 (MQ=255)
           ttttACTTCGGAAGCCGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAg                                                                                                                                         >  1:262366/1‑149 (MQ=255)
           ttttACTTCGGAAGCCGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAg                                                                                                                                         >  5:37218/1‑149 (MQ=255)
             ttagttagGATGCCGAAGTCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATa                                                                                                                                          <  6:224743/146‑1 (MQ=255)
              tACTTCGGAAGCCGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCa                                                                                                                                       <  6:141406/148‑1 (MQ=255)
               aCTTCGGAAGCCGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCAtt                                                                                                                                     >  8:84140/1‑149 (MQ=255)
                cTTCGGAAGCCGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCAtt                                                                                                                                     <  1:263370/148‑1 (MQ=255)
                   cGGAAGCCGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGaa                                                                                                                                  <  2:4741/148‑1 (MQ=255)
                     gAAGCCGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAGTTCGgcgaagc                                                                                                                                            <  5:235874/136‑5 (MQ=255)
                     gAAGCCGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAGTTCGgcgaagc                                                                                                                                            >  6:235874/1‑132 (MQ=255)
                         ccGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCg                                                                                                                           >  6:300086/1‑149 (MQ=255)
                         ccGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCg                                                                                                                           <  2:30450/149‑1 (MQ=255)
                           gAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCTTTGAACAGCTCGgt                                                                                                                         <  2:193945/149‑1 (MQ=255)
                               gCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTg                                                                                                                     <  5:300086/149‑1 (MQ=255)
                               gCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTg                                                                                                                     >  3:210409/1‑149 (MQ=255)
                                 tGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGATCGGTGTTgg                                                                                                                    >  1:221168/1‑148 (MQ=255)
                                     cttcttCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACg                                                                                                               >  5:181964/1‑149 (MQ=255)
                                         ttcACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAg                                                                                                           <  3:221581/149‑1 (MQ=255)
                                           cACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCa                                                                                                            <  2:262366/146‑1 (MQ=255)
                                            aCCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGtt                                                                                                        >  5:175892/1‑149 (MQ=255)
                                            aCCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGGGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTTGCAACGGCAGGt                                                                                                         <  6:181964/148‑1 (MQ=255)
                                             ccGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTg                                                                                                                     <  6:216203/135‑1 (MQ=255)
                                             ccGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTg                                                                                                                     >  5:216203/1‑135 (MQ=255)
                                              cGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTc                                                                                                      >  2:265732/1‑149 (MQ=255)
                                              cGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTc                                                                                                      <  2:10879/149‑1 (MQ=255)
                                               ggAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTc                                                                                                      >  3:301085/1‑148 (MQ=255)
                                                gagaACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCa                                                                                                    <  4:301085/149‑1 (MQ=255)
                                                     cTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTccc                                                                                               <  4:210409/149‑1 (MQ=255)
                                                        ggCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTcc                                                                                                >  7:115303/1‑145 (MQ=255)
                                                        ggCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTcc                                                                                                <  8:115303/145‑1 (MQ=255)
                                                         gCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGt                                                                                           >  8:209424/1‑149 (MQ=255)
                                                           aaCCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCg                                                                                         >  2:2474/1‑149 (MQ=255)
                                                                  tGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGGTTCCAGTCCCAAGTCGGCTTGcc                                                                                  >  3:255982/1‑149 (MQ=255)
                                                                    gTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAg                                                                                >  6:1218/1‑149 (MQ=255)
                                                                       aCGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAGCCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCtt                                                                             >  4:228379/1‑149 (MQ=255)
                                                                              tCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGcagc                                                                      >  4:23348/1‑149 (MQ=255)
                                                                                aCCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGcagc                                                                      >  2:280156/1‑147 (MQ=255)
                                                                                 cccAGACGCGCGTAGCGGTCGGCTTTGACGGTGTTGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGc                                                                               >  5:53133/1‑137 (MQ=255)
                                                                                 cccAGACGCGCGTAGCGGTCGGCTTTGACGGTGTTGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGc                                                                               <  6:53133/137‑1 (MQ=255)
                                                                                 cccAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGcagccag                                                                   >  2:20095/1‑149 (MQ=255)
                                                                                  ccAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGcagccagc                                                                  >  1:57743/1‑149 (MQ=255)
                                                                                       cgcgcgTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCttt                                                             <  6:175892/149‑1 (MQ=255)
                                                                                          gcgTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCt                                                                                      >  6:210948/1‑121 (MQ=255)
                                                                                              aGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGcc                                                      >  5:35709/1‑149 (MQ=255)
                                                                                               gCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGcc                                                      <  1:20095/148‑1 (MQ=255)
                                                                                                  gTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTATCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGa                                                  >  5:189077/1‑149 (MQ=255)
                                                                                                   tCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGAt                                                 <  7:209424/149‑1 (MQ=255)
                                                                                                      gCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAAc                                              <  5:1218/149‑1 (MQ=255)
                                                                                                       cTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGTTAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACa                                             <  2:19703/149‑1 (MQ=255)
                                                                                                              ggTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTgg                                       >  5:262955/1‑148 (MQ=255)
                                                                                                               gTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGaa                                     <  1:280156/149‑1 (MQ=255)
                                                                                                               gTGACGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCCGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGaa                                     <  7:283410/149‑1 (MQ=255)
                                                                                                                tGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAAt                                    >  3:238119/1‑149 (MQ=255)
                                                                                                                 gATGGTGTAACCCTGCGAAACCGCTTCGCCTAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGTCAACGGCATGTTTACAGTCACAAATAGGCTTGCAAGCTTAGGCATCAAGCTCTTTTACAGCCTGGATAACAGACTGGAAt                                    >  4:194610/1‑148 (MQ=255)
                                                                                                                  aTGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATtc                                  >  7:240492/1‑149 (MQ=255)
                                                                                                                   tGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTACCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGAt                                                 <  1:44443/133‑1 (MQ=255)
                                                                                                                   tGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTACCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGAt                                                 >  2:44443/1‑133 (MQ=255)
                                                                                                                      tgtAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGc                                           >  4:72715/1‑136 (MQ=255)
                                                                                                                      tgtAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGc                                           <  3:72715/136‑1 (MQ=255)
                                                                                                                      tgtAACACTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCATTGACAGCGTGGATAACGGCCTGGAATtcgtca                              >  4:277864/1‑148 (MQ=255)
                                                                                                                          aCCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTggg                          <  3:251886/149‑1 (MQ=255)
                                                                                                                           cccTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGTATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGc                         <  6:37218/149‑1 (MQ=255)
                                                                                                                           cccTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGc                         >  3:113039/1‑149 (MQ=255)
                                                                                                                              tGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCttt                                                             >  3:22100/1‑110 (MQ=255)
                                                                                                                              tGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCttt                                                             <  4:22100/110‑1 (MQ=255)
                                                                                                                              tGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGaa                      >  6:180527/1‑149 (MQ=255)
                                                                                                                                  aaaCCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAg                   >  4:223019/1‑148 (MQ=255)
                                                                                                                                     ccGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGaa                                     <  1:10279/127‑1 (MQ=255)
                                                                                                                                     ccGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGaa                                     >  2:10279/1‑127 (MQ=255)
                                                                                                                                     ccGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATtcgtc                               <  4:87720/133‑1 (MQ=255)
                                                                                                                                     ccGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATtcgtc                               >  3:87720/1‑133 (MQ=255)
                                                                                                                                     ccGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAGTACg               >  3:51814/1‑149 (MQ=255)
                                                                                                                                       tccgcGCCGAGTTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAGTACGGc             >  3:90021/5‑149 (MQ=255)
                                                                                                                                       gCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGc                                                       >  3:41288/1‑107 (MQ=255)
                                                                                                                                       gCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGc                                                       <  4:41288/107‑1 (MQ=255)
                                                                                                                                       gCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAGTACGGc             <  4:113039/149‑1 (MQ=255)
                                                                                                                                       gCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAGTACGGc             <  1:2474/149‑1 (MQ=255)
                                                                                                                                             ccGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAGTACGGCACCCAg       >  2:7453/1‑149 (MQ=255)
                                                                                                                                              cGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAGTACGGCACCCAg       >  6:307288/1‑148 (MQ=255)
                                                                                                                                                 aCTCGGCGCGGATAGCATTGAACAGCTCGGGGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAGTACGGCACCCAGCATc   <  6:311385/149‑1 (MQ=255)
                                                                                                                                                  cTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAGTACGGCACCCGGCATCt  >  5:64376/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCGAAGATCTCTTTTACTTCGGAAGCCGAAGGCTGGCCTTCTTCACCGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAGTACGGCACCCAGCATCT  >  NC_002947/5353221‑5353515

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