Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 650,848 A→T 9.6% intergenic (‑137/+43) 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,8480AT 9.6% 227.9 / 9.6 80intergenic (‑137/+43)aroQ‑I/dsbD‑Itype II 3‑dehydroquinate dehydratase/thiol/disulfide interchange protein 1
Reads supporting (aligned to +/- strand):  ref base A (37/35);  new base T (2/6);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCGCAG  >  NC_002947/650703‑650994
                                                                                                                                                 |                                                                                                                                                    
cGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATc                                                                                                                                                   <  1:335046/149‑1 (MQ=255)
cGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATc                                                                                                                                                   <  7:105206/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                                                                                                                                               <  6:48659/148‑1 (MQ=255)
     ccATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCGGt                                                                                                                                               <  2:432546/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                                                                            <  1:2490/149‑1 (MQ=255)
                                                                     aaCGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGt                                                                            >  4:82994/1‑149 (MQ=255)
                                                                           aaGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAAt                                                                      <  8:264681/149‑1 (MQ=255)
                                                                               gtCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCGGTCAc                                                                                                                                            <  3:338282/77‑1 (MQ=37)
                                                                               gtCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCGGTCAc                                                                                                                                            >  4:338282/1‑77 (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                                  >  6:383055/1‑149 (MQ=255)
                                                                                                               cGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCAc                                  >  1:416635/1‑149 (MQ=255)
                                                                                                               cGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCAc                                  >  8:79276/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)
                                                                                                                                                aaaTCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCgcag  >  2:358259/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTA‑‑TAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCGCAG  >  NC_002947/650703‑650994

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: