Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 797,509 G→A 6.4% intergenic (+33/‑44) PP_0682 → / → PP_0683 inner membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947797,5090GA 6.4% 188.9 / 2.0 62intergenic (+33/‑44)PP_0682/PP_0683inner membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/26);  new base A (2/2);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

CATGGATGCGCAGAACATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCATCAGCTCAAT  >  NC_002947/797362‑797657
                                                                                                                                                   |                                                                                                                                                    
cATGGATGCGCAGAACATCTCCACCGATTACCTGATGTGGTACGTGATCTTCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCCCACTGCAACTCCTGTggggg                                                                                                                                                     >  4:332672/1‑149 (MQ=255)
    gATGCGCAGAACATCTCCACCGATTACCTGATGTGGTACGTGATCTTCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCCCACTGCAACTCCTGTGGgggcgg                                                                                                                                                  <  3:332672/148‑1 (MQ=255)
    gATGCGCAGAACATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGgggcggg                                                                                                                                                 <  4:81722/149‑1 (MQ=255)
         gcAGAACATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGc                                                                                                                                            <  7:129515/149‑1 (MQ=255)
          cAGAACATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGcc                                                                                                                                           >  4:195839/1‑149 (MQ=255)
           aGAACATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGccc                                                                                                                                          <  2:337109/149‑1 (MQ=255)
             aaCATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGc                                                                                                                                        >  2:125830/1‑149 (MQ=255)
                     cACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCca                                                                                                                                   >  6:354743/1‑146 (MQ=255)
                     cACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCca                                                                                                                                   >  4:345523/1‑146 (MQ=255)
                             aacTGATGTGGTACGTGATCATCCACATGACCTTCGTGGGCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGa                                                                                                                         <  3:50029/146‑1 (MQ=255)
                             aCCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGa                                                                                                                         <  4:98389/148‑1 (MQ=255)
                              ccTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAgg                                                                                                                       <  5:130390/149‑1 (MQ=255)
                              ccTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAg                                                                                                                        <  8:87182/148‑1 (MQ=255)
                               cTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGa                                                                                                                         <  7:119408/146‑1 (MQ=255)
                               cTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGa                                                                                                                         >  8:119408/1‑146 (MQ=255)
                               cTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGc                                                                                                                      >  6:3377/1‑149 (MQ=255)
                               cTGATGTGGTACGTGATCATCCACATGACCTTCGTGGGCTCGGCGTTCTGCATGGGCTAGCTGGACAAGCTGACCAAGCTCTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAggg                                                                                                                      <  6:11683/149‑2 (MQ=255)
                                tGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCATGCCCGCCCCcaca                                                                                                                                 >  7:348962/1‑137 (MQ=255)
                                tGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCATGCCCGCCCCcaca                                                                                                                                 <  8:348962/137‑1 (MQ=255)
                                   tgtgGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCaa                                                                                                                  <  3:345523/149‑1 (MQ=255)
                                        tACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGAc                                                                                                             <  3:195839/149‑1 (MQ=255)
                                        tACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGAc                                                                                                             >  3:350445/1‑149 (MQ=255)
                                           gTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAg                                                                                                          >  1:332011/1‑149 (MQ=255)
                                              atcatcCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTgg                                                                                                       >  6:113107/1‑149 (MQ=255)
                                                 atcCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAACCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCg                                                                                                     <  2:318835/148‑1 (MQ=255)
                                                  tcCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCgt                                                                                                    >  6:160061/1‑148 (MQ=255)
                                                  tcCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCAt                                                                                                                                              <  8:111150/106‑1 (MQ=255)
                                                  tcCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCAt                                                                                                                                              >  7:111150/1‑106 (MQ=255)
                                                   ccACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCgt                                                                                                    >  5:312983/1‑147 (MQ=255)
                                                     acaTGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCgtcagt                                                                                                >  8:53203/1‑149 (MQ=255)
                                                         gACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGa                                                                                                                         >  8:188744/1‑120 (MQ=255)
                                                         gACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGa                                                                                                                         <  7:188744/120‑1 (MQ=255)
                                                                gTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTgc                                                                                      <  6:312983/148‑1 (MQ=255)
                                                                gTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTgc                                                                                      <  4:350445/148‑1 (MQ=255)
                                                                 tGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCcaca                                                                                                                                 <  7:280344/104‑1 (MQ=255)
                                                                 tGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCcaca                                                                                                                                 >  8:280344/1‑104 (MQ=255)
                                                                    tctcGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCgtca                                                                                                  >  4:123703/1‑132 (MQ=255)
                                                                    tctcGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCgtca                                                                                                  <  3:123703/132‑1 (MQ=255)
                                                                    tctcGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCAtt                                                                                 <  5:160061/149‑1 (MQ=255)
                                                                      tcGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGgggcggg                                                                                                                                                 >  7:294916/1‑83 (MQ=255)
                                                                      tcGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGgggcggg                                                                                                                                                 <  8:294916/83‑1 (MQ=255)
                                                                      tcGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATAGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATtct                                                                               >  7:288572/1‑149 (MQ=255)
                                                                        ggCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTcc                                                                             <  1:125830/149‑1 (MQ=255)
                                                                        ggCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTcc                                                                             <  5:354743/149‑1 (MQ=255)
                                                                           gTTATGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCATAACTGCCACTACTGAGGGGGCGGGCATGCCCGCTCCCACAGATAACGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTATCCCta                                                                          >  5:258741/1‑148 (MQ=255)
                                                                                 cATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGt                                                                     <  5:3377/148‑1 (MQ=255)
                                                                                 cATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTa                                                                    <  2:332011/149‑1 (MQ=255)
                                                                                          ccTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAAt                                                            >  6:351005/1‑148 (MQ=255)
                                                                                           cTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGa                                                          <  5:351005/149‑1 (MQ=255)
                                                                                             ggACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCAtt                                                                                 >  4:100143/1‑124 (MQ=255)
                                                                                              gACATGGTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGc                                                        <  7:130036/148‑1 (MQ=255)
                                                                                                   gCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATTCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCAt                                                  >  6:49077/1‑149 (MQ=255)
                                                                                                   gCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCAt                                                  >  1:239364/1‑149 (MQ=255)
                                                                                                       aCCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGt                                              >  2:216316/1‑149 (MQ=255)
                                                                                                                  aGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATcg                                   <  5:113107/149‑1 (MQ=255)
                                                                                                                  aGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATcg                                   >  8:226506/1‑149 (MQ=255)
                                                                                                                      cccccAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGa                                <  2:239364/148‑1 (MQ=255)
                                                                                                                      cccccAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAg                               >  8:23050/1‑149 (MQ=255)
                                                                                                                       ccccAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAgg                              >  1:243493/1‑149 (MQ=255)
                                                                                                                        cccAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGt                             >  1:107832/1‑149 (MQ=255)
                                                                                                                               cTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATcgc                                  >  6:35851/1‑137 (MQ=255)
                                                                                                                               cTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATcgc                                  <  5:35851/137‑1 (MQ=255)
                                                                                                                               cTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCa                      >  2:238572/1‑149 (MQ=255)
                                                                                                                               cTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCa                      >  7:68422/1‑149 (MQ=255)
                                                                                                                                  aCTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATg                    >  3:328231/1‑148 (MQ=255)
                                                                                                                                      cAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCt               >  7:148780/1‑149 (MQ=255)
                                                                                                                                            cTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCAt           >  3:193848/1‑147 (MQ=255)
                                                                                                                                                   ggcgggcTTGGCCGCTCCCACAGATATCGCGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCATCAGCTCAAt  <  8:68422/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CATGGATGCGCAGAACATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCATCAGCTCAAT  >  NC_002947/797362‑797657

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: