Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,127,616 T→C 8.6% Q90R (CAG→CGG)  PP_5375 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,127,6160TC 8.6% 159.9 / 4.8 59Q90R (CAG→CGG) PP_5375LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (28/25);  new base C (2/3);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

GCGTAATCCTTCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCAAGAAGGCGAAC  >  NC_002947/6127472‑6127761
                                                                                                                                                  |                                                                                                                                                          
gCGTAATCCTTCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGc                                                                                                                                                         <  8:28348/148‑1 (MQ=255)
gCGTAATCCTTCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGc                                                                                                                                                         <  2:214820/148‑1 (MQ=255)
      tCCTTCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATc                                                                                                                                                 <  8:104709/149‑1 (MQ=255)
              aTGCTCGTGGGCGATGATCCGCAGGTCGAGGTTTTCATGGGCACTCTCGAATTCCGACACCTTGGGCAGTAAGCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTgg                                                                                                                                  <  7:206795/149‑1 (MQ=255)
              aTGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTgg                                                                                                                                  >  4:262484/1‑149 (MQ=255)
              aTGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACCGGCGATATACGC‑‑‑‑‑‑‑CAGTgg                                                                                                                                  <  5:299584/149‑1 (MQ=255)
                 aTCGTGGGCGATGATCCGCAGGGCGATGTTTTCATGGGCACTCTGTAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACCGGCGATATACGC‑‑‑‑‑‑‑CAGTGgcaa                                                                                                                               <  6:140135/149‑3 (MQ=255)
                  tCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGCGACT‑‑TGTTGGTTGCCTTGCCACCGGCGATATACGC‑‑‑‑‑‑‑CAGTGgcaag                                                                                                                              <  1:117419/149‑4 (MQ=255)
                              atCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGt                                                                                                                                                       >  1:89925/1‑120 (MQ=255)
                              atCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGt                                                                                                                                                       <  2:89925/120‑1 (MQ=255)
                               tCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAgc                                                                                                                <  6:41517/149‑1 (MQ=255)
                                ccGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAgca                                                                                                               >  3:137547/1‑149 (MQ=255)
                                  gCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGa                                                                                                                    >  4:277547/1‑143 (MQ=255)
                                  gCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGa                                                                                                                    <  3:277547/143‑1 (MQ=255)
                                         gATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTg                                                                                                                                   >  5:59598/1‑121 (MQ=255)
                                         gATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTg                                                                                                                                   <  6:59598/121‑1 (MQ=255)
                                             ttttCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAATGCATGCTTCAcctg                                                                                                   <  3:103191/148‑4 (MQ=255)
                                              tttCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCggg                                                                                                                            >  1:289221/1‑123 (MQ=255)
                                              tttCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCggg                                                                                                                            <  2:289221/123‑1 (MQ=255)
                                              tttCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCAcc                                                                                                 >  2:112071/1‑149 (MQ=255)
                                                tCATGGGCACTCTGGAATTCCGCCACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGgc                                                                                                                                 >  4:67107/1‑116 (MQ=255)
                                                tCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTTTATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGt                                                                                               >  3:58632/1‑149 (MQ=255)
                                                tCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTg                                                                                                                                   >  4:43868/1‑114 (MQ=255)
                                                tCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGgc                                                                                                                                 >  6:185047/1‑116 (MQ=255)
                                                          tctGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCAcc                                                                                                 >  7:116330/1‑137 (MQ=255)
                                                          tctGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCAcc                                                                                                 <  8:116330/137‑1 (MQ=255)
                                                          tctGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACGGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCa                                                                                                       <  5:192143/131‑1 (MQ=255)
                                                          tctGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACCGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCa                                                                                                       >  6:192143/1‑131 (MQ=255)
                                                          tctGGAATTCCGACACCGTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATa                                                                                      >  4:37657/1‑148 (MQ=255)
                                                            tGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACTGAGGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATAc                                                                                     >  1:278197/1‑149 (MQ=255)
                                                                        acCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTTTATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTggg                                                                       <  4:58632/149‑1 (MQ=255)
                                                                          cTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCg                                                                                           >  8:252591/1‑127 (MQ=255)
                                                                           ttGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTgggg                                                                      <  8:12556/147‑1 (MQ=255)
                                                                           ttGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTgggg                                                                      >  7:12556/1‑147 (MQ=255)
                                                                                cAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAgg                                                                >  6:196495/1‑148 (MQ=255)
                                                                                       cacaGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGc                                                         >  4:88836/1‑148 (MQ=255)
                                                                                        acaGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTa                                                       <  3:37657/149‑1 (MQ=255)
                                                                                           gCGACGCCATGGGGGTTTTGGGGGCAACGGTGTTT‑‑TTTTTGTTTCCTTGCCATTGGGTTATTTCCCCCTTTTTCGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACAt                                                           <  1:112071/149‑1 (MQ=255)
                                                                                             gACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAggg                                                               <  2:64723/136‑1 (MQ=255)
                                                                                             gACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGggt                                                  <  5:196495/149‑1 (MQ=255)
                                                                                                  cATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGACACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTAACGCACCGTTTAGTAATACTTCGCGCAGGTGGGGGTCAGGGACATGCTACTGGTGGTAc                                             >  1:276515/1‑149 (MQ=255)
                                                                                                         gTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGaa                                      >  1:65911/1‑149 (MQ=255)
                                                                                                            cTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATa                                    <  3:88836/148‑1 (MQ=255)
                                                                                                               gtggCAACGGTGACT‑‑TGTTGGATGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCAtt                                                                                                         >  5:226399/1‑76 (MQ=255)
                                                                                                               gtggCAACGGTGACT‑‑TGTTGGATGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCAtt                                                                                                         <  6:226399/76‑1 (MQ=255)
                                                                                                                tggCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGcc                               <  2:278197/149‑1 (MQ=255)
                                                                                                                  gCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTg                                 <  7:331607/145‑1 (MQ=255)
                                                                                                                  gCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTg                                 >  8:331607/1‑145 (MQ=255)
                                                                                                                  gCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCtt                             >  8:257263/1‑149 (MQ=255)
                                                                                                                    aaCGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGaa                                      >  4:63896/1‑138 (MQ=255)
                                                                                                                    aaCGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGaa                                      <  3:63896/138‑1 (MQ=255)
                                                                                                                      cGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTg                                 <  4:91311/141‑1 (MQ=255)
                                                                                                                      cGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTg                                 >  3:91311/1‑141 (MQ=255)
                                                                                                                        gTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGAc                                                             <  2:226578/111‑1 (MQ=255)
                                                                                                                        gTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGAc                                                             >  1:226578/1‑111 (MQ=255)
                                                                                                                             t‑‑tgttgGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGggtg                                                 <  2:11989/118‑1 (MQ=255)
                                                                                                                             t‑‑tgttgGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGggtg                                                 >  1:11989/1‑118 (MQ=255)
                                                                                                                                           gcgccACCGGCGATATACGC‑‑‑‑‑‑‑CAGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTc                 >  5:20571/3‑138 (MQ=255)
                                                                                                                                               cACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCAAGAAGGCGAAc  >  2:264214/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                          
GCGTAATCCTTCACATGATCGTGGGCGATGATCCGCAGGTCGATGTTTTCATGGGCACTCTGGAATTCCGACACCTTGGGCAGTAACCACAGCGACGCCATGGCGGTACTGGTGGCAACGGTGACT‑‑TGTTGGTTGCCTTGCCACTGGCGTATATCGC‑‑‑‑‑‑‑CGGTGGCGCGGGCCAGGTGAAGCATGCATTCACGCACCGTTTCGTAATACTGCGCGCCGGTGGGGGTCAGGGACATGCTACGGGTGGTACGGGTGAATAGTGCCTTGCCCAGGTAGTCCTCAAGAAGGCGAAC  >  NC_002947/6127472‑6127761

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: