Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,756,127 A→T 5.9% H523L (CAT→CTT)  PP_3319 → GGDEF domain‑containing sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,756,1270AT 5.9% 206.3 / 2.9 68H523L (CAT→CTT) PP_3319GGDEF domain‑containing sensory box protein
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base T (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

ATCGCGTTTGCGAGCACAGAAGCCTACGTAATGACCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGTTGCGCGACCTATGTGATCAGGGCTTCAACGGTATCTCGAT  >  NC_002947/3755986‑3756274
                                                                                                                                             |                                                                                                                                                      
aTCGCGTTTGCGAGCACAGAAGCCTACGTAATGACCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGt                                                                                                                                                  >  6:357012/1‑148 (MQ=255)
        tGCGAGCACAGAAGCCTACGTAATGACCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGtt                                                                                                                                         >  2:221996/1‑149 (MQ=255)
          cGAGCACAGAAGCCTACGTAATGACCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTc                                                                                                                                        <  4:282782/148‑1 (MQ=255)
           gAGCACAGAAGCCTACGTAATGACCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCAt                                                                                                                                      <  4:294355/149‑1 (MQ=255)
               acaGAAGCCTACGTAATGACCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCAtctct                                                                                                                                  >  8:50011/1‑149 (MQ=255)
                 aGAAGCCTACGTAATGACCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAggg                                                                                                                                           <  1:24133/138‑1 (MQ=255)
                 aGAAGCCTACGTAATGACCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAggg                                                                                                                                           >  2:24133/1‑138 (MQ=255)
                    aGCCTACGTAATGACCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAg                                                                                                                              <  7:50011/148‑1 (MQ=255)
                       cTACGTAATGACCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTag                                                                                                                                               >  4:177491/1‑128 (MQ=255)
                       cTACGTAATGACCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTag                                                                                                                                               <  3:177491/128‑1 (MQ=255)
                       cTACGTAATGACCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGt                                                                                                                          >  4:200246/1‑149 (MQ=255)
                       cTACGTAATGACCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCATGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCAGGTTAGAGGGTTCATCTCTCCAGCGGt                                                                                                                          >  6:248642/1‑149 (MQ=255)
                        tACGTAATGACCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAAACGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGtt                                                                                                                         >  5:33363/1‑149 (MQ=255)
                                 aCCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCg                                                                                                                >  8:95434/1‑149 (MQ=255)
                                 aCCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCg                                                                                                                >  7:142683/1‑149 (MQ=255)
                                  ccTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAggg                                                                                                                                           <  2:353890/121‑1 (MQ=255)
                                  ccTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAggg                                                                                                                                           >  1:353890/1‑121 (MQ=255)
                                  ccTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCgg                                                                                                               <  1:221996/149‑1 (MQ=255)
                                      cAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTg                                                                                                             <  8:108412/147‑1 (MQ=255)
                                       aaaaGGCTATCAAAACTCAATCGATCAGCCTGCATTAACAACCGCAAATACAACCGCGCACCCGTCGGGTGGTCGCGATCGAACCGTTGTTGCGCAGTCAGCCTCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCCGTGGCa                                                                                                          >  3:225363/1‑149 (MQ=255)
                                         aaGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCa                                                                                                          >  1:53426/1‑147 (MQ=255)
                                           ggCTATCAAAACTCAATCGTTCAGCCTGCTTTCTCAACCGCAACTAGTCGCGCGCTCGCGTCGGGTGGTCGGGATCGTAGCGTTGTAGCGCTGGGAGCATCCGGATAGAGGGTTCATCTCTCCAGCGGTTTTCATAGCGGTGGCAGAgg                                                                                                      <  6:69948/149‑1 (MQ=255)
                                              tATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGa                                                                                                     >  7:117830/1‑147 (MQ=255)
                                              tATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGa                                                                                                     >  5:186994/1‑147 (MQ=255)
                                                tCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGt                                                                                                 >  4:80841/1‑149 (MQ=255)
                                                  aaaaCTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTgg                                                                                               >  1:168119/1‑149 (MQ=255)
                                                  aaaaCTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTg                                                                                                <  5:357012/148‑1 (MQ=255)
                                                         aaTCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCaaa                                                                                          <  4:314982/147‑1 (MQ=255)
                                                              ttCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATccc                                                                                   <  8:117830/149‑1 (MQ=255)
                                                                     tGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCAAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAg                                                                                                  >  1:18398/1‑127 (MQ=255)
                                                                     tGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCAAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAg                                                                                                  <  2:18398/127‑1 (MQ=255)
                                                                      gCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCGGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGcc                                                                           >  8:186687/1‑149 (MQ=255)
                                                                      gCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGcc                                                                           >  2:266675/1‑149 (MQ=255)
                                                                         ttATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCttt                                                                        <  3:80841/149‑1 (MQ=255)
                                                                         ttATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCttt                                                                        <  3:200246/149‑1 (MQ=255)
                                                                         ttATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCttt                                                                        <  2:168119/149‑1 (MQ=255)
                                                                           aTCTACCGCAAATAGCAGCGCGCACGCGGCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATGCGGTTAGAGTGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTg                                                                       <  4:225363/148‑1 (MQ=255)
                                                                            tCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTggg                                                                     <  7:95434/149‑1 (MQ=255)
                                                                            tCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTgg                                                                      <  2:53426/148‑1 (MQ=255)
                                                                            tCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTgg                                                                      <  7:238638/148‑1 (MQ=255)
                                                                            tCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTgg                                                                      <  6:186994/148‑1 (MQ=255)
                                                                                 cGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCAGCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTc                                                                <  6:8972/149‑1 (MQ=255)
                                                                                  gCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGCGCCTTTGGGTACTc                                                                >  4:237938/1‑148 (MQ=255)
                                                                                   cAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCTTCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCg                                                               >  3:156606/1‑148 (MQ=255)
                                                                                   cAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCTTCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGa                                                              >  3:162033/1‑149 (MQ=255)
                                                                                          aaGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACt                                                                 >  3:353249/1‑139 (MQ=255)
                                                                                          aaGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACt                                                                 <  4:353249/139‑1 (MQ=255)
                                                                                            gcgcgcACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑Tg                                                  >  3:232399/1‑149 (MQ=255)
                                                                                                    gcgTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGt                                          >  4:240111/1‑149 (MQ=255)
                                                                                                       tCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGTTgc                                       <  1:67132/149‑1 (MQ=255)
                                                                                                          ggtggtCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGtt                                         >  3:213168/1‑144 (MQ=255)
                                                                                                          ggtggtCGGGATCGAAGCCTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGtt                                         <  4:213168/144‑1 (MQ=255)
                                                                                                           gtggtCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTg                                                                                                >  6:15580/1‑91 (MQ=255)
                                                                                                           gtggtCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTg                                                                                                <  5:15580/91‑1 (MQ=255)
                                                                                                                cGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGTTGCGCGACCTAt                              >  2:193019/1‑149 (MQ=255)
                                                                                                                   gATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATAc                                                           <  3:131108/120‑1 (MQ=255)
                                                                                                                   gATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATAc                                                           >  4:131108/1‑120 (MQ=255)
                                                                                                                        aaGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGTTGCGCGACCTATGTGATCAg                      >  8:14502/1‑149 (MQ=255)
                                                                                                                        aaGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATAGCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGTTGCGCGACCTATGTGATCAg                      <  6:33363/149‑1 (MQ=255)
                                                                                                                          gCGTTGTTGCGCTGGGAGCTTCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGTTGCGCGACCTATGTGATCAggg                    <  4:162033/149‑1 (MQ=255)
                                                                                                                          gCGTTGTTGCGCTGGGAGCTTCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGTTGCGCGACCTATGTGATCAggg                    <  4:156606/149‑1 (MQ=255)
                                                                                                                                ttgCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTTGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGTTGCGCGACCTATGTGATCAGGGCTTCa               <  7:186687/148‑1 (MQ=255)
                                                                                                                                ttgCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGTTGCGCGACCTATGTGATCAGGGCTTCa               <  1:266675/148‑1 (MQ=255)
                                                                                                                                   cgcTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGTTGCGCGACCTATGTGATCAGGGCTTCAACgg           <  3:240111/149‑1 (MQ=255)
                                                                                                                                    gcTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATACCTCTGAGCCTTTGGGTACTCGATACGGCGAAAGGCGCGCAGTTCCGCGACATATGTGATCAGGGATTCAAc             >  1:300714/1‑149 (MQ=255)
                                                                                                                                        ggAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGTTGCGCGACCTAtgtg                           <  2:206450/128‑1 (MQ=255)
                                                                                                                                        ggAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGTTGCGCGACCTAtgtg                           >  1:206450/1‑128 (MQ=255)
                                                                                                                                            cATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGTTGCGCGACCTATGTGATCAGGGCTTCAACGGTATCTCGAt  <  8:142683/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                      
ATCGCGTTTGCGAGCACAGAAGCCTACGTAATGACCTGCAAAAGGCTATCAAAACTCAATCGTTCAGCCTGCATTATCAACCGCAAATAGAAGCGCGCACGCGTCGGGTGGTCGGGATCGAAGCGTTGTTGCGCTGGGAGCATCCGGTTAGAGGGTTCATCTCTCCAGCGGTTTTCATACCGGTGGCAGAGGACAGTGGTCAAATTATCCCTCTGAGCCTTTGGGTACTCGATACGGCG‑‑‑TGCGCGCAGTTGCGCGACCTATGTGATCAGGGCTTCAACGGTATCTCGAT  >  NC_002947/3755986‑3756274

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