Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,965,245 T→C 9.4% D38G (GAC→GGC)  fleQ ← transcriptional regulator FleQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,965,2450TC 9.4% 176.0 / 7.2 64D38G (GAC→GGC) fleQtranscriptional regulator FleQ
Reads supporting (aligned to +/- strand):  ref base T (34/24);  new base C (2/4);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

ACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTCCCTAGTCATATCGGCGCCCTG  >  NC_002947/4965100‑4965392
                                                                                                                                                 |                                                                                                                                                    
aCCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCAt                                                                                                                                                   <  3:321607/149‑1 (MQ=255)
      aaGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGCCAtgactgg                                                                                                                                             >  8:13105/1‑149 (MQ=255)
         cAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGCCATGACTGGcag                                                                                                                                          <  5:301000/149‑2 (MQ=255)
         cAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGCCATGACTGGcag                                                                                                                                          <  7:13105/149‑2 (MQ=255)
            aaCCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGCCATGACTGGcagcag                                                                                                                                       <  5:129165/149‑5 (MQ=255)
                ggaaggaaCTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAg                                                                                                                                  <  2:117996/149‑1 (MQ=255)
                   aggaaCTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGtt                                                                                                                               <  6:117705/149‑1 (MQ=255)
                     gaaCTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGtttt                                                                                                                             >  6:135317/1‑149 (MQ=255)
                             cccAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGcc                                                                                                                      >  7:173529/1‑148 (MQ=255)
                              ccAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGc                                                                                                                       <  8:174297/146‑1 (MQ=255)
                              ccAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGc                                                                                                                       >  7:174297/1‑146 (MQ=255)
                                     gCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGtt                                                                                                             <  6:175453/149‑1 (MQ=255)
                                        aCTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAg                                                                                                          >  1:162009/1‑149 (MQ=255)
                                              ttAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGcaccac                                                                                                    >  2:370292/1‑149 (MQ=255)
                                                    agcCCAAGCAGATTGCCCGGAGCATTGACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAg                                                                                              <  2:205136/149‑1 (MQ=255)
                                                    agcCCAAGCAGATAGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAg                                                                                              <  3:90242/149‑1 (MQ=255)
                                                        cAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTcg                                                                                          >  6:377047/1‑149 (MQ=255)
                                                         aaGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTcgc                                                                                         >  7:299254/1‑149 (MQ=255)
                                                         aaGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTcgc                                                                                         >  8:200375/1‑149 (MQ=255)
                                                              gATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCgcggcgg                                                                                    <  2:277247/149‑1 (MQ=255)
                                                               aTTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCgcggcggc                                                                                   >  6:80876/1‑149 (MQ=255)
                                                                    ccGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGgcgcgcg                                                                              >  8:36543/1‑149 (MQ=255)
                                                                    ccGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGgcgcgcg                                                                              >  3:363767/1‑149 (MQ=255)
                                                                         gCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTgtc                                                                         <  7:110924/149‑1 (MQ=255)
                                                                         gCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTgtc                                                                         <  5:377047/149‑1 (MQ=255)
                                                                           aTTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACACGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTgtcgt                                                                       >  4:103221/1‑149 (MQ=255)
                                                                            ttCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCgcgg                                                                                       <  4:125895/132‑1 (MQ=255)
                                                                            ttCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCgcgg                                                                                       >  3:125895/1‑132 (MQ=255)
                                                                            ttCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTgtcgt                                                                       >  1:341638/1‑148 (MQ=255)
                                                                            ttCACGGTGCCGATGAGCACGCACAGCACTTCGCGACTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTgtcgt                                                                       >  7:303293/1‑148 (MQ=255)
                                                                              cACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCAt                                                                    >  5:279606/1‑149 (MQ=255)
                                                                                     cctaTGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCTCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATcag                                                             <  3:127336/146‑1 (MQ=255)
                                                                                     ccGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATcag                                                             >  7:264451/1‑149 (MQ=255)
                                                                                     ccGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATca                                                              <  5:135317/148‑1 (MQ=255)
                                                                                       gATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGAGGCGGCGCGCGCTGTCGTCATCGATcagca                                                           <  8:173529/149‑1 (MQ=255)
                                                                                        aTGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATcagcag                                                          <  4:363767/149‑1 (MQ=255)
                                                                                            gcacgcacAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATc                                                      <  2:162009/149‑1 (MQ=255)
                                                                                            gcacgcacAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATc                                                      >  6:10420/1‑149 (MQ=255)
                                                                                               cgcacAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTg                                                   >  5:282316/1‑149 (MQ=255)
                                                                                                gcacAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTgg                                                  <  8:299254/149‑1 (MQ=255)
                                                                                                 cacaGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCgcggcggc                                                                                   <  3:268014/115‑1 (MQ=255)
                                                                                                 cacaGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCgcggcggc                                                                                   >  4:268014/1‑115 (MQ=255)
                                                                                                  acaGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGtt                                                >  8:119474/1‑149 (MQ=255)
                                                                                                                   gAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCaa                                 >  7:198815/1‑147 (MQ=255)
                                                                                                                   gAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCaa                                 <  8:303293/147‑1 (MQ=255)
                                                                                                                   gAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATa                               <  7:36543/149‑1 (MQ=255)
                                                                                                                        cAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGGGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAAc                          <  5:20914/149‑1 (MQ=255)
                                                                                                                         aaTGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACt                         >  3:214063/1‑149 (MQ=255)
                                                                                                                             gCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTccc                     >  7:118297/1‑149 (MQ=255)
                                                                                                                               tCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCGAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTccc                     >  1:307083/1‑147 (MQ=255)
                                                                                                                                cGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTCCCTAg                  <  8:264451/149‑1 (MQ=255)
                                                                                                                                cGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTCCCTAg                  >  3:358470/1‑149 (MQ=255)
                                                                                                                                    cGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTCCCTAGTCat              >  5:207498/1‑149 (MQ=255)
                                                                                                                                       ctgctgCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACGACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTCCCTAGTCATATc           >  1:102937/1‑149 (MQ=255)
                                                                                                                                        gtgtgtcAGCCATGACTGGCAGCAGGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTCCCTa                   >  6:285211/7‑141 (MQ=255)
                                                                                                                                        gtgtgtcAGCCATGACTGGCAGCAGGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTCCCTa                   <  5:285211/135‑1 (MQ=255)
                                                                                                                                         gctgcCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTCCCTAGTCATATCgg         >  5:256560/1‑149 (MQ=255)
                                                                                                                                           tgcCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTCCCTAGTCATATCGgcg       >  4:15801/1‑149 (MQ=255)
                                                                                                                                            gcCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTCCCTAGTCATATCGgcgc      >  5:376055/1‑149 (MQ=255)
                                                                                                                                              cAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGc                             >  6:5469/1‑124 (MQ=255)
                                                                                                                                              cAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGc                             <  5:5469/124‑1 (MQ=255)
                                                                                                                                               aGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTCCCTAGTCATATCGGCGCCCt   <  2:102937/149‑1 (MQ=255)
                                                                                                                                                gTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCCCGCCACATGCAATAGCAACTTCCCTAGTCATATCGGCGCCCTg  >  4:230439/1‑149 (MQ=255)
                                                                                                                                                gTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTACCTAGTCATATCGGCGCCCTg  >  2:122200/1‑149 (MQ=255)
                                                                                                                                                 tCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTCCCTAGTCATATCGGCGCCCTg  >  3:360832/1‑148 (MQ=255)
                                                                                                                                                 tCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTCCCTAGTCATATCGGCGCCCTg  >  1:170950/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTGGTTTCACGCCACATGCAATAGCAACTTCCCTAGTCATATCGGCGCCCTG  >  NC_002947/4965100‑4965392

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: