Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 650,845 T→A 9.0% intergenic (‑134/+46) aroQ‑I ← / ← dsbD‑I type II 3‑dehydroquinate dehydratase/thiol/disulfide interchange protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947650,8450TA 9.0% 217.9 / 10.1 78intergenic (‑134/+46)aroQ‑I/dsbD‑Itype II 3‑dehydroquinate dehydratase/thiol/disulfide interchange protein 1
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base A (2/5);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

AGTAGCGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCGC  >  NC_002947/650698‑650992
                                                                                                                                                   |                                                                                                                                                     
aGTAGCGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGATCGGATAttt                                                                                                                                                      <  1:253095/149‑1 (MQ=255)
 gTAGCGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATAttt                                                                                                                                                      >  6:115979/1‑148 (MQ=255)
 gTAGCGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCACGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATAttt                                                                                                                                                      >  1:287045/1‑148 (MQ=255)
     cGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATc                                                                                                                                                 <  7:105206/149‑1 (MQ=255)
     cGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATc                                                                                                                                                 <  1:335046/149‑1 (MQ=255)
       ttGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCg                                                                                                                                               <  5:115979/149‑1 (MQ=255)
        ttccATCTGCGTATTCCTTGGATTTGTATGGCAATTCGACACAGCTCGGCGAGTGAGCCGTATAGCAACGGCAAGTGTCCAGTTCACGGCAATAGCCAGCACGATGTCCGAGGTTCGAGCGAAGTATTTGACCTGATATTTAAATCCg                                                                                                                                               <  4:298174/146‑1 (MQ=255)
         gCCATCTGCGGATTCCTTGGATTTGTAGGGCATTTCGACTCAGCGCGGCGAGTGTGGCGTAAAGCTACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCgg                                                                                                                                              <  5:232296/148‑1 (MQ=255)
         gCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCgg                                                                                                                                              <  2:398298/148‑1 (MQ=255)
         gCCATCTGCGGATTCCTTGGATTGGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCgg                                                                                                                                              <  8:188319/148‑1 (MQ=255)
          ccATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCGGt                                                                                                                                             <  2:432546/148‑1 (MQ=255)
          ccATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCGGt                                                                                                                                             <  6:48659/148‑1 (MQ=255)
           cATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCa                                                                                                                                           >  6:3651/1‑149 (MQ=255)
                  ggATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCAttt                                                                                                                                    <  5:3651/149‑1 (MQ=255)
                    aTTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCa                                                                                                                                  <  6:235892/149‑1 (MQ=255)
                     ttCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCAGTCACTCATTTCAg                                                                                                                                 >  1:409492/1‑149 (MQ=255)
                        cTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGAt                                                                                                                              >  6:381154/1‑149 (MQ=255)
                         ttGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGAtt                                                                                                                             >  1:309855/1‑149 (MQ=255)
                           ggATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGAtttt                                                                                                                           >  4:430428/1‑149 (MQ=255)
                            gATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGAttttt                                                                                                                          >  8:74740/1‑149 (MQ=255)
                              tttGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTcc                                                                                                                        >  6:35945/1‑149 (MQ=255)
                                tGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCg                                                                                                                       <  1:230540/148‑1 (MQ=255)
                                    gggCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTc                                                                                                                  >  3:74484/1‑149 (MQ=255)
                                     ggCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGACCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTcg                                                                                                                 >  1:138055/1‑149 (MQ=255)
                                        aaTTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATCGCCAGCACGATGTCCGAGGGTCGCGCGAAGTATGTGACCCGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCgg                                                                                                              <  3:156946/149‑1 (MQ=255)
                                             gACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGt                                                                                                         >  1:221896/1‑149 (MQ=255)
                                                 caGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTta‑‑ta                                                                                                   <  2:309855/149‑1 (MQ=255)
                                                  aGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCATTAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCgg                                                                                                              >  4:113656/1‑139 (MQ=255)
                                                  aGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCgg                                                                                                              <  3:113656/139‑1 (MQ=255)
                                                   gcgcGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGa                                                                                                 >  3:117626/1‑149 (MQ=255)
                                                      cggcgAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCg                                                                                              <  5:381154/149‑1 (MQ=255)
                                                        gcgAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTACCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGAc                                                                                            >  7:264681/1‑149 (MQ=255)
                                                               tgCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCcg                                                                                     <  5:431372/149‑1 (MQ=255)
                                                                gCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCcgc                                                                                    >  4:52298/1‑149 (MQ=255)
                                                                  cGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCcgcg                                                                                   >  1:190192/1‑148 (MQ=255)
                                                                     aaaGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCAcg                                                                               >  3:192784/1‑149 (MQ=255)
                                                                         cAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCAcgcgcg                                                                           >  4:168574/1‑149 (MQ=255)
                                                                          aaCGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGt                                                                          >  4:82994/1‑149 (MQ=255)
                                                                          aaCGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGt                                                                          <  1:2490/149‑1 (MQ=255)
                                                                                aaGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAAt                                                                    <  8:264681/149‑1 (MQ=255)
                                                                                    gtCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCGGTCAc                                                                                                                                          >  4:338282/1‑77 (MQ=37)
                                                                                    gtCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCGGTCAc                                                                                                                                          <  3:338282/77‑1 (MQ=37)
                                                                                     tCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGaa                                                                                                           >  1:250442/1‑107 (MQ=255)
                                                                                     tCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGaa                                                                                                           <  2:250442/107‑1 (MQ=255)
                                                                                     tCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTgcg                                                               <  2:190192/149‑1 (MQ=255)
                                                                                       cAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTgcgc                                                              <  2:7481/148‑1 (MQ=255)
                                                                                       cAGCTCCCGGCAATATCCAGCACGATGTCCGAGGTTCGCTCGAAGTATGTGACCGGATATTTATATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGATTCTgagc                                                              <  2:287045/148‑1 (MQ=255)
                                                                                         gTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGaa                                                           <  1:49445/149‑1 (MQ=255)
                                                                                            cccGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGATGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAAtt                                                         >  5:84927/1‑148 (MQ=255)
                                                                                            cccGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAAtt                                                         >  4:29585/1‑148 (MQ=255)
                                                                                            cccGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAAtt                                                         <  2:138055/148‑1 (MQ=255)
                                                                                             ccGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCa                                                       >  6:15741/1‑149 (MQ=255)
                                                                                              cGGCAATAGCCAGCACGATGTCCGATGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCAc                                                      <  4:12412/149‑1 (MQ=255)
                                                                                                gCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCg                                                    <  2:221896/149‑1 (MQ=255)
                                                                                                  aaTAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCgg                                                   >  6:19670/1‑148 (MQ=255)
                                                                                                       ccAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTcgc                                                                                                                <  8:328931/84‑1 (MQ=255)
                                                                                                        cAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTAt                                            <  3:82994/149‑1 (MQ=255)
                                                                                                         aGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAAtt                                                         >  6:360074/1‑135 (MQ=255)
                                                                                                         aGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAAtt                                                         <  5:360074/135‑1 (MQ=255)
                                                                                                         aGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTAtc                                           >  4:195135/1‑149 (MQ=255)
                                                                                                          gCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGGTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTAtct                                          <  1:254511/149‑1 (MQ=255)
                                                                                                           cACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTAtctc                                         >  6:78648/1‑149 (MQ=255)
                                                                                                             cGATGTCCGAGGTGCGCGCGAAGTTTGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATTTAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTAtctc                                         <  1:126153/149‑1 (MQ=255)
                                                                                                                  tCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAAt                                   <  7:74740/148‑1 (MQ=255)
                                                                                                                    cGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTT‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCAAATCAc                                >  7:268692/1‑149 (MQ=255)
                                                                                                                    cGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCAc                                >  8:79276/1‑149 (MQ=255)
                                                                                                                    cGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCAc                                >  6:383055/1‑149 (MQ=255)
                                                                                                                    cGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCAc                                >  1:416635/1‑149 (MQ=255)
                                                                                                                         ttCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGa                            >  1:370147/1‑148 (MQ=255)
                                                                                                                           cgcgcgAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATc                          <  5:78648/148‑1 (MQ=255)
                                                                                                                             ctcgAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGc                       <  5:35945/147‑1 (MQ=255)
                                                                                                                                      tgtgACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAAttt               <  6:84927/148‑1 (MQ=255)
                                                                                                                                      tgtgACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTc              >  3:278668/1‑149 (MQ=255)
                                                                                                                                        tgACCGGATTTAAATATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGc            >  1:272268/1‑149 (MQ=255)
                                                                                                                                        tgACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGc            >  3:191160/1‑149 (MQ=255)
                                                                                                                                        tgACCGGATATTTAAAACCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGc            <  3:430428/149‑1 (MQ=255)
                                                                                                                                         gACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGc            <  4:278668/148‑1 (MQ=255)
                                                                                                                                         gACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCt           >  8:360280/1‑149 (MQ=255)
                                                                                                                                          aCCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTg          >  3:122457/1‑149 (MQ=255)
                                                                                                                                            cGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGcc        >  2:178545/1‑149 (MQ=255)
                                                                                                                                                 atTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCg   >  1:263541/1‑149 (MQ=255)
                                                                                                                                                  tttAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCgc  <  1:175225/149‑1 (MQ=255)
                                                                                                                                                  tttAAATCCAGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCgc  <  2:409492/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGTAGCGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCGC  >  NC_002947/650698‑650992

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: