Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,502,358 G→T 6.3% intergenic (+150/‑14) yccA → / → xanP FtsH‑mediated proteolysis substrate/modulator/xanthine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,502,3580GT 6.3% 199.6 / 3.8 64intergenic (+150/‑14)yccA/xanPFtsH‑mediated proteolysis substrate/modulator/xanthine permease
Reads supporting (aligned to +/- strand):  ref base G (34/26);  new base T (2/2);  total (36/28)
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.72e-01

GATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTGGCAGCCCTGCAGCACCTGCTGGCGATCATCGTGC  >  NC_002947/4502211‑4502501
                                                                                                                                                   |                                                                                                                                                
gattgCTGTACCGCAAGAAGCGCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCggg                                                                                                                                                  <  6:176295/148‑1 (MQ=255)
  ttgCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGggcg                                                                                                                                               >  2:110711/1‑149 (MQ=255)
       gTACCGCAAGAAGCCCGCTTCGGCCGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGATCGTGCCGTCATCCCGTTTTCACTACAGCTATCAGCCATAACCCCGAAGAATGCGCTCCTTTTTTCTTCCTGGGCGGCttt                                                                                                                                          >  8:219959/1‑149 (MQ=255)
          ccGCTAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTGATCCAGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCa                                                                                                                                        <  8:107994/148‑1 (MQ=255)
          ccGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAg                                                                                                                                       >  4:19767/1‑149 (MQ=255)
              aagaagCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAAt                                                                                                                                   <  6:114260/149‑1 (MQ=255)
                       gCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTcaca                                                                                                                          >  2:68971/1‑149 (MQ=255)
                       gCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGACGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTcaca                                                                                                                          <  5:171343/149‑1 (MQ=255)
                        cTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCacac                                                                                                                         <  1:68971/149‑1 (MQ=255)
                            ggcggGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGaa                                                                                                                      >  7:287910/1‑148 (MQ=255)
                               gggCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAAcaca                                                                                                                  <  8:214992/149‑1 (MQ=255)
                                  cTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGcc                                                                                                               <  8:166327/149‑1 (MQ=255)
                                      tttGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAgg                                                                                                           >  6:126364/1‑149 (MQ=255)
                                               gcTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGgc                                                                                                  >  3:79909/1‑149 (MQ=255)
                                                cTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGgcg                                                                                                 >  4:99943/1‑149 (MQ=255)
                                                         ggTTGATCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAAt                                                                                        <  7:219959/149‑1 (MQ=255)
                                                          gTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAAt                                                                                        >  1:155909/1‑148 (MQ=255)
                                                          gTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATg                                                                                       <  8:287910/149‑1 (MQ=255)
                                                                cTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTg                                                                                  <  5:7875/148‑1 (MQ=255)
                                                                cTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTg                                                                                  >  6:7875/1‑148 (MQ=255)
                                                                 tGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAggg                                                                                                    >  5:27829/1‑129 (MQ=255)
                                                                 tGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAggg                                                                                                    <  6:27829/129‑1 (MQ=255)
                                                                 tGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGt                                                                                >  4:28811/1‑149 (MQ=255)
                                                                  gCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTg                                                                               >  2:293873/1‑149 (MQ=255)
                                                                  gCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTg                                                                               >  4:191828/1‑149 (MQ=255)
                                                                     ggATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCtt                                                                            <  5:126364/149‑1 (MQ=255)
                                                                         gCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCCCACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTc                                                                        >  4:32603/1‑149 (MQ=255)
                                                                           tCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCACGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACTAACACATCACAGGCGCAGGGGCGCCTGCCAATGAACTcttgcttgttatg                                                                      >  4:234442/1‑144 (MQ=255)
                                                                            cGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTCCTTGGTCTg                                                                      >  7:76335/1‑148 (MQ=255)
                                                                             gTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTgga                                                                    <  3:191828/149‑1 (MQ=255)
                                                                             gTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTgg                                                                     <  1:110711/148‑1 (MQ=255)
                                                                              tGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCATCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTgga                                                                    <  3:234442/148‑1 (MQ=255)
                                                                                ccGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTggagga                                                                 >  1:31519/1‑149 (MQ=255)
                                                                                 cGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGAc                                                                <  3:99943/149‑1 (MQ=255)
                                                                                     aTCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGc                                                            >  7:193564/1‑149 (MQ=255)
                                                                                       cccGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCATCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCa                                                          >  4:261720/1‑149 (MQ=255)
                                                                                                  tACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGAt                                               >  4:260833/1‑149 (MQ=255)
                                                                                                   aCCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATc                                              >  2:288417/1‑149 (MQ=255)
                                                                                                    ccGCTATCCGCCATCACCCCGTAGACTGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCg                                             <  8:258787/149‑1 (MQ=255)
                                                                                                    ccGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGCTCg                                             >  8:225931/1‑149 (MQ=255)
                                                                                                     cGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCgg                                            >  8:171621/1‑149 (MQ=255)
                                                                                                          tCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCct                                         <  1:293873/147‑1 (MQ=255)
                                                                                                             gCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTgg                                    >  8:113476/1‑149 (MQ=255)
                                                                                                             gCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTgg                                    <  8:76335/149‑1 (MQ=255)
                                                                                                               cATCACCCCGTAGAATGCGCTCCTTTTTTCATCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTGGCa                                  <  3:261720/149‑1 (MQ=255)
                                                                                                                       cGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTgga                                                                    <  5:78928/107‑1 (MQ=255)
                                                                                                                       cGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTgga                                                                    >  6:78928/1‑107 (MQ=255)
                                                                                                                        gTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTGGCAGCCCTgcag                        >  5:89459/1‑149 (MQ=255)
                                                                                                                         tAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTGGCAGCCCTgcagc                       <  3:28811/149‑1 (MQ=255)
                                                                                                                                cgcTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTGGCAGCCCTGCAGCACctgct                >  3:109459/1‑149 (MQ=255)
                                                                                                                                 gcTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTGGCAGCCCTGCAGCACctgctg               >  4:39130/1‑149 (MQ=255)
                                                                                                                                   tCCTTTTTTCTTCCGGTGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTgg                                    <  2:47146/127‑1 (MQ=255)
                                                                                                                                   tCCTTTTTTCTTCCGGTGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTgg                                    >  1:47146/1‑127 (MQ=255)
                                                                                                                                   tCCTTTTTTCTTCCGGTGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTgg                                    >  1:51708/1‑127 (MQ=255)
                                                                                                                                   tCCTTTTTTCTTCCGGTGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTgg                                    <  2:51708/127‑1 (MQ=255)
                                                                                                                                   tCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTGGCAGCCCTGCAGCACCTGCTGGc             >  3:141966/1‑149 (MQ=255)
                                                                                                                                   tCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTGGCAGCCCTGCAGCACCTGCTGGc             >  2:164129/1‑149 (MQ=255)
                                                                                                                                   tCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTGGCAGCCCTGCAGCACCTGCTGGc             >  1:121636/1‑149 (MQ=255)
                                                                                                                                     cTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGtt                                                  <  2:301388/111‑1 (MQ=255)
                                                                                                                                     cTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGtt                                                  >  1:301388/1‑111 (MQ=255)
                                                                                                                                     cTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTGGCAGCCCTGCAGCACCTGCTGGCg            <  7:225931/148‑1 (MQ=255)
                                                                                                                                      ttttttCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTGGCAGCCCTGCAGCACCTGCTGGCGat          >  8:260901/1‑149 (MQ=255)
                                                                                                                                      ttttttCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCAAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTGGCAGCCCTGCAGCACCTGCTGGCGa           >  5:292395/1‑148 (MQ=255)
                                                                                                                                           tcttcCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTGGCAGCCCTGCAGCACCTGCTGGCGatcatc      <  8:193564/148‑1 (MQ=255)
                                                                                                                                               ccGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTGGCAGCCCTGCAGCACCTGCTGGCGATCATCGTGc  <  2:31519/148‑1 (MQ=255)
                                                                                                                                                cGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTGGCAGCCCTCGCAGCACCTGCTGGCGATCATCGTGc  <  1:164129/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACACGAACACAGCCCAGGCGCAGGGGCGCCTGCCAATGAACTGGTGCTTGGTCTGGAGGACAAGCCACGGCTGTTGATCGGCCTGCTGGCAGCCCTGCAGCACCTGCTGGCGATCATCGTGC  >  NC_002947/4502211‑4502501

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: