Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,352,182 C→T 9.5% intergenic (+72/+70) PP_2947 → / ← PP_2948 transcriptional regulator MvaT/GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,352,1820CT 9.5% 170.9 / 7.1 63intergenic (+72/+70)PP_2947/PP_2948transcriptional regulator MvaT/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (31/26);  new base T (2/4);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

AGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGAGGAAATCA  >  NC_002947/3352034‑3352330
                                                                                                                                                     |                                                                                                                                                     
aGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGcc                                                                                                                                                       >  7:193027/1‑149 (MQ=255)
    tGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGcctcct                                                                                                                                                   <  8:167675/149‑1 (MQ=255)
       cAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGa                                                                                                                                                <  5:239973/149‑1 (MQ=255)
        aaCCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAGGCGACGCTGTCGCTTCCCCGaa                                                                                                                                               <  5:118520/149‑1 (MQ=255)
        aaCCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGa                                                                                                                                                >  2:220428/1‑148 (MQ=255)
          ccACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAGGCGACGCTGTCGCTTCCTCGaaaa                                                                                                                                             <  8:132982/149‑2 (MQ=255)
           cacaAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGaaa                                                                                                                                              >  5:146805/1‑147 (MQ=255)
            acaAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACa                                                                                                                                            <  1:63703/149‑1 (MQ=255)
               aGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGc                                                                                                                                        >  5:350075/1‑149 (MQ=255)
                        aGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATctt                                                                                                                              >  5:240013/1‑147 (MQ=255)
                          gAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGTGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGATACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgtgg                                                                                                                            <  8:91558/149‑1 (MQ=255)
                           aaTGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgtggg                                                                                                                           <  1:144534/149‑1 (MQ=255)
                               gAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGtgggtggg                                                                                                                       >  7:251552/1‑149 (MQ=255)
                               gAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGtgggtggg                                                                                                                       >  7:31570/1‑149 (MQ=255)
                                 aaGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCg                                                                                                                     <  6:217145/149‑1 (MQ=255)
                                      cAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAAt                                                                                                                 >  6:354943/1‑148 (MQ=255)
                                             gggCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGAc                                                                                                         <  4:34961/149‑1 (MQ=255)
                                               gCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGt                                                                                                       <  4:19663/149‑1 (MQ=255)
                                               gCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGt                                                                                                       >  1:368740/1‑149 (MQ=255)
                                                 agaagaAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGc                                                                                                             >  1:249391/1‑141 (MQ=255)
                                                 agaagaAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGc                                                                                                             <  2:249391/141‑1 (MQ=255)
                                                     gaagTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTaa                                                                                                 <  7:290016/149‑1 (MQ=255)
                                                      aagTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGAGGAATGAGCTGACGTGGTTaaa                                                                                                >  1:374984/1‑149 (MQ=255)
                                                         tGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAg                                                                                             >  7:32293/1‑149 (MQ=255)
                                                           gAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGa                                                                                           <  7:70059/149‑1 (MQ=255)
                                                                       tCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACt                                                                               <  8:251552/149‑1 (MQ=255)
                                                                          tGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAg                                                                            >  6:66969/1‑149 (MQ=255)
                                                                                 tgCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACCGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCa                                                                      <  6:350075/148‑1 (MQ=255)
                                                                                       tGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCg                                                               <  2:33663/149‑1 (MQ=255)
                                                                                           gCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGttt                                                            >  3:63783/1‑148 (MQ=255)
                                                                                           gCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGttt                                                            >  5:69301/1‑148 (MQ=255)
                                                                                           gCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGttt                                                            >  1:352025/1‑148 (MQ=255)
                                                                                            cccTCATTGCGCTATCTGAAACGACTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGTATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCa                                                          <  1:354190/149‑1 (MQ=255)
                                                                                             ccTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACa                                                                                              >  7:345171/1‑112 (MQ=255)
                                                                                             ccTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACa                                                                                              <  8:345171/112‑1 (MQ=255)
                                                                                                 aTTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACt                                                     >  3:70416/1‑149 (MQ=255)
                                                                                                    gcgcTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCg                                                  >  6:76490/1‑149 (MQ=255)
                                                                                                     cgcTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCg                                                  >  2:19105/1‑148 (MQ=255)
                                                                                                      gcTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGt                                                 >  7:356604/1‑148 (MQ=255)
                                                                                                         aTCTGAAACGCCTTGTTCTTTTTCGAGGAAGAGACAGCGTCGCCTCCTCGAAACGGGGCGTTCATCGTGGGTGTGCGGTATGAGAGGACGTGGTTATACAGTATGTTTCGGCACTAAGCGGTCAGTGTACGTTGGAATACttgatatat                                             >  8:341893/1‑142 (MQ=255)
                                                                                                            tGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAAtttt                                          >  2:100798/1‑149 (MQ=255)
                                                                                                              aaaCGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTgg                                        <  6:146805/149‑1 (MQ=255)
                                                                                                              aaaCGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTg                                         <  4:63783/148‑1 (MQ=255)
                                                                                                                   ccTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgtgg                                                                                                                            <  4:120698/60‑1 (MQ=255)
                                                                                                                   ccTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgtgg                                                                                                                            >  3:120698/1‑60 (MQ=255)
                                                                                                                   ccTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCa                                   <  2:120906/149‑1 (MQ=255)
                                                                                                                     ttGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGc                                 <  2:374984/149‑1 (MQ=255)
                                                                                                                            tttttCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTg                           >  4:366347/1‑149 (MQ=255)
                                                                                                                            tgtttCGAGGAGGCGACGCTGTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTaa                                                                             >  8:351680/3‑99 (MQ=37)
                                                                                                                            tgtttCGAGGAGGCGACGCTGTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTaa                                                                             <  7:351680/97‑1 (MQ=37)
                                                                                                                             ttttCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGCTGGGCGGAATGAGCTGAAGTGGTTAAAAAGGAAGTTTAGCCAATAAGCGCTCAGTGTCCGTTGAAATACTTCGTAAATTTTGGCATCAGCCGATTgg                          >  8:240880/1‑149 (MQ=255)
                                                                                                                              tttCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTg                           <  8:32293/147‑1 (MQ=255)
                                                                                                                                tCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGc                        <  5:76490/148‑1 (MQ=255)
                                                                                                                                tCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGc                        <  5:354943/148‑1 (MQ=255)
                                                                                                                                  gaggagGCGACGCTGTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGa                      >  5:127100/1‑148 (MQ=255)
                                                                                                                                  gAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGa                      >  2:262188/1‑148 (MQ=255)
                                                                                                                                  gAGGAAGCGACAGCGTCGCCTACTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACACGAAGTTTCGAAACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTCGCATCAGCCGATTGGGCGa                      >  5:137391/1‑148 (MQ=255)
                                                                                                                                      aaGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTg                 <  7:280939/149‑1 (MQ=255)
                                                                                                                                          gACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAAc             >  3:5156/1‑149 (MQ=255)
                                                                                                                                           aCAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACg            >  4:379000/1‑149 (MQ=255)
                                                                                                                                             gctgTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACgga          <  6:127100/146‑1 (MQ=255)
                                                                                                                                              gCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACggag         <  1:220428/149‑1 (MQ=255)
                                                                                                                                                gTCTCCTCCTCTAATCAGGGCGTTCATCGTCTGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAATCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTTGCATCAGCCGATTGTGCGACATTGAAACggagg        <  7:240880/148‑1 (MQ=255)
                                                                                                                                                 tCGCCTCCTCGAAACAGGGCGTTCATGCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACggagg        >  1:203081/1‑148 (MQ=255)
                                                                                                                                                 tCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACggagga       >  3:341175/1‑148 (MQ=255)
                                                                                                                                                    cctcctCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGAGGAAATc   <  8:31570/149‑1 (MQ=255)
                                                                                                                                                    cctcctCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGAGGAAATc   <  1:298987/149‑1 (MQ=255)
                                                                                                                                                     ctcctcGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGAGGAAATCa  >  5:101882/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGAGGAAATCA  >  NC_002947/3352034‑3352330

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: