Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,948,139 G→A 9.1% 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,1390GA 9.1% 181.1 / 5.3 66A112V (GCC→GTC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base A (4/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

ACGCGACCGGTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAG  >  NC_002947/5947992‑5948281
                                                                                                                                                   |                                                                                                                                              
aCGCGACCGGTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGc                                                                                                                                               <  6:284883/149‑1 (MQ=255)
  gcgACCGGTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCg                                                                                                                                              >  5:282603/1‑148 (MQ=255)
  gcgACCGGTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCg                                                                                                                                              >  3:356290/1‑148 (MQ=255)
  gcgACCGGTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACg                                                                                                                                              >  1:193071/1‑148 (MQ=255)
     aCCGGTCAAGTCCTCGGTGGAGCCGTTACCGGACTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAg                                                                                                                                           >  8:175498/1‑148 (MQ=255)
       cGGTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTaa                                                                                                                                        <  2:40657/149‑1 (MQ=255)
         gTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCATGGCGAAGTAAc                                                                                                                                       >  7:372434/1‑148 (MQ=255)
         gTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACg                                                                                                                                      <  3:221497/149‑1 (MQ=255)
          tCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGc                                                                                                                                     >  2:20541/1‑149 (MQ=255)
            aaGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTc                                                                                                                                   >  5:190303/1‑149 (MQ=255)
            aaGTCCTCGGTGGACCCCTTACCGGCCTCCAACGTCACGCGCTTGATCGGGAACAGCGGCGTCAGGTTTACGAACAGGATCATGTTCTTCGCGTTTTCCGGGCGCTCGAAGTTGATGGTATCGCCCTTCAGCAGGGCGAAGTAACGCTc                                                                                                                                   >  6:274376/1‑149 (MQ=255)
             aGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTcc                                                                                                                                  <  3:333637/149‑1 (MQ=255)
                ccTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTcccc                                                                                                                                <  3:106269/148‑1 (MQ=255)
                  ccGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAggg                                                                                                                                                <  6:60256/129‑1 (MQ=255)
                           ccGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTgg                                                                                                                    >  5:134283/1‑149 (MQ=255)
                             gTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGc                                                                                                                  <  2:209284/149‑1 (MQ=255)
                             gTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGc                                                                                                                  <  3:198575/149‑1 (MQ=255)
                              ttACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCg                                                                                                                 >  2:14681/1‑149 (MQ=255)
                                  cGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATc                                                                                                             >  2:35848/1‑149 (MQ=255)
                                  cGGCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATc                                                                                                             <  2:193071/149‑1 (MQ=255)
                                    gCCTCCATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCtt                                                                                                           <  5:350171/149‑1 (MQ=255)
                                       tgcATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGCTCTTGTTGTTCGCGTTGTCCGGGCGGTCGAAGTTGATGGTATCGGCCGTCAGCAGGGCGAAGTAGCGCTCGCCTGCCTGCGGTGGCGGGAGCTTGcc                                                                                                        >  4:305900/3‑149 (MQ=255)
                                        ccATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCcg                                                                                                       >  4:344215/1‑149 (MQ=255)
                                         cATCTTCAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCcg                                                                                                       <  6:282603/148‑1 (MQ=255)
                                               cAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGAt                                                                                                  <  4:45955/147‑1 (MQ=255)
                                               cAGGCGCTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATgg                                                                                                >  3:29900/1‑149 (MQ=255)
                                                  gcgcTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGgtgt                                                                                             >  3:388092/1‑149 (MQ=255)
                                                  gcgcTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCGTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACGTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGgtgt                                                                                             >  4:89229/1‑149 (MQ=255)
                                                     cTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTCTCGc                                                                                          >  3:177659/1‑149 (MQ=255)
                                                     cTTGTTCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTCTCGc                                                                                          <  7:117849/149‑1 (MQ=255)
                                                          tCGGGAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTa                                                                                     >  7:172553/1‑149 (MQ=255)
                                                             ggAACAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCTCATCTTGCCGACGATGGTGTCGCCGGTACGc                                                                                  <  1:35848/149‑1 (MQ=255)
                                                               aaCAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAg                                                                                <  8:172553/149‑1 (MQ=255)
                                                                  aGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTCTCGCCGGTACGCAGGtt                                                                             <  3:344215/149‑1 (MQ=255)
                                                                     ggcgTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGaa                                                                          <  4:356290/149‑1 (MQ=255)
                                                                             gTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGa                                                                  <  6:100265/149‑1 (MQ=255)
                                                                                   gAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATctgg                                                             <  1:20541/148‑1 (MQ=255)
                                                                                   gAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATctgg                                                             <  5:172893/148‑1 (MQ=255)
                                                                                   gAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTCTCGCCGGTACGCAGGTTGaaa                                                                         <  5:316227/136‑1 (MQ=255)
                                                                                   gAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTCTCGCCGGTACGCAGGTTGaaa                                                                         >  6:316227/1‑136 (MQ=255)
                                                                                             ttGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCgaga                                                   >  1:159802/1‑148 (MQ=255)
                                                                                              tGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGCAGTAATGCTCCGGTTCGTTCGCTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTTCGCAGGTTGATACGGCGGATCTGGCTGGGCGAGACa                                                 >  5:89379/1‑149 (MQ=255)
                                                                                               ggtcttcGAGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCTGCAGGGCGAAGTAAGGCTCCCCTTCCTTCTGTGGGCGTTTCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGAc                                                  <  5:274376/145‑1 (MQ=255)
                                                                                               gTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGAGCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCga                                                     <  6:190303/144‑1 (MQ=255)
                                                                                                ttcttaGCTTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCGCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATa                                               <  8:372434/149‑1 (MQ=255)
                                                                                                   ttcGCGTTTTCCGGGCGGTCGAGGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGaa                                            >  5:255012/1‑148 (MQ=255)
                                                                                                   ttcGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGa                                             >  4:193487/1‑148 (MQ=255)
                                                                                                   ttcGCGGTGTCCGGGCGGTCGATGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAg                                              <  6:255012/147‑1 (MQ=255)
                                                                                                      gcgTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATc                                         <  7:175498/149‑1 (MQ=255)
                                                                                                              cGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTgg                                   <  5:401304/147‑1 (MQ=255)
                                                                                                              cGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTgg                                   >  6:401304/1‑147 (MQ=255)
                                                                                                               gggCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGggcc                                <  4:388092/149‑1 (MQ=255)
                                                                                                                 gCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGggccgg                              >  1:45070/1‑149 (MQ=255)
                                                                                                                   ggTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGggccggcc                            <  4:29900/149‑1 (MQ=255)
                                                                                                                     tagAAGTTGATGGTATAGACCTTCAGCAGGGAGAAGTAACGCTCACCTTCCTTCTTTGGGAGGATCTTGCCGACGATGGTCTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCAGGGCGAGCAATAGATATCGTCTTGGCCGGACAg                          >  7:188573/3‑149 (MQ=255)
                                                                                                                     tCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCga                                                     <  2:227804/122‑1 (MQ=255)
                                                                                                                     tCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCga                                                     >  1:227804/1‑122 (MQ=255)
                                                                                                                     tCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAg                          >  2:150655/1‑149 (MQ=255)
                                                                                                                            tGATGTTATCGACCTTCAGCAGGTCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTTCCGACGTTGGTCTCGCCGGTACGCATTTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTTGGCCGGCCAGGTATGa                    <  8:188573/148‑1 (MQ=255)
                                                                                                                             gATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACg                  >  8:292795/1‑149 (MQ=255)
                                                                                                                               tGGTATCGCCCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAACCGGCGGATCTGGCTGGGCGAGACATAGCTATCGTCTGGGCCGGCCAGGTAGGACgga                >  6:101758/1‑147 (MQ=255)
                                                                                                                               tGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCa                <  1:150655/149‑1 (MQ=255)
                                                                                                                                   aTCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAg            >  8:92273/1‑149 (MQ=255)
                                                                                                                                      gACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGcc          <  3:89229/148‑1 (MQ=255)
                                                                                                                                           tCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACg     >  5:52887/1‑148 (MQ=255)
                                                                                                                                             agcagGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAg  <  6:134283/149‑1 (MQ=255)
                                                                                                                                             agcagGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTCTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCa   <  4:177659/148‑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: ‑