Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,948,139 G→A 50.0% A112V (GCC→GTC)  rho ← transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,948,1390GA50.0% ‑5.2 / 89.1 64A112V (GCC→GTC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base G (15/17);  new base A (17/15);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

ACGCGACCGGTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAG  >  NC_002947/5947992‑5948281
                                                                                                                                                   |                                                                                                                                               
aCGCGACCGGTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGc                                                                                                                                                <  3:178414/149‑1 (MQ=255)
   cgACCGGTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGaa                                                                                                                                             <  3:14430/149‑1 (MQ=255)
         gTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGaa                                                                                                                                             <  4:304168/143‑1 (MQ=255)
         gTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGaa                                                                                                                                             >  3:304168/1‑143 (MQ=255)
            aaGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTaa                                                                                                                                         <  6:303169/144‑1 (MQ=255)
            aaGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTaa                                                                                                                                         >  5:303169/1‑144 (MQ=255)
            aaGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTc                                                                                                                                    <  3:321611/149‑1 (MQ=255)
              gTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTccc                                                                                                                                  >  2:124794/1‑149 (MQ=255)
               tCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTcccc                                                                                                                                 <  3:162269/149‑1 (MQ=255)
                 cTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCcctt                                                                                                                               >  8:210790/1‑149 (MQ=255)
                 cTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGGTTTACGGGCGGTCGAAGTTGATGGTATCGACCTGAAGCAGGGCGAAGTAACGCAcccaat                                                                                                                               >  8:256422/1‑146 (MQ=255)
                    ggTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCccttc                                                                                                                              <  7:256422/147‑1 (MQ=255)
                     gTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAg                                                                                                                                            <  5:374737/132‑1 (MQ=255)
                     gTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAg                                                                                                                                            >  6:374737/1‑132 (MQ=255)
                        gAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCg                                                                                                                         <  5:66951/148‑1 (MQ=255)
                          gCCGTTACCGGCCTCCATCTTCAGGCTCTTGTTCGGGAACAGCTGCGTCAGGGTTTCGAACAGGATCTTGATCTTCGCGTTTTCCTGGCGGTCGAAGTTGATGCTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTg                                                                                                                      <  7:278126/149‑1 (MQ=255)
                                  cGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATc                                                                                                              >  4:177645/1‑149 (MQ=255)
                                  cGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATc                                                                                                              <  4:273631/149‑1 (MQ=255)
                                       tCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGcc                                                                                                         >  2:96340/1‑149 (MQ=255)
                                       tCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGcc                                                                                                         >  7:307687/1‑149 (MQ=255)
                                          aTCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTa                                                                                                                                          <  1:82783/113‑1 (MQ=255)
                                          aTCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTa                                                                                                                                          >  2:82783/1‑113 (MQ=255)
                                               cAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATg                                                                                                  <  1:124794/148‑1 (MQ=255)
                                                 ggCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACg                                                                                                                                       >  6:28306/1‑109 (MQ=255)
                                                             ggAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACg                                                                                    >  3:265576/1‑148 (MQ=255)
                                                              gAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTATTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACg                                                                                    >  1:365929/1‑147 (MQ=255)
                                                                       cgTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAAc                                                                         <  3:279163/149‑1 (MQ=255)
                                                                        gTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAAcg                                                                        >  5:20212/1‑149 (MQ=255)
                                                                            ggTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAAcggcg                                                                     >  3:151743/1‑148 (MQ=255)
                                                                                         gATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGc                                                       >  7:176447/1‑148 (MQ=255)
                                                                                               gTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACa                                                 >  2:233007/1‑148 (MQ=255)
                                                                                               gTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACa                                                 >  2:142986/1‑148 (MQ=255)
                                                                                                   ttcGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCgag                                                    >  7:306965/1‑141 (MQ=255)
                                                                                                   ttcGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCgag                                                    <  8:306965/141‑1 (MQ=255)
                                                                                                     cgcgTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGat                                            <  1:233007/147‑1 (MQ=255)
                                                                                                     cgcgTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGata                                           <  8:307687/148‑1 (MQ=255)
                                                                                                           ttCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCt                                     <  7:210790/148‑1 (MQ=255)
                                                                                                                ggCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGggccg                               <  1:142986/149‑1 (MQ=255)
                                                                                                                ggCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGggccg                               <  4:151743/149‑1 (MQ=255)
                                                                                                                    gTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGggccggc                             >  4:252609/1‑147 (MQ=255)
                                                                                                                     tCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACa                                                 <  5:283350/126‑1 (MQ=255)
                                                                                                                     tCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACa                                                 >  6:283350/1‑126 (MQ=255)
                                                                                                                      cGAAGTTGATGGTATCGACCTGCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAgg                         <  4:199579/149‑1 (MQ=255)
                                                                                                                       gAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGt                        <  8:176447/149‑1 (MQ=255)
                                                                                                                         aGTTGATGGTGTCGCCCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTGGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAg                      >  4:167077/1‑149 (MQ=255)
                                                                                                                             gATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGggccggc                             <  6:90347/138‑1 (MQ=255)
                                                                                                                             gATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGggccggc                             >  5:90347/1‑138 (MQ=255)
                                                                                                                             gATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACg                  >  5:155239/1‑149 (MQ=255)
                                                                                                                             gATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACg                  >  5:84736/1‑149 (MQ=255)
                                                                                                                             gATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGACGGCCAGGTAGTACg                  >  2:274840/1‑149 (MQ=255)
                                                                                                                              atcgTATCGACATTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTAGTTGCGCGGATCTTGCCGACCATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGAACTGGCTGGGCGAGACATAGATATAGTCTTGTCCGCCCAGGAAGGACg                  >  4:295964/4‑148 (MQ=255)
                                                                                                                                ggTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAg                                              >  4:265289/1‑118 (MQ=255)
                                                                                                                                ggTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAg                                              <  3:265289/118‑1 (MQ=255)
                                                                                                                                 gTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATc              <  3:252609/149‑1 (MQ=255)
                                                                                                                                  tATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCAt               >  3:225671/1‑147 (MQ=255)
                                                                                                                                   aTCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGGCGAGGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAg            >  2:156172/1‑149 (MQ=255)
                                                                                                                                      gACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGcc          <  3:177645/148‑1 (MQ=255)
                                                                                                                                      gACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGt                                       <  6:357077/119‑1 (MQ=255)
                                                                                                                                      gACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGt                                       >  5:357077/1‑119 (MQ=255)
                                                                                                                                        ccTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGcc          <  2:365929/146‑1 (MQ=255)
                                                                                                                                         cTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCTAAc      <  3:167077/149‑1 (MQ=255)
                                                                                                                                              gcagGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCCTAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAg  <  1:96340/149‑1 (MQ=255)
                                                                                                                                               cagGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAg  >  3:68365/1‑147 (MQ=255)
                                                                                                                                               cagGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAg  <  4:68365/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ACGCGACCGGTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAG  >  NC_002947/5947992‑5948281

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: