Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 233,873 G→C 8.9% intergenic (+61/+2) PP_0180 → / ← PP_0181 cytochrome c family protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947233,8730GC 8.9% 189.8 / 7.8 65intergenic (+61/+2)PP_0180/PP_0181cytochrome c family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/31);  new base C (4/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAGGCGTTGCGCAGCCGCGATCAGG  >  NC_002947/233732‑234016
                                                                                                                                             |                                                                                                                                               
cAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCTACTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAggc                                                                                                                                          <  8:124473/149‑1 (MQ=255)
 aaGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGGGTCAGATATTGCGCAAGCCTGATGCTCAggct                                                                                                                                         >  4:275556/1‑149 (MQ=255)
           cTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCCCAAGGACCCCGCAAGGTCTAAGATCTTTCGCAAACCAGATTCTCAGGCTGGCTGTGgcg                                                                                                                               >  6:178800/1‑149 (MQ=255)
           cTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGGCCGCGCAAGGTCTCAGATATTTCGCCCGCCGGAGCCTCAggcttgcgcaagat                                                                                                                               >  6:97672/1‑142 (MQ=255)
                 ggCCATCGTGTTCCTTTACGGGCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCTGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTga                                                                                                                         <  5:167029/149‑1 (MQ=255)
                  gCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTga                                                                                                                         <  6:63397/148‑1 (MQ=255)
                  gCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTg                                                                                                                          <  3:90387/147‑1 (MQ=255)
                    cATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTGCCACAGGGACCGCGCAAGGTCTCAGATATTTCGCAAGCCAGATGCTCAGGCTGGCTGTGGCGCCGTGAGaa                                                                                                                      >  2:93257/1‑149 (MQ=255)
                    cATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGaa                                                                                                                      <  2:272608/149‑1 (MQ=255)
                    cATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGaa                                                                                                                      <  8:236579/149‑1 (MQ=255)
                          gtTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCt                                                                                                                >  6:84085/1‑149 (MQ=255)
                           ttCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTg                                                                                                               <  5:178800/149‑1 (MQ=255)
                             ccTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTGCCACACGGACCCCGCAAGGGGGGAGATATTGCGCAAACCAGATGATCAGGCGGGGGGTGGCGCCGTGAGAATTTTGTGTa                                                                                                             >  2:186097/1‑149 (MQ=255)
                              cTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACATGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTa                                                                                                             <  7:223400/148‑1 (MQ=255)
                              cTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCCAGCCAGATGCTAAGGCTGGCGGTGGCGCCTTGAGAATGTTCTGTAg                                                                                                            >  6:259415/1‑149 (MQ=255)
                              cTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCCGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAg                                                                                                            >  4:242141/1‑149 (MQ=255)
                              cTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTg                                                                                                               <  2:72073/146‑1 (MQ=255)
                              cTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTg                                                                                                               >  1:72073/1‑146 (MQ=255)
                              cTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTa                                                                                                             <  5:84085/148‑1 (MQ=255)
                               tttACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAg                                                                                                            >  1:94068/1‑148 (MQ=255)
                                                      gAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTc                                                                                    <  8:143252/149‑1 (MQ=255)
                                                          aGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGa                                                                                <  1:93257/149‑1 (MQ=255)
                                                          aGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGa                                                                                <  3:78754/149‑1 (MQ=255)
                                                                 tgctgcGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAgg                                                                          <  7:110074/148‑1 (MQ=255)
                                                                      cGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGt                                                                     >  5:79846/1‑148 (MQ=255)
                                                                           aaCTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTa                                                               <  5:259415/149‑1 (MQ=255)
                                                                           aaCTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTa                                                               <  3:242141/149‑1 (MQ=255)
                                                                             ctgactgGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACCGCGTCATGACCCAGGCTGGTCAGGTAgc                                                             <  1:186097/149‑1 (MQ=255)
                                                                                      gCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCAGCCATCCg                                                    >  6:100085/1‑149 (MQ=255)
                                                                                         tGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGcc                                                                                                                             >  4:192384/1‑73 (MQ=255)
                                                                                         tGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGcc                                                                                                                             <  3:192384/73‑1 (MQ=255)
                                                                                                tCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGt                                          <  3:275556/149‑1 (MQ=255)
                                                                                                  ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGt                                          >  6:280238/1‑147 (MQ=255)
                                                                                                  ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTc                                         <  2:94068/148‑1 (MQ=255)
                                                                                                  ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCa                                        <  7:184811/149‑1 (MQ=255)
                                                                                                  ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCa                                        >  1:229146/1‑149 (MQ=255)
                                                                                                   cacaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg                                       >  1:258261/1‑149 (MQ=255)
                                                                                                     caGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGcc                                     <  6:79846/149‑1 (MQ=255)
                                                                                                          aCCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg                                       <  7:27229/142‑1 (MQ=255)
                                                                                                          aCCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg                                       >  8:27229/1‑142 (MQ=255)
                                                                                                           ttgcgcAAGATCTGAGACCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTtct                               >  5:306969/3‑149 (MQ=255)
                                                                                                           ttgcgcAAGATCTGAGACCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTtct                               >  4:293981/3‑149 (MQ=255)
                                                                                                                cAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCtt                                                                                                                           >  4:104115/1‑52 (MQ=255)
                                                                                                                cAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCtt                                                                                                                           <  3:104115/52‑1 (MQ=255)
                                                                                                                   ggTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAgg                       >  1:321414/1‑149 (MQ=255)
                                                                                                                     tctcAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTg                                                                                                      <  5:126817/68‑1 (MQ=255)
                                                                                                                     tctcAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTg                                                                                                      >  6:126817/1‑68 (MQ=255)
                                                                                                                     tctcAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAAc                                                                                          <  6:158576/80‑1 (MQ=255)
                                                                                                                     tctcAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAAc                                                                                          >  5:158576/1‑80 (MQ=255)
                                                                                                                       tcAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGc                                      <  3:202609/130‑1 (MQ=255)
                                                                                                                       tcAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGc                                      >  4:202609/1‑130 (MQ=255)
                                                                                                                       tcAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAg                        >  6:302886/1‑144 (MQ=255)
                                                                                                                       tcAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAg                        <  5:302886/144‑1 (MQ=255)
                                                                                                                       tcAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAGGCGtt                   >  8:69421/1‑149 (MQ=255)
                                                                                                                       tcAGATCTTGCGCAAGCCTGATGCTCAGACTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAGGCGtt                   >  2:260979/1‑149 (MQ=255)
                                                                                                                              ttGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAGGCGTTGCGCAGc            >  7:156133/1‑149 (MQ=255)
                                                                                                                                gcgcAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAGGCGTTGCGCa              <  2:321414/145‑1 (MQ=255)
                                                                                                                                 agccaGCCTGAGCATCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAg                                                                  >  7:61734/5‑92 (MQ=255)
                                                                                                                                 agccaGCCTGAGCATCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAg                                                                  <  8:61734/88‑1 (MQ=255)
                                                                                                                                 acccaGCCTGAGCATCAGGCTCGCTGTGGCGCCATCAGAATGTACACAACCAATTGGACACTATCAACGGCGACATGACCCAGGCAGGTCAGGTAGCCGCCATCCGGCAGAACGCTCAGCCCCTTCTCATGCAGCCGTTGCGCAGCcgc         >  7:238995/5‑149 (MQ=255)
                                                                                                                                  gccaGCCTGAGCATCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAGGCGTTg                  <  2:225827/136‑1 (MQ=255)
                                                                                                                                  gccaGCCTGAGCATCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAGGCGTTg                  >  1:225827/4‑139 (MQ=255)
                                                                                                                                  gcAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAGGCGTTGCGCAGCcgc         >  4:290842/1‑148 (MQ=255)
                                                                                                                                      gCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAGGCGTTGCGCAGCCGCGATca    <  2:258261/149‑1 (MQ=255)
                                                                                                                                        cTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAGGCGTTGCGCAGCCGCGATcagg  <  2:229146/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                               
CAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAGGCGTTGCGCAGCCGCGATCAGG  >  NC_002947/233732‑234016

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑