Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,974,391 T→A 6.1% N28Y (AAC→TAC)  PP_3503 ← sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,974,3910TA 6.1% 197.2 / 2.6 65N28Y (AAC→TAC) PP_3503sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base A (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GCAATCCATTCGGTGTTGCTGCGCTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTTGTGGCGCCCGTTTAAC  >  NC_002947/3974243‑3974538
                                                                                                                                                      |                                                                                                                                                   
gCAATCCATTCGGTGTTGCTGCGCTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGtt                                                                                                                                                     >  5:183524/1‑149 (MQ=255)
 cAATCCATTCGGTGTTGCTGCGCTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTg                                                                                                                                                    >  8:169478/1‑149 (MQ=255)
   aTCCATTCGGTGTTGCTGCGCGTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTgcg                                                                                                                                                  <  2:323565/149‑1 (MQ=255)
        ttCGGTGTTGCTGCGCTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGc                                                                                                                                             <  2:179303/149‑1 (MQ=255)
           ggTGTTGCTGCGCTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAg                                                                                                                                          <  5:327378/149‑1 (MQ=255)
            gtgtTGCTGCGCTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAgg                                                                                                                                         >  4:158518/1‑149 (MQ=255)
               ttGCTGCGCTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTa                                                                                                                                       <  1:239956/149‑1 (MQ=255)
                      gcTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCgg                                                                                                                               <  6:183524/149‑1 (MQ=255)
                       cTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGt                                                                                                                              >  7:155232/1‑149 (MQ=255)
                           aTCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCACGCGCAATAGCCCGACATCGCACGTGTGGTCCCGCGCAGCCCCGAGCATACAGCTGTCCACATCCCCGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGc                                                                                                                           >  2:297989/1‑148 (MQ=255)
                           aTCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGc                                                                                                                           >  2:293541/1‑148 (MQ=255)
                                gTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAAtcgt                                                                                                                     >  7:27012/1‑149 (MQ=255)
                                   cTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAAtcgtcg                                                                                                                   >  6:268869/1‑148 (MQ=255)
                                   cTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAAtcgtcg                                                                                                                   >  7:258679/1‑148 (MQ=255)
                                   cTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAAtcgtcg                                                                                                                   >  4:79123/1‑148 (MQ=255)
                                       aCCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTg                                                                                                                             >  6:194236/1‑134 (MQ=255)
                                       aCCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTg                                                                                                                             <  5:194236/134‑1 (MQ=255)
                                       aCCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAggg                                                                                                              <  8:258679/149‑1 (MQ=255)
                                         cgcgTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGt                                                                                                            >  7:75486/1‑149 (MQ=255)
                                            gTCCGGGTGCCTCAGGTGCGAAGTCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTcga                                                                                                         <  1:126975/149‑1 (MQ=255)
                                                  gTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGAt                                                                                                     >  6:100791/1‑147 (MQ=255)
                                                  gTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATaa                                                                                                   <  7:42481/149‑1 (MQ=255)
                                                    gCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAg                                                                                                  >  6:9335/1‑148 (MQ=255)
                                                       tCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATaa                                                                                                   >  8:70231/1‑144 (MQ=255)
                                                       tCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATaa                                                                                                   <  7:70231/144‑1 (MQ=255)
                                                       tCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCaaa                                                                                              >  7:193299/1‑149 (MQ=255)
                                                       tCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCaaa                                                                                              >  1:202466/1‑149 (MQ=255)
                                                         aGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACg                                                                                            >  6:212191/1‑149 (MQ=255)
                                                            tGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCg                                                                                         <  5:9335/149‑1 (MQ=255)
                                                             gCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGt                                                                                        <  5:100791/149‑1 (MQ=255)
                                                              cGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGtt                                                                                       >  7:145639/1‑149 (MQ=255)
                                                              cGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGtt                                                                                       >  1:19143/1‑149 (MQ=255)
                                                                   ccTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCCGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGGCGACGATAAGCAAACGGCGTTGAGCg                                                                                  <  1:297989/149‑1 (MQ=255)
                                                                   ccTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCg                                                                                  <  1:293541/149‑1 (MQ=255)
                                                                          gcgcAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGgct                                                                           >  5:893/1‑149 (MQ=255)
                                                                               aTAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCtgctgc                                                                      <  8:155232/149‑1 (MQ=255)
                                                                                                 tgtgATCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGCGCGGCAGGCTGCTGCATGACCGATCCTTGgcgc                                                    >  4:48658/1‑149 (MQ=255)
                                                                                                  gtgGTCCAGCGCAGCCCCGAGCATACAGCTGTCCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGgc                                                      >  1:203440/1‑147 (MQ=255)
                                                                                                  gtgGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGgc                                                      <  2:203440/146‑1 (MQ=255)
                                                                                                      tCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGtt                                               >  2:192383/1‑149 (MQ=255)
                                                                                                      tCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGtt                                               >  2:73330/1‑149 (MQ=255)
                                                                                                         aGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGttttt                                            <  3:79123/149‑1 (MQ=255)
                                                                                                          gcgcAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTg                                           >  1:119345/1‑149 (MQ=255)
                                                                                                          gcgcAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTg                                           >  8:262567/1‑149 (MQ=255)
                                                                                                                  cccGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGtttttt                                    <  1:192383/149‑1 (MQ=255)
                                                                                                                    cGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACCAGCGGTGGCAACCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGa                                  <  7:169478/149‑1 (MQ=255)
                                                                                                                     gAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGaa                                 <  6:893/149‑1 (MQ=255)
                                                                                                                            cAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTaaaa                           <  8:193299/148‑1 (MQ=255)
                                                                                                                            cAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTaaaa                           <  5:188580/148‑1 (MQ=255)
                                                                                                                               cTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGc                       <  8:75486/149‑1 (MQ=255)
                                                                                                                                 gTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGaa                                 <  6:2347/137‑1 (MQ=255)
                                                                                                                                 gTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGaa                                 >  5:2347/1‑137 (MQ=255)
                                                                                                                                     cacaTCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGtt                  >  1:47838/1‑149 (MQ=255)
                                                                                                                                     cacaTCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGtt                  >  4:313499/1‑149 (MQ=255)
                                                                                                                                       caTCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTtgt                <  2:202466/149‑1 (MQ=255)
                                                                                                                                        aTCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTtgtg               >  6:324322/1‑149 (MQ=255)
                                                                                                                                        aTCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTtgtg               <  5:268869/149‑1 (MQ=255)
                                                                                                                                          cctggCCTGCGCAACGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGaa                                 <  8:275410/125‑1 (MQ=255)
                                                                                                                                          cctggCCTGCGCAACGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGaa                                 >  7:275410/5‑129 (MQ=255)
                                                                                                                                           ctggCCTGCGCAACGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGtttt                                      >  7:113447/4‑123 (MQ=255)
                                                                                                                                           ctggCCTGCGCAACGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGtttt                                      <  8:113447/120‑1 (MQ=255)
                                                                                                                                             aGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTTGTGGCGcc          <  7:262567/149‑1 (MQ=255)
                                                                                                                                               ccTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGtttttt                                    >  2:213482/1‑121 (MQ=255)
                                                                                                                                               ccTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGtttttt                                    <  1:213482/121‑1 (MQ=255)
                                                                                                                                               ccTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTTGTGGCGCCCg        <  5:324322/149‑1 (MQ=255)
                                                                                                                                                    gTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTTGTGGCGCCCGTTTaa   >  7:79836/1‑149 (MQ=255)
                                                                                                                                                     ttGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTTGTGGCGCCCGTTTAAc  <  1:73330/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GCAATCCATTCGGTGTTGCTGCGCTTGATCATGTCCTTGACCGCGTCCGGGTGCCTCAGGTGCGAAGCCTGCAAGCGCAATAGCCCGACATCGCACGTGTGGTCCAGCGCAGCCCCGAGCATACAGCTGTCCACATCCCAGCCTGCGTTGCGCAGGCCAGGTAACAAGCGGTGGCAATCGTCGCAGGGGTCGACGATAAGCAAACGGCGTTGAGCGGCAGGCTGCTGCATGACCGATCCTTGGCGCCAGTTTTTGGTTTTTTGAAGTAAAACAGCAAGTTGTGGCGCCCGTTTAAC  >  NC_002947/3974243‑3974538

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: