Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,356,500 T→A 6.2% R116R (CGA→CGT)  PP_2952 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,356,5000TA 6.2% 192.9 / 2.6 64R116R (CGA→CGT) PP_2952LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base A (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAGGGCC  >  NC_002947/3356354‑3356644
                                                                                                                                                  |                                                                                                                                                
aCCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCg                                                                                                                                                <  5:311235/149‑1 (MQ=255)
 ccTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCgg                                                                                                                                               <  3:99666/149‑1 (MQ=255)
  cTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGc                                                                                                                                              >  1:367749/1‑149 (MQ=255)
        gTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACg                                                                                                                                        >  4:56773/1‑149 (MQ=255)
        gTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACg                                                                                                                                        >  5:45382/1‑149 (MQ=255)
           ccAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGcc                                                                                                                                      <  2:147017/148‑1 (MQ=255)
            cAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGt                                                                                                                                    <  8:250974/149‑1 (MQ=255)
              aGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGa                                                                                                                                  <  3:56773/149‑1 (MQ=255)
                      gTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAggg                                                                                                                           >  2:318615/1‑148 (MQ=255)
                           aGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGgcagcag                                                                                                                     >  5:175759/1‑149 (MQ=255)
                           aGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGgcagcag                                                                                                                     >  5:7275/1‑149 (MQ=255)
                           aGGGCTATATAGTTGGCGTCGGCGTTCAGGCAATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGgcagcag                                                                                                                     >  1:230550/1‑149 (MQ=255)
                                  tataGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCgg                                                                                                                                               >  4:16138/1‑116 (MQ=255)
                                  tataGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCgg                                                                                                                                               <  3:16138/116‑1 (MQ=255)
                                  tataGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTg                                                                                                              >  1:7378/1‑149 (MQ=255)
                                  tataGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTg                                                                                                              <  4:373442/149‑1 (MQ=255)
                                      gTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTgcgc                                                                                                           <  2:103004/148‑1 (MQ=255)
                                         ggCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGcccc                                                                                                        <  2:271272/148‑1 (MQ=255)
                                          gCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGccc                                                                                                         >  1:249333/1‑146 (MQ=255)
                                                gCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGa                                                                                                <  1:44114/149‑1 (MQ=255)
                                                       ggCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGcgcg                                                                                         <  1:190814/149‑1 (MQ=255)
                                                        gCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGccc                                                                                                         <  2:94469/132‑1 (MQ=255)
                                                        gCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGccc                                                                                                         >  1:94469/1‑132 (MQ=255)
                                                          cATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGcc                                                                                      >  1:317871/1‑149 (MQ=255)
                                                          cATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGc                                                                                       >  1:35527/1‑148 (MQ=255)
                                                            tAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCAc                                                                                    <  3:285277/149‑1 (MQ=255)
                                                                  aTGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCaa                                                                              <  2:249333/149‑1 (MQ=255)
                                                                        atatCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTcc                                                                        <  2:230550/149‑1 (MQ=255)
                                                                           tCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCt                                                                     >  5:326948/1‑149 (MQ=255)
                                                                             tACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGtt                                                                                                               <  4:301902/105‑1 (MQ=255)
                                                                             tACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGtt                                                                                                               >  3:301902/1‑105 (MQ=255)
                                                                             tACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCt                                                                     >  4:217406/1‑147 (MQ=255)
                                                                               ccGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCt                                                                     <  1:318615/145‑1 (MQ=255)
                                                                               ccGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCt                                                                     <  2:317871/145‑1 (MQ=255)
                                                                               ccGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAg                                                                  <  2:352052/148‑1 (MQ=255)
                                                                                  gATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTgg                                                              >  1:295775/1‑149 (MQ=255)
                                                                                     cgcgAAAAAGATCGCCGATCTGGTCGGAGAGAAACAAGCGCACGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGCAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGgcg                                                           >  3:224375/1‑149 (MQ=255)
                                                                                          aaaaGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTc                                                      >  1:35174/1‑149 (MQ=255)
                                                                                          aaaaGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTc                                                      <  2:35174/149‑1 (MQ=255)
                                                                                               tttcgCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCt                                                    <  5:234778/142‑1 (MQ=255)
                                                                                                tCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACg                                                >  8:204391/1‑149 (MQ=255)
                                                                                                tCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACg                                                <  2:185772/149‑1 (MQ=255)
                                                                                                 cGGCGATCTGGTCGGAGAGAAACAAGCGCTGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACg                                                >  6:89685/1‑148 (MQ=255)
                                                                                                    cGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAgg                                            <  2:7378/149‑1 (MQ=255)
                                                                                                         tGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTc                                       >  4:340350/1‑149 (MQ=255)
                                                                                                           gTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTgg                                                              >  4:202961/1‑124 (MQ=255)
                                                                                                           gTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTc                                       <  8:278009/147‑1 (MQ=255)
                                                                                                           gTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTc                                       <  2:367749/147‑1 (MQ=255)
                                                                                                             cGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCgcgc                                   >  4:16348/1‑149 (MQ=255)
                                                                                                              ggAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGcc                                  >  4:156303/1‑149 (MQ=255)
                                                                                                                   gaAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGc                             <  6:45382/149‑1 (MQ=255)
                                                                                                                     caCAAGCGCAGGTTCAACTGTGGTTGCTGTCGGCGAAACGCCGTGAGCCAGGGGAGCTGCACGTTGCGCCCCATATCGGAGGGCGCGGCGACCTGCAAGGTTCCGCTGAGCTGGGCGCTTTCGTGACGCAGGTTGTCGCGACCGTGCCt                           <  4:224375/148‑1 (MQ=255)
                                                                                                                      aCAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGTGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGc                             <  6:122401/146‑1 (MQ=255)
                                                                                                                      aCAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCt                           >  3:42302/1‑148 (MQ=255)
                                                                                                                      aCAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGAGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGc                             >  5:122401/1‑146 (MQ=255)
                                                                                                                        aaGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAg                        <  8:153609/149‑1 (MQ=255)
                                                                                                                         aGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAgg                       <  2:295775/149‑1 (MQ=255)
                                                                                                                                gTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCa                <  7:204391/149‑1 (MQ=255)
                                                                                                                                       tgtgGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGAt         >  5:49203/1‑149 (MQ=255)
                                                                                                                                            aTGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAggg    <  6:326948/149‑1 (MQ=255)
                                                                                                                                              gCTCACGGCGTTTCGCCGAGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAgg                       <  3:61497/128‑1 (MQ=255)
                                                                                                                                              gCTCACGGCGTTTCGCCGAGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAgg                       >  4:61497/1‑128 (MQ=255)
                                                                                                                                              gCTCACGGCGTTTCGCCGAGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATc        <  6:309445/143‑1 (MQ=255)
                                                                                                                                              gCTCACGGCGTTTCGCCGAGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATc        >  5:309445/1‑143 (MQ=255)
                                                                                                                                                tctcGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAGGGcc  >  7:225572/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAGGGCC  >  NC_002947/3356354‑3356644

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: ‑