Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,352,174 A→G 8.0% intergenic (+64/+78) 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,1740AG 8.0% 173.9 / 4.9 62intergenic (+64/+78)PP_2947/PP_2948transcriptional regulator MvaT/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (30/27);  new base G (2/3);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGA  >  NC_002947/3352032‑3352322
                                                                                                                                              |                                                                                                                                                    
cAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCg                                                                                                                                                >  7:197320/1‑149 (MQ=255)
          aaCCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGa                                                                                                                                       >  6:11836/1‑148 (MQ=255)
            ccACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAGGCGACGCTGTCGCTTCCTCGaaaa                                                                                                                                    <  8:220522/149‑2 (MQ=255)
             cacaAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGaaa                                                                                                                                     >  3:328031/1‑147 (MQ=255)
                aaGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCTGAAGAAGTGGAATCCTGGGTGTCCTGTTTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAggg                                                                                                                                >  6:367674/1‑149 (MQ=255)
                      ctCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCa                                                                                                                          >  6:283041/1‑149 (MQ=255)
                       tCAAGGAATGGAAAGCCCAGTACGGGGCAGNAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCAt                                                                                                                         >  1:28629/1‑149 (MQ=255)
                                       ccAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAAt                                                                                                         <  3:243414/149‑1 (MQ=255)
                                           tACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGc                                                                                                     <  5:240704/149‑1 (MQ=255)
                                              ggggCAGAAGAAGGGGACTCCTGGGGGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGa                                                                                                  <  8:197320/149‑1 (MQ=255)
                                                ggCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACg                                                                                                >  2:119878/1‑149 (MQ=255)
                                                 gCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGt                                                                                               >  7:71568/1‑149 (MQ=255)
                                                   agaagaAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTgg                                                                                             <  4:328031/149‑1 (MQ=255)
                                                    gaagaaGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGt                                                                                            >  5:167459/1‑149 (MQ=255)
                                                     aagaagTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGt                                                                                            >  4:220149/1‑148 (MQ=255)
                                                         agTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTaaa                                                                                        <  5:283041/148‑1 (MQ=255)
                                                            ggAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAg                                                                                     >  5:96026/1‑148 (MQ=255)
                                                             gAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgt                                                                                                                      <  4:229071/114‑1 (MQ=255)
                                                             gAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgt                                                                                                                      >  3:229071/1‑114 (MQ=255)
                                                             gAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGa                                                                                   <  5:11836/149‑1 (MQ=255)
                                                                         tCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCAc                                                                        <  2:28629/148‑1 (MQ=255)
                                                                         tCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCAc                                                                        <  3:220149/148‑1 (MQ=255)
                                                                          ccTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACt                                                                       >  2:5083/1‑148 (MQ=255)
                                                                           cTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACt                                                                       >  5:96353/1‑147 (MQ=255)
                                                                             gATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAgc                                                                   >  8:150669/1‑149 (MQ=255)
                                                                               tttGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGGGGAAt                                                                                                         >  1:215704/1‑109 (MQ=255)
                                                                               tttGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAAt                                                                                                         <  2:215704/109‑1 (MQ=255)
                                                                                   tgCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAATCGCTCa                                                              <  1:81429/148‑1 (MQ=255)
                                                                                   tgCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCa                                                              <  6:96353/148‑1 (MQ=255)
                                                                                       gATGTCGCCCTCATTGCGCTATCTGAAAGGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTc                                                         >  1:201203/1‑149 (MQ=255)
                                                                                       gATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTc                                                         >  7:269199/1‑149 (MQ=255)
                                                                                       gATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAACGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTc                                                         >  1:98091/1‑149 (MQ=255)
                                                                                        aTGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCa                                                              >  6:264297/1‑143 (MQ=255)
                                                                                        aTGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCa                                                              <  5:264297/143‑1 (MQ=255)
                                                                                         tGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTATAGAGGAAGCGACAGCTTCGCCTCCTCGAAACTGGGCGTTCATCGTGGGTGTTCGGAATGAGCTTACGTGGTTAATCAGGAAGTTTCGCCACTAAGCGCTCAGTGTcc                                                        <  2:298425/148‑1 (MQ=255)
                                                                                             gCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGttt                                                    >  1:30629/1‑148 (MQ=255)
                                                                                             gCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTc                                                   <  6:96026/149‑1 (MQ=255)
                                                                                              cccTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTc                                                   >  1:119877/1‑148 (MQ=255)
                                                                                               ccTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCaa                                                 >  7:24515/1‑149 (MQ=255)
                                                                                                cTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGc                                                                                                              <  2:351088/87‑1 (MQ=255)
                                                                                                cTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGc                                                                                                              >  1:351088/1‑87 (MQ=255)
                                                                                                cTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCaa                                                 >  5:223932/1‑148 (MQ=255)
                                                                                                        gcTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGt                                         >  4:251704/1‑148 (MQ=255)
                                                                                                                aaaCGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGaa                                                                                  <  4:110249/99‑1 (MQ=255)
                                                                                                                aaaCGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGaa                                                                                  >  3:110249/1‑99 (MQ=255)
                                                                                                                 aaCGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCCTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGc                               <  5:367674/149‑1 (MQ=255)
                                                                                                                   cGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCAt                             <  6:223932/149‑1 (MQ=255)
                                                                                                                   cGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCAt                             <  2:119877/149‑1 (MQ=255)
                                                                                                                               tttCGAGGAGGCGACGCTGTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCa                                                              <  7:264649/104‑1 (MQ=38)
                                                                                                                               tttCGAGGAGGCGACGCTGTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCa                                                              >  8:264649/1‑104 (MQ=38)
                                                                                                                               tttCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTgg                  <  8:269199/148‑1 (MQ=255)
                                                                                                                               tttCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTgg                  <  1:5083/148‑1 (MQ=255)
                                                                                                                                ttCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTggg                 <  8:24515/148‑1 (MQ=255)
                                                                                                                                ttCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGc                >  6:210661/1‑149 (MQ=255)
                                                                                                                                  cGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGACACTAAGCGCTCAGTGTCCGTTTCAATACATCGTAAATTTAGGCATCAGCCGATTGGGCGa              >  5:151258/1‑149 (MQ=255)
                                                                                                                                    aggaggCGACGCTGTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTc                                                         <  7:113012/104‑1 (MQ=38)
                                                                                                                                    aggaggCGACGCTGTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTc                                                         >  8:113012/1‑104 (MQ=38)
                                                                                                                                    aGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACa            <  2:30629/149‑1 (MQ=255)
                                                                                                                                          cGACAGCGTCGCCTCCTCTAAACAGGGCGTTCATCGTGGGTCTGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGaaa      <  1:119878/149‑1 (MQ=255)
                                                                                                                                          cGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGaaa      <  6:167459/149‑1 (MQ=255)
                                                                                                                                           gACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAAc     >  8:262615/1‑149 (MQ=255)
                                                                                                                                             cAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACgg   <  7:150669/149‑1 (MQ=255)
                                                                                                                                              aGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACgga  <  2:98091/149‑1 (MQ=255)
                                                                                                                                              aGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACgga  <  2:201203/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGA  >  NC_002947/3352032‑3352322

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