Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 726,621 C→T 8.9% E45K (GAA→AAA)  PP_0620 ← GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947726,6210CT 8.9% 189.6 / 7.0 67E45K (GAA→AAA) PP_0620GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (23/38);  new base T (3/3);  total (26/41)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.60e-01

AGTGCAGGTAGTTGCGCAGTGCGCGGCTGGCTGGTTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCTTCTGGCGCTTGCGTTTGGGTTTGAGCATTTCGCTATTTTTCA  >  NC_002947/726474‑726767
                                                                                                                                                   |                                                                                                                                                  
aGTGCAGGTAGTTGCGCAGTGCGCGGCTGGCTGGTGCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCg                                                                                                                                                   <  3:104661/149‑1 (MQ=255)
    cAGGTAGTTGCGCAGTGCGCGGCTGGCTGGTTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAAt                                                                                                                                               <  7:65910/149‑1 (MQ=255)
        tAGTTGCGCAGTGCGCGGCTGGCTGGTTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTc                                                                                                                                           <  2:122344/149‑1 (MQ=255)
           ttGCGCAGTGCGCGGCTGGCTGGTTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTTCGAATGCTCCaa                                                                                                                                        <  2:325652/149‑2 (MQ=255)
            tGCGCAGTGCGCGGCTGGCTGGTTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATc                                                                                                                                       <  8:13473/149‑1 (MQ=255)
             gcgcAGTGCGCGGCTGGCTGGTTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATc                                                                                                                                       <  8:21167/148‑1 (MQ=255)
              cgcAGTGCGCGGCTGGCTGGTTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCa                                                                                                                                      <  4:259914/148‑1 (MQ=255)
              cgcAGTGCGCGGCTGGCTGGTTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCa                                                                                                                                      <  8:296305/148‑1 (MQ=255)
               gcAGTGCGCGGCTGGCTGGTTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGc                                                                                                                                    >  4:142261/1‑149 (MQ=255)
                    gcgcGGCTGGCTGGTTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCg                                                                                                                               >  5:88336/1‑149 (MQ=255)
                        ggctggctGGGTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCg                                                                                                                           <  7:122700/149‑1 (MQ=255)
                          ctggctggTTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGt                                                                                                                         <  7:291322/149‑1 (MQ=255)
                             gctggTTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTagg                                                                                                                      >  8:55617/1‑149 (MQ=255)
                              ctggTTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTaggc                                                                                                                     <  4:66542/149‑1 (MQ=255)
                               tggTTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCAt                                                                                                                                        >  4:28096/1‑129 (MQ=255)
                               tggTTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCAt                                                                                                                                        <  3:28096/129‑1 (MQ=255)
                                  ttCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGc                                                                                                                                             <  6:188380/121‑1 (MQ=255)
                                     ggTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGt                                                                                                              >  1:122708/1‑149 (MQ=255)
                                      gTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAACTGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTc                                                                                                             >  4:286990/1‑149 (MQ=255)
                                          ccACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGACGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTaggcag                                                                                                                   <  3:233191/139‑1 (MQ=255)
                                          ccACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGACGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTaggcag                                                                                                                   >  4:233191/1‑139 (MQ=255)
                                           cACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCgg                                                                                                        >  6:286277/1‑149 (MQ=255)
                                           cACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCgg                                                                                                        >  7:78889/1‑149 (MQ=255)
                                            aCTTCGTTCAGGTAGGCGTCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCggg                                                                                                       <  2:311532/149‑1 (MQ=255)
                                                cGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTg                                                                                                    <  4:71223/148‑1 (MQ=255)
                                                 gTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCa                                                                                                  <  3:260668/149‑1 (MQ=255)
                                                  ttCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAg                                                                                                 >  8:216846/1‑149 (MQ=255)
                                                   tCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGG‑‑GGTGCCCTTGGAGCATTCg                                                                                                                                                   <  4:199757/96‑1 (MQ=255)
                                                   tCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGG‑‑GGTGCCCTTGGAGCATTCg                                                                                                                                                   >  3:199757/1‑96 (MQ=255)
                                                    cAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAgc                                                                                                <  1:58034/148‑1 (MQ=255)
                                                       gTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCtt                                                                                            >  4:271828/1‑149 (MQ=255)
                                                        tAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTc                                                                                           <  7:81173/149‑1 (MQ=255)
                                                               ccccGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGAc                                                                                     >  1:163653/1‑148 (MQ=255)
                                                                 ccGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGAt                                                                                  <  5:286277/149‑1 (MQ=255)
                                                                    gTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGa                                                                               <  6:34493/149‑1 (MQ=255)
                                                                        ggCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAgcgc                                                                           <  7:165676/149‑1 (MQ=255)
                                                                          ccGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCtt                                                                         <  6:88336/149‑1 (MQ=255)
                                                                           cGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTg                                                                        <  5:17125/149‑1 (MQ=255)
                                                                                 gggTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGAt                                                                  <  8:78889/149‑1 (MQ=255)
                                                                                  ggTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGAtt                                                                 <  2:163653/149‑1 (MQ=255)
                                                                                  ggTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGAtt                                                                 <  7:55617/149‑1 (MQ=255)
                                                                                  ggTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGAtt                                                                 >  3:320968/1‑149 (MQ=255)
                                                                                   gTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTTCGAATGCTCCaag                                                                                                                                       <  5:271012/78‑3 (MQ=38)
                                                                                   gTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTTCGAATGCTCCaag                                                                                                                                       >  6:271012/1‑76 (MQ=38)
                                                                                     gCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTcga                                                              <  8:81067/149‑1 (MQ=255)
                                                                                      cTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTcga                                                              <  4:320968/148‑1 (MQ=255)
                                                                                         aCCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGAt                                                          >  3:26129/1‑149 (MQ=255)
                                                                                           cAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATca                                                        >  5:161366/1‑149 (MQ=255)
                                                                                                       ttcgaGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCt                                            >  5:272165/1‑149 (MQ=255)
                                                                                                       ttcgaGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCt                                            >  3:247573/1‑149 (MQ=255)
                                                                                                       ttcgaGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGAGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACa                                               >  7:161669/1‑146 (MQ=255)
                                                                                                       ttcgaGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGAGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACa                                               <  8:161669/146‑1 (MQ=255)
                                                                                                           aGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCTCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCTTCTg                                        <  4:158019/149‑1 (MQ=255)
                                                                                                               ccTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCTTCTGgcgc                                    <  3:271828/149‑1 (MQ=255)
                                                                                                                         ttCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCTTCTGGCGCTTGCGTTTgg                          >  2:44866/1‑149 (MQ=255)
                                                                                                                           cgcgGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCTTCTGGCGCTTGCGTTTgg                          <  2:314729/147‑1 (MQ=255)
                                                                                                                            gcgggcggTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCTTCTGGCGCTTGCGTTTGGGtt                       <  7:331130/149‑1 (MQ=255)
                                                                                                                               ggcggTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCTTCTGGCGCTTGCGTTTGGGTTTGa                    <  2:122708/149‑1 (MQ=255)
                                                                                                                               ggcggTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCTTCTGGCGCTTGCGTTTGGGTTTGa                    <  7:216846/149‑1 (MQ=255)
                                                                                                                                    tGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCTTCTGGCGCTTGCGTTTGGGTTTGAGCAtt               >  4:216118/1‑149 (MQ=255)
                                                                                                                                      cccTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCTTCTGGCGCTTGCGTTTGGGTTTGAGCAttt              <  3:216118/148‑1 (MQ=255)
                                                                                                                                       tgatGGAGCATTTCGAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCTTCTGGCGCTTGCGTTTGGGTTTGAGCATTTCg            >  2:179025/4‑149 (MQ=255)
                                                                                                                                           ggAGCATTTCGAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCTTCTGGCGCTTGCGTTTGGGTTTGAGCAttt              <  8:120012/143‑1 (MQ=255)
                                                                                                                                           ggAGCATTTCGAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCTTCTGGCGCTTGCGTTTGGGTTTGAGCAttt              >  7:120012/1‑143 (MQ=255)
                                                                                                                                                 ttCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTTACTGATTCGACGATCACATCGGACAGCTTCTGGCGCTTGCGTTTGGGTTTGAGCATTTCGCTATTTTTCa  <  3:142261/149‑1 (MQ=255)
                                                                                                                                                 ttCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTCACTGATTCGACGATCACATCGGACAGCTTCTGGCGCTTGCGTTTGGGTTTGAGCATTTCGCTATTTTTCa  >  7:244464/1‑149 (MQ=255)
                                                                                                                                                  tCGTAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCTTCTGGCGCTTGCGTTTGGGTTTGAGCATTTCGCTATTTTTCa  >  6:248160/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGTGCAGGTAGTTGCGCAGTGCGCGGCTGGCTGGTTCGGTGCCCACTTCGTTCAGGTAGGCGCCCCCGGTCGGGCCGGTGCGGGTGCTCACCAGGCCTTCGACTTCGAGCGCCTTGAGGGCTTCGCGGGCGGTGCCCTTGGAGCATTCGAAATGCTCCATCAGCTCCCGCTCGGTAGGCAGGCGGTCGCCGGGGTGCAGCGCTTCGGTGACGATCGAGCGCTTGACTGATTCGACGATCACATCGGACAGCTTCTGGCGCTTGCGTTTGGGTTTGAGCATTTCGCTATTTTTCA  >  NC_002947/726474‑726767

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