Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 797,509 G→A 7.1% 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 7.1% 206.9 / 3.1 70intergenic (+33/‑44)PP_0682/PP_0683inner membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/31);  new base A (2/3);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

ATGGATGCGCAGAACATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCATCAGCTCAA  >  NC_002947/797363‑797656
                                                                                                                                                  |                                                                                                                                                    
aTGGATGCGCAGAACATCTCCACCGATTACCTGATGTGATACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGgggc                                                                                                                                                   >  3:169626/1‑149 (MQ=255)
       cgcAGAACATCTCCCCCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCAt                                                                                                                                             <  3:210634/148‑1 (MQ=255)
           gAACATCTCCCCCGATAACCTGATGTGGTACGTGATCATCCACCTGACCTTCGTGGGCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCg                                                                                                                                        <  3:185204/149‑1 (MQ=255)
            aaCATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGc                                                                                                                                       >  3:266052/1‑149 (MQ=255)
              cATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCATGCCCGccc                                                                                                                                     <  5:306472/149‑3 (MQ=255)
                  tcCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCATGCCCGCCCCca                                                                                                                                  >  2:213248/1‑148 (MQ=255)
                   ccACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTccc                                                                                                                                   <  4:276586/146‑1 (MQ=255)
                   ccACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTccc                                                                                                                                   >  3:276586/1‑146 (MQ=255)
                   ccACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCca                                                                                                                                  >  6:276411/1‑147 (MQ=255)
                    cACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGTGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCcaca                                                                                                                                <  2:138902/148‑1 (MQ=255)
                    cACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGTGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCca                                                                                                                                  >  1:138902/1‑146 (MQ=255)
                    cACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCca                                                                                                                                  >  1:196819/1‑146 (MQ=255)
                    cACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCca                                                                                                                                  >  3:26362/1‑146 (MQ=255)
                      ccGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCATGCCCGCCCCCACAgga                                                                                                                             <  1:213248/149‑3 (MQ=255)
                          ttACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCg                                                                                                                         <  4:266052/149‑1 (MQ=255)
                             ccTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAg                                                                                                                       <  2:259638/148‑1 (MQ=255)
                             ccTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCAt                                                                                                                                             >  3:226436/1‑126 (MQ=255)
                             ccTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCAt                                                                                                                                             <  4:226436/126‑1 (MQ=255)
                              ctcaTGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTCGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGCGGCGGGCATGcc                                                                                                                                          >  3:39402/4‑128 (MQ=255)
                              cTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAgg                                                                                                                      >  6:289104/1‑148 (MQ=255)
                               tGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTcc                                                                                                                                    <  3:294437/133‑1 (MQ=255)
                               tGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTcc                                                                                                                                    >  4:294437/1‑133 (MQ=255)
                                     ggTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGCGCGGGCATGCCCGCTCCCACAGATATCGa                                                                                                                        <  1:24453/140‑1 (MQ=255)
                                     ggTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGa                                                                                                                        >  2:24453/1‑139 (MQ=255)
                                      gTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCATCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTcc                                                                                                                                    >  6:322680/1‑126 (MQ=255)
                                      gTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCATCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTcc                                                                                                                                    <  5:322680/126‑1 (MQ=255)
                                         cGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGa                                                                                                          <  5:276411/149‑1 (MQ=255)
                                         cGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGa                                                                                                          <  4:169626/149‑1 (MQ=255)
                                             atcatcCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTg                                                                                                       <  5:289104/148‑1 (MQ=255)
                                                  ccACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACAAAGCACTGAGTCCCACCAGGCCTCACTGCAACTCCTGTGGTGGCGGGCATTCCCGCTACCACAGATATCGAGGCTCAACTGACGAGTGGCgt                                                                                                   >  5:292533/1‑147 (MQ=255)
                                                   cacaTGAGATTCGTGGTCTCGGCGTTCTGCATTGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGAAACTCCTGTGGGGGCGGGAATGCCCGATCCCACAGATAACGATGCTCAACTGACGAATGGCgtcag                                                                                                >  6:292445/1‑149 (MQ=255)
                                                    acaTGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCgtcagt                                                                                               <  2:247590/149‑1 (MQ=255)
                                                      aTGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCgtcagtca                                                                                             >  8:243874/1‑149 (MQ=255)
                                                       tGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCAt                                                                                            <  1:73560/149‑1 (MQ=255)
                                                                   tctcGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGa                                                                                                          >  1:296298/1‑123 (MQ=255)
                                                                   tctcGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGa                                                                                                          <  2:296298/123‑1 (MQ=255)
                                                                    ctcGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATtc                                                                               >  6:134240/1‑149 (MQ=255)
                                                                     tcGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTCACGAGTGGCGTCAGTCATGCCTTGCGCATtct                                                                              >  8:102673/1‑149 (MQ=255)
                                                                       ggCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATtctc                                                                             <  7:102673/148‑1 (MQ=255)
                                                                       ggCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATtctc                                                                             <  4:274975/148‑1 (MQ=255)
                                                                       ggCGGTCTGCAGGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCACGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATtctc                                                                             <  8:297305/148‑1 (MQ=255)
                                                                          gTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCtt                                                                         >  8:87628/1‑149 (MQ=255)
                                                                                cATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCt                                                                      <  5:134240/146‑1 (MQ=255)
                                                                                    ggCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACaa                                                                >  6:45946/1‑148 (MQ=255)
                                                                                         ccTGGACAAGCTGACCAAGCACTGAGTCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAAt                                                           >  3:124304/1‑147 (MQ=255)
                                                                                           tcgACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTACCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAg                                                        >  8:170927/3‑149 (MQ=255)
                                                                                            ggACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGc                                                       >  5:171262/1‑149 (MQ=255)
                                                                                               cAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGACctc                                                     >  2:305359/1‑148 (MQ=255)
                                                                                                  gCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGa                                                         >  1:157401/1‑141 (MQ=255)
                                                                                                  gCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGa                                                         <  2:157401/141‑1 (MQ=255)
                                                                                                  gCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGct                                                      <  1:34569/144‑1 (MQ=255)
                                                                                                  gCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGct                                                      >  2:34569/1‑144 (MQ=255)
                                                                                                        cAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTgg                                           <  4:26362/149‑1 (MQ=255)
                                                                                                        cAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTgg                                           >  1:95077/1‑149 (MQ=255)
                                                                                                          aGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGtt                                         >  7:110708/1‑149 (MQ=255)
                                                                                                             aCTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTccc                                      <  4:246112/149‑1 (MQ=255)
                                                                                                             aCTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTcc                                       >  3:246112/1‑148 (MQ=255)
                                                                                                               tGAGTCCCCCAATGCCTCACTGCAACTCCTGTGGTGGGGGTCATGCCCGCTCCCACATATATCTAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATg                                                          <  3:118965/127‑1 (MQ=255)
                                                                                                                 tttCCGCCAAGGCCTCACTGCATCTCCTGTGGGGGCGGGCATTCCCGCTCCCACAGATATCGATGCTCAACTTACGAGTGGGGTCAGTCATGCCTTGCTCATTCTCCCTTCCTTTAGAAAAAATGAGCTCTCATGCGTGGTTCCCAta                                   <  5:292445/146‑2 (MQ=255)
                                                                                                                 aGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATc                                   <  6:171262/148‑1 (MQ=255)
                                                                                                                 aGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATc                                   <  6:175971/148‑1 (MQ=255)
                                                                                                                  gTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATcgc                                 >  2:212645/1‑149 (MQ=255)
                                                                                                                        ccAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTg                           <  5:45946/149‑1 (MQ=255)
                                                                                                                            gCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGctc                        >  7:163346/1‑148 (MQ=255)
                                                                                                                             ccTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCtctc                      <  4:124304/149‑1 (MQ=255)
                                                                                                                                 aCTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATg                   >  6:121786/1‑148 (MQ=255)
                                                                                                                                        cTCCTGTGGGGGCGGGCATGCCCGCTCCCACATATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCa           <  2:196819/149‑1 (MQ=255)
                                                                                                                                        cTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCa           >  6:293929/1‑149 (MQ=255)
                                                                                                                                         tCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCAt          >  1:311100/1‑149 (MQ=255)
                                                                                                                                           cTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCAt          <  2:95077/147‑1 (MQ=255)
                                                                                                                                            tgtgGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCATCAg       >  2:307955/1‑149 (MQ=255)
                                                                                                                                               gggggCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATg                   >  2:98024/1‑134 (MQ=255)
                                                                                                                                               gggggCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATg                   <  1:98024/134‑1 (MQ=255)
                                                                                                                                                 gggcgggcATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCATCAGCTCaa  <  7:87628/149‑1 (MQ=255)
                                                                                                                                                  ggcgggcATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCATCAGCTCaa  >  5:270527/1‑148 (MQ=255)
                                                                                                                                                  ggcgggcATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCATCAGCTCaa  >  4:27508/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATGGATGCGCAGAACATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCATCAGCTCAA  >  NC_002947/797363‑797656

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: