Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,269,389 C→G 6.5% D15H (GAT→CAT)  PP_1109 ← GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,269,3890CG 6.5% 235.2 / 7.4 77D15H (GAT→CAT) PP_1109GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (33/39);  new base G (2/3);  total (35/42)
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.16e-01

AAGTAGCCTCGCCCGGGTGATGCCGTCAACACACCTTCGCTGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTCGTTTGCAAATGGCCAGCACATGCCAA  >  NC_002947/1269241‑1269520
                                                                                                                                                     |                                                                                                                                   
aaGTAGCCTCGCCCGGGTGATGCCGTCAACACACCTTCGCTGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATc                                                                                                                                     <  8:48837/149‑1 (MQ=255)
   tAGCCTCGCCCGGGTGATGCCGTCAACACACCTTCGCTGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAg                                                                                                                                  <  1:116804/149‑1 (MQ=255)
   tAGCCTCGCCCGGGTGATGCCGTCAACACACCTTCGCTGACCAGCAGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAg                                                                                                                                  >  2:174272/1‑149 (MQ=255)
     accgcGCCCGGGTGATGCCGTCAACACACCTTCGCTGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGc                                                                                                                                <  5:43889/145‑1 (MQ=255)
      ccTCGCCCGGGTGATGCCGTCAACACACCTTCGCTGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCg                                                                                                                               >  5:240737/1‑149 (MQ=255)
          gCCCGGGTGATGCCGTCAACACACCTTCGCTGCCGAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGGCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGGATGCAGGCgag                                                                                                                             <  3:43244/147‑3 (MQ=255)
          gCCCGGGTGATGCCGTCAACACACCTTCGCTGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCgttgt                                                                                                                           >  7:30962/1‑149 (MQ=255)
           cccGGGTGATGCCGTCAACACACCTTCGCTGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGGATGCAGGCgagtgg                                                                                                                          <  8:66056/149‑6 (MQ=255)
             cGGGTGATGCCGTCAACACACCTTCGCTGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTc                                                                                                                        <  6:240737/149‑1 (MQ=255)
              gggTGATGCCGTCAACACACCTTCGCTGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTc                                                                                                                        >  4:280705/1‑148 (MQ=255)
                    tGCCGTCAACACAGCTTCGCTGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGGGGTCAGCTTGCGGATGGAGGGCATGCGTTGGCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAgg                                                                                                                  <  2:247575/148‑1 (MQ=255)
                         tCAACACACCTTCGCTGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGGATGCAGGCgagt                                                                                                                            <  4:222401/133‑4 (MQ=255)
                         tCAACACACCTTCGCTGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGGATGCAGGCgagt                                                                                                                            >  3:222401/1‑130 (MQ=255)
                              acacCTTCGCTGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCgttgt                                                                                                                           >  2:233744/1‑129 (MQ=255)
                               cacCTTCGCTGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTg                                                                                                      <  2:316996/149‑1 (MQ=255)
                                   ttCGCTGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAgg                                                                                                  >  7:216657/1‑149 (MQ=255)
                                        tGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGAt                                                                                              >  3:58018/1‑148 (MQ=255)
                                          aCCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATc                                                                                             <  4:58018/147‑1 (MQ=255)
                                           ccAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATg                                                                                          <  1:174272/149‑1 (MQ=255)
                                            cAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATgg                                                                                         >  8:231075/1‑149 (MQ=255)
                                            cAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATgg                                                                                         <  3:280705/149‑1 (MQ=255)
                                            cAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATgg                                                                                         <  1:198755/149‑1 (MQ=255)
                                             aGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGc                                                                                        <  3:294502/149‑1 (MQ=255)
                                               cTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAg                                                                                      <  5:317856/149‑1 (MQ=255)
                                                  ggCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGt                                                                                    >  1:91421/1‑149 (MQ=255)
                                                            gACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGtt                                                                         <  2:91938/149‑1 (MQ=255)
                                                            gACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGtt                                                                         <  8:30962/149‑1 (MQ=255)
                                                             acacGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTa                                                                        <  8:197213/149‑1 (MQ=255)
                                                             acacGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTa                                                                        <  6:42482/149‑1 (MQ=255)
                                                                 gacgTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGaa                                                                     >  3:200908/1‑148 (MQ=255)
                                                                 gacgTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGaa                                                                     >  4:274936/1‑148 (MQ=255)
                                                                    gtcgtGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGaaaaa                                                                  >  3:184929/1‑148 (MQ=255)
                                                                    gtcgtGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGaaaaa                                                                  >  7:234712/1‑148 (MQ=255)
                                                                       gtGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGaaaaaaaa                                                               <  3:274936/148‑1 (MQ=255)
                                                                        tGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCCTTGCCCGGGCAACCACTCCCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAACGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAg                                                              >  6:172203/1‑148 (MQ=255)
                                                                               tGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTa                                                                        >  3:313348/1‑131 (MQ=255)
                                                                               tGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTa                                                                        <  4:313348/131‑1 (MQ=255)
                                                                                 aaGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAg                                                                                                                                  >  8:456/1‑71 (MQ=255)
                                                                                 aaGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAg                                                                                                                                  <  7:456/71‑1 (MQ=255)
                                                                                   gTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTaa                                                  <  3:261525/149‑1 (MQ=255)
                                                                                   gTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTaa                                                  <  1:219511/149‑1 (MQ=255)
                                                                                    tCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAAc                                                        >  3:45375/1‑142 (MQ=255)
                                                                                    tCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAAc                                                        <  4:45375/142‑1 (MQ=255)
                                                                                    tCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAAc                                                 <  7:231075/149‑1 (MQ=255)
                                                                                         tGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAg                                                              >  2:35194/1‑131 (MQ=255)
                                                                                         tGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAg                                                              <  1:35194/131‑1 (MQ=255)
                                                                                           cAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACgg                                               <  4:301190/144‑1 (MQ=255)
                                                                                           cAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACgg                                               >  3:301190/1‑144 (MQ=255)
                                                                                              cGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCtt                                       <  1:232378/149‑1 (MQ=255)
                                                                                                 tCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTc                                    >  4:119897/1‑149 (MQ=255)
                                                                                                   aGCTTTCGGATTGATGGGATGCGTTGCCCGGGCATCCACTCTCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGGAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGa                                  <  1:193108/149‑1 (MQ=255)
                                                                                                      ttGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGa                                  >  3:58805/1‑146 (MQ=255)
                                                                                                      ttGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGa                                  >  3:63064/1‑146 (MQ=255)
                                                                                                      ttGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGa                                  <  4:58805/146‑1 (MQ=255)
                                                                                                      ttGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGa                                  <  4:63064/146‑1 (MQ=255)
                                                                                                      ttGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTg                               <  8:216657/149‑1 (MQ=255)
                                                                                                         cGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTc                            <  2:91421/149‑1 (MQ=255)
                                                                                                             aTGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTCGtt                         >  3:50196/1‑149 (MQ=255)
                                                                                                              tGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTCGttt                        <  8:234712/149‑1 (MQ=255)
                                                                                                                gAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTCGTTTg                       <  4:50196/148‑1 (MQ=255)
                                                                                                                 aGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTTGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTCGTTTGc                      >  5:301246/1‑148 (MQ=255)
                                                                                                                          gTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTc                                    <  3:270751/125‑1 (MQ=255)
                                                                                                                          gTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTc                                    >  4:270751/1‑125 (MQ=255)
                                                                                                                          gTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTCGTTTGCAAATGGCCAg            >  2:137353/1‑149 (MQ=255)
                                                                                                                             gCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGATACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTCGTTTGCAAATGGCCAGca          >  5:22404/1‑148 (MQ=255)
                                                                                                                             gCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTCGTTTGCAAATGGc               <  6:55249/143‑1 (MQ=255)
                                                                                                                             gCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTCGTTTGCAAATGGc               >  5:55249/1‑143 (MQ=255)
                                                                                                                              cccGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTCGTTTGCAAATGGCCAGcaca        >  8:258535/1‑149 (MQ=255)
                                                                                                                                 gggCTACCACTCGCCTGAATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCATGGGATCATGGCAGGTCCGAAAGCGTTAGAAAAATAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTCGTTTGCAAATGGCCAGCACATGc     <  5:64989/149‑1 (MQ=255)
                                                                                                                                 gggCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTCGTTTGCAAATGGCCAGCACATGc     <  7:258535/149‑1 (MQ=255)
                                                                                                                                  ggCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTCGTTTGCAAATGGCCa             >  6:27696/1‑140 (MQ=255)
                                                                                                                                  ggCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTCGTTTGCAAATGGCCa             <  5:27696/140‑1 (MQ=255)
                                                                                                                                    cAACCACTCGCCTGCATCCAGGCGTTGTCTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTcc                                                                                  <  4:130713/69‑1 (MQ=255)
                                                                                                                                    cAACCACTCGCCTGCATCCAGGCGTTGTCTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTcc                                                                                  >  3:130713/1‑69 (MQ=255)
                                                                                                                                     gaaacaacGCCTGGATGCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTCGTTTGCAAATGGCCAGCACATGCCaa  >  6:48181/8‑148 (MQ=255)
                                                                                                                                        cACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGaa                                                                     >  6:325239/1‑78 (MQ=255)
                                                                                                                                        cACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGaa                                                                     <  5:325239/78‑1 (MQ=255)
                                                                                                                                           tCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGaaaaaaaa                                                               >  8:298099/1‑81 (MQ=255)
                                                                                                                                                     |                                                                                                                                   
AAGTAGCCTCGCCCGGGTGATGCCGTCAACACACCTTCGCTGACCAGCTGGGCGTAGGCCGACACGACGTCGTGGTAACTGAAGTCCACTGCAGCGGTCAGCTTGCGGATGGAGGGCATGCGTTGCCCGGGCAACCACTCGCCTGCATCCAGGCGTTGTTTCAAAAGGTGATGAAGGCTGCAGGGGATCATGGCAGGTCCGAAACCGTTAGAAAAAAAAGCCTAACAGTTAACGGAAATGCTTGTCGACTGTTCGTTTGCAAATGGCCAGCACATGCCAA  >  NC_002947/1269241‑1269520

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: