Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,784,033 T→C 5.3% Q15R (CAG→CGG) ‡ glnD ← uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,784,0330TC 5.3% 223.0 / 2.3 76Q15R (CAG→CGG) ‡glnDuridylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base C (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

CAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGTGGC  >  NC_002947/1783885‑1784180
                                                                                                                                                    |                                                                                                                                                   
cAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCt                                                                                                                                                     >  6:115624/1‑149 (MQ=255)
 aGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTg                                                                                                                                                    <  5:115624/149‑1 (MQ=255)
               gaGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGt                                                                                                                                      >  3:215046/1‑149 (MQ=255)
                 gCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACTGTc                                                                                                                                     >  1:113851/1‑148 (MQ=255)
                    aGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTAGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGCGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGaa                                                                                                                                  >  8:319163/1‑148 (MQ=255)
                     ggCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGaa                                                                                                                                  >  6:344950/1‑147 (MQ=255)
                             ccTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCgggg                                                                                                                        <  5:184210/149‑1 (MQ=255)
                              cTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCTATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGt                                                                                                                       >  6:26038/1‑149 (MQ=255)
                              cTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTg                                                                                                                                                    >  6:359758/1‑120 (MQ=255)
                              cTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTg                                                                                                                                                    <  5:359758/120‑1 (MQ=255)
                              cTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGaa                                                                                                                                                 >  7:163178/1‑123 (MQ=255)
                              cTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGaa                                                                                                                                                 <  8:163178/123‑1 (MQ=255)
                              cTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCgggg                                                                                                                        >  5:55464/1‑148 (MQ=255)
                              cTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCgggg                                                                                                                        <  6:55464/148‑1 (MQ=255)
                                  aTCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAc                                                                                                                   <  2:272259/149‑1 (MQ=255)
                                   tCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCgggg                                                                                                                        >  1:66680/1‑143 (MQ=255)
                                   tCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCgggg                                                                                                                        <  2:66680/143‑1 (MQ=255)
                                      gCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGt                                                                                                                       <  2:232678/141‑1 (MQ=255)
                                      gCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGt                                                                                                                       >  1:232678/1‑141 (MQ=255)
                                       cGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTggg                                                                                                              >  8:178945/1‑149 (MQ=255)
                                             gATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATcg                                                                                                        >  8:110577/1‑149 (MQ=255)
                                              aTGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCAc                                                                                                                                         >  4:253884/1‑115 (MQ=255)
                                              aTGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCAc                                                                                                                                         <  3:253884/115‑1 (MQ=255)
                                               tGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCAGGCGGAACTgg                                                                                                                                             >  3:7242/1‑110 (MQ=255)
                                               tGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCAGGCGGAACTgg                                                                                                                                             <  4:7242/110‑1 (MQ=255)
                                                   gCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGaa                                                                                                                                                 <  5:1870/102‑1 (MQ=255)
                                                   gCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGaa                                                                                                                                                 >  6:1870/1‑102 (MQ=255)
                                                       ccGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTc                                                                                              >  2:178604/1‑149 (MQ=255)
                                                         ggCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCt                                                                                            <  7:178945/149‑1 (MQ=255)
                                                                    gTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTg                                                                                 >  7:109675/1‑149 (MQ=255)
                                                                    gTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTg                                                                                 >  6:350223/1‑149 (MQ=255)
                                                                    ctttGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCCCGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTg                                                                                 >  7:227193/2‑149 (MQ=255)
                                                                       tatCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGc                                                                               <  4:116118/146‑1 (MQ=255)
                                                                           gtgCACCTGGCGGGCCAGGCGGATGGCTTTGTTGAAGGCGGGGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGaa                                                                                                                                  <  6:198059/91‑1 (MQ=255)
                                                                           gAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCTGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGaa                                                                                                                                  >  5:198059/1‑93 (MQ=255)
                                                                            aGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCg                                                                         <  4:215046/149‑1 (MQ=255)
                                                                                 cTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTg                                                                                           <  4:32731/126‑1 (MQ=255)
                                                                                 cTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTg                                                                                           >  3:32731/1‑126 (MQ=255)
                                                                                 cTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTg                                                                    >  6:349228/1‑149 (MQ=255)
                                                                                 cTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTg                                                                    >  1:305268/1‑149 (MQ=255)
                                                                                          cAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCggg                                                                       >  1:48338/1‑137 (MQ=255)
                                                                                          cAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCggg                                                                       <  2:48338/137‑1 (MQ=255)
                                                                                           aGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCg                                                                                                                           <  1:243087/84‑1 (MQ=255)
                                                                                           aGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCg                                                                                                                           >  2:243087/1‑84 (MQ=255)
                                                                                             gCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGa                                                        <  8:227193/149‑1 (MQ=255)
                                                                                                    gCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTgcgc                                                 <  2:305268/149‑1 (MQ=255)
                                                                                                    gCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGGTCCGCCGGGAACTGGCCACGGTCGAACAGCTCGGGGTCCCCCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGGCGGGGTGCGGTCAGGCCGAAGTgcg                                                  <  5:358752/148‑1 (MQ=255)
                                                                                                      tttCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCgg                                               >  6:96947/1‑149 (MQ=255)
                                                                                                         cTTGAAGGCGGCGATGGTGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTg                                                                                 <  2:220566/112‑1 (MQ=255)
                                                                                                          ttGAAGGGGTCGATGGGTCTCGCCTTGAGGGCCAGTTCCGCATGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTTGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGca                                                                              <  6:136209/114‑2 (MQ=255)
                                                                                                          ttGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGcc                                                                              >  5:136209/1‑114 (MQ=255)
                                                                                                            gacgtcggcgATGGGGCTCGCCTTGAGGGCCAGGTCCGCCTGGAACGGGCCACGGTCGTACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGGGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATg                                           <  5:344950/142‑1 (MQ=255)
                                                                                                            gAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGAt                                            <  1:188488/146‑1 (MQ=255)
                                                                                                            gAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGgt                                         >  7:5932/1‑149 (MQ=255)
                                                                                                             aaGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGgtg                                        >  8:31522/1‑149 (MQ=255)
                                                                                                             aaGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGgt                                         <  5:96947/148‑1 (MQ=255)
                                                                                                             aaGGCGGCGATGGGGCGCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCGCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGgt                                         <  7:319163/148‑1 (MQ=255)
                                                                                                              aGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGgtgt                                       >  4:240665/1‑149 (MQ=255)
                                                                                                               ggcggcGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTgtc                                      <  5:350223/149‑1 (MQ=255)
                                                                                                                          ggCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTgcg                             >  3:6127/1‑147 (MQ=255)
                                                                                                                           gCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGgt                                         >  4:121319/1‑134 (MQ=255)
                                                                                                                            cTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTgtcgt                                    >  6:110822/1‑138 (MQ=255)
                                                                                                                            cTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTgtcgt                                    <  5:110822/138‑1 (MQ=255)
                                                                                                                            cTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTgtcgt                                    >  2:169289/1‑138 (MQ=255)
                                                                                                                            cTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTgtcgt                                    <  1:169289/138‑1 (MQ=255)
                                                                                                                               gCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAggg                        <  7:31522/147‑1 (MQ=255)
                                                                                                                                 gaccgtggCCAGTTCCAGGCGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTg                      >  1:144938/7‑147 (MQ=255)
                                                                                                                                 cTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGaa                    <  3:240665/149‑1 (MQ=255)
                                                                                                                                        gCCAGTTCCAGGCGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGAtctc              <  2:144938/148‑1 (MQ=255)
                                                                                                                                         ccAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCg             >  4:25536/1‑148 (MQ=255)
                                                                                                                                         ccAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCg             >  1:21286/1‑148 (MQ=255)
                                                                                                                                                cGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGt     >  4:126281/1‑149 (MQ=255)
                                                                                                                                                cGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGt     <  3:25536/149‑1 (MQ=255)
                                                                                                                                                 gCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGTg    <  5:349228/149‑1 (MQ=255)
                                                                                                                                                  ccTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGTgg   <  3:126281/149‑1 (MQ=255)
                                                                                                                                                  ccTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGTgg   <  8:5932/149‑1 (MQ=255)
                                                                                                                                                   cTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGTGGc  >  6:20139/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGTGGC  >  NC_002947/1783885‑1784180

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