Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,356,514 T→A 5.1% T112S (ACG→TCG)  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,5140TA 5.1% 242.9 / 2.3 79T112S (ACG→TCG) PP_2952LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (33/42);  new base A (2/2);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

CAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAGGGCCGAGGGCAGGTACTGCTC  >  NC_002947/3356366‑3356661
                                                                                                                                                    |                                                                                                                                                   
cAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGt                                                                                                                                                     >  8:223222/1‑149 (MQ=255)
 aaGGTGCCAGTGGCAGGGCTATATAGTAGGCGTCGGCGTTCAGGCCAAAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGACGGAGAGAAACAACCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGt                                                                                                                                                     >  2:265420/1‑148 (MQ=255)
  aGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGa                                                                                                                                                   <  1:190674/149‑1 (MQ=255)
  aGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGa                                                                                                                                                   <  7:185423/149‑1 (MQ=255)
  aGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCAGt                                                                                                                                                     >  6:310023/1‑147 (MQ=255)
             gCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGAGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGgcag                                                                                                                                         <  8:352143/148‑1 (MQ=255)
                  gCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCAc                                                                                                                                   >  1:345104/1‑149 (MQ=255)
                    tatataGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCAGGt                                                                                                                                 >  4:352390/1‑149 (MQ=255)
                      tataGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGtt                                                                                                                                <  5:216062/148‑1 (MQ=255)
                      tataGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGtt                                                                                                                                <  2:318530/148‑1 (MQ=255)
                      tataGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTg                                                                                                                               >  6:98125/1‑149 (MQ=255)
                        taGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTgcg                                                                                                                             >  5:193435/1‑149 (MQ=255)
                           ttGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGGCGTGAGCCAGGGCAGCAGCACGTTGCGccc                                                                                                                          <  1:30439/149‑1 (MQ=255)
                           ttGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGccc                                                                                                                          <  3:370546/149‑1 (MQ=255)
                           ttGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGccc                                                                                                                          >  2:154905/1‑149 (MQ=255)
                           ttGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGccc                                                                                                                          <  2:79738/149‑1 (MQ=255)
                             ggCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGccc                                                                                                                          >  3:99875/1‑147 (MQ=255)
                                gTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCTGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGccc                                                                                                                          >  6:8274/1‑144 (MQ=255)
                                  cGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGccc                                                                                                                          >  6:290686/1‑142 (MQ=255)
                                      gTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAACAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGACACGCCGTGAGCGAGGGCAGCAGCACGTTGCTCCCCAGATCGGAgg                                                                                                               <  2:169304/149‑1 (MQ=255)
                                         cAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCTGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATGGGAGGgc                                                                                                             <  6:52119/148‑1 (MQ=255)
                                         cAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGgcg                                                                                                            >  5:32421/1‑149 (MQ=255)
                                          aGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGgcgc                                                                                                           >  2:172842/1‑149 (MQ=255)
                                             ccATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCTGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCgg                                                                                                         <  4:266244/148‑1 (MQ=255)
                                               aTAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCa                                                                                                      <  7:223222/149‑1 (MQ=255)
                                                    ggATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGc                                                                                                 <  3:155997/149‑1 (MQ=255)
                                                     gATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGc                                                                                                 >  1:253745/1‑148 (MQ=255)
                                                      aTGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCaa                                                                                               <  6:65317/149‑1 (MQ=255)
                                                       tGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAAc                                                                                              <  4:167527/149‑1 (MQ=255)
                                                         gCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCCCCTGCAACGt                                                                                            >  6:214925/1‑149 (MQ=255)
                                                                 tACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCt                                                                                      >  3:225172/1‑147 (MQ=255)
                                                                 tACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCt                                                                                      >  1:80537/1‑147 (MQ=255)
                                                                 tACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGa                                                                                    <  4:99875/149‑1 (MQ=255)
                                                                    cGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCt                                                                                 <  2:345104/149‑1 (MQ=255)
                                                                     ggATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTg                                                                                >  8:52228/1‑149 (MQ=255)
                                                                     ggATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTg                                                                                >  2:111686/1‑149 (MQ=255)
                                                                       aTCGCGAAAAAGAGCGGCGATCTGGTCGGAGAGAAACCAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTggg                                                                              <  6:221218/149‑1 (MQ=255)
                                                                        tCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGgc                                                                             <  5:251054/149‑1 (MQ=255)
                                                                        tCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGgc                                                                             >  6:251054/1‑149 (MQ=255)
                                                                        tCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGgc                                                                             <  1:172842/149‑1 (MQ=255)
                                                                           cgAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGGCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCt                                                                          <  5:265327/149‑1 (MQ=255)
                                                                            gAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCtt                                                                         >  3:101450/1‑149 (MQ=255)
                                                                              aaaaGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCAGCTGCAACGTTCCGCTGAGCTGGGCGCTTTc                                                                       <  1:34643/149‑1 (MQ=255)
                                                                                 aGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTg                                                                    <  4:101450/149‑1 (MQ=255)
                                                                                   aTCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGAc                                                                  <  5:98125/149‑1 (MQ=255)
                                                                                     cGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACg                                                                 >  5:42146/1‑148 (MQ=255)
                                                                                              ggTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTcg                                                       >  2:215139/1‑149 (MQ=255)
                                                                                               gTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGt                                                            <  5:290686/143‑1 (MQ=255)
                                                                                               gTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTcgc                                                      >  5:88311/1‑149 (MQ=255)
                                                                                               gTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTc                                                        <  6:42146/147‑1 (MQ=255)
                                                                                                  ggAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCt                                                                                      <  8:208538/114‑1 (MQ=255)
                                                                                                   gagagaAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGccc                                                  >  2:242577/1‑149 (MQ=255)
                                                                                                    agagaAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCg                                                 <  2:53876/149‑1 (MQ=255)
                                                                                                          aCAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCt                                            >  8:88192/1‑148 (MQ=255)
                                                                                                          aCAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTg                                           <  7:337311/149‑1 (MQ=255)
                                                                                                            aaGCGGAGGTTCAACTGTGGAAGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCAGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAg                                         <  1:265420/149‑1 (MQ=255)
                                                                                                            aaGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAg                                         >  2:188802/1‑149 (MQ=255)
                                                                                                                   ggTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCa                                                               >  4:275689/1‑120 (MQ=255)
                                                                                                                   ggTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCa                                                               <  3:275689/120‑1 (MQ=255)
                                                                                                                       cAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAgg                               <  5:8274/148‑1 (MQ=255)
                                                                                                                            gtgGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGAt                          >  7:108406/1‑148 (MQ=255)
                                                                                                                              ggATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAg                       >  4:206764/1‑149 (MQ=255)
                                                                                                                               gATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAg                                                                                                                       <  4:297701/52‑1 (MQ=255)
                                                                                                                               gATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAg                                                                                                                       >  3:297701/1‑52 (MQ=255)
                                                                                                                               gATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAgg                      <  8:132124/149‑1 (MQ=255)
                                                                                                                               gATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAgg                      <  6:32421/149‑1 (MQ=255)
                                                                                                                               ctggCTCACGGCGTTTCGCCGAGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAgg                                        >  2:119164/4‑131 (MQ=255)
                                                                                                                               ctggCTCACGGCGTTTCGCCGAGAGCCAGGGCAGCAGCACGTTGCGCCCCAGAACGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAgg                                        <  1:119164/128‑1 (MQ=255)
                                                                                                                                tggCTCACGGCGTTTCGCCGAGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTc                                                        >  6:51274/3‑114 (MQ=255)
                                                                                                                                tggCTCACGGCGTTTCGCCGAGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTc                                                        <  5:51274/112‑1 (MQ=255)
                                                                                                                                aTGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAggg                     <  7:88192/149‑1 (MQ=255)
                                                                                                                                 ggCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAGGGc                    >  6:62428/2‑149 (MQ=255)
                                                                                                                                      tcGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAGGGCCGAgg               >  1:204373/1‑149 (MQ=255)
                                                                                                                                        ggCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAGGGCCGAgg               <  2:253745/147‑1 (MQ=255)
                                                                                                                                        ggCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAGGGCCGAg                <  1:215139/146‑1 (MQ=255)
                                                                                                                                        ggCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAGGGCCGAGGGc             <  2:87952/149‑1 (MQ=255)
                                                                                                                                         gCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAGGGCCGAggg              >  5:79007/1‑147 (MQ=255)
                                                                                                                                           gAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAGGGCCGAGGGCAgg          >  2:177834/1‑149 (MQ=255)
                                                                                                                                           gAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAGGGCCGAGGGCAgg          >  7:176306/1‑149 (MQ=255)
                                                                                                                                           gAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAGGGCCGAGGGCAgg          <  3:206764/149‑1 (MQ=255)
                                                                                                                                            aaaCGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAGGGCCGAGGGCAGGt         <  7:52228/149‑1 (MQ=255)
                                                                                                                                                   gTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAGGGCCGAGGGCAGGTACTGCTc  <  6:88311/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCATAGCGGATGGCGATATCTACCGGATCGCGAAAAAGATCGGCGATCTGGTCGGAGAGAAACAAGCGCAGGTTCAACTGTGGATGCTCTCGGCGAAACGCCGTGAGCCAGGGCAGCAGCACGTTGCGCCCCAGATCGGAGGGCGCGGCCACCTGCAACGTTCCGCTGAGCTGGGCGCTTTCCTGACGCAGGTTGTCGCGCCCGTGCCTGAGGGCTTCCAGGGCGATCAGGGCCGAGGGCAGGTACTGCTC  >  NC_002947/3356366‑3356661

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