Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 233,874 C→A 8.8% intergenic (+62/+1) 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,8740CA 8.8% 202.1 / 6.9 65intergenic (+62/+1)PP_0180/PP_0181cytochrome c family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base A (4/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAGGCGTTGCGCAGCCGCGATCAGGGCAGG  >  NC_002947/233732‑234021
                                                                                                                                              |                                                                                                                                                   
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                                                                                                                           <  8:236579/149‑1 (MQ=255)
                    cATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGaa                                                                                                                           <  2:272608/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                                             >  1:229146/1‑149 (MQ=255)
                                                                                                  ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCa                                             <  7:184811/149‑1 (MQ=255)
                                                                                                   cacaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg                                            >  1:258261/1‑149 (MQ=255)
                                                                                                     caGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGcc                                          <  6:79846/149‑1 (MQ=255)
                                                                                                          aCCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg                                            >  8:27229/1‑142 (MQ=255)
                                                                                                          aCCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg                                            <  7:27229/142‑1 (MQ=255)
                                                                                                           ttgcgcAAGATCTGAGACCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTtct                                    >  4:293981/3‑149 (MQ=255)
                                                                                                           ttgcgcAAGATCTGAGACCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTtct                                    >  5:306969/3‑149 (MQ=255)
                                                                                                                cAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCtt                                                                                                                                <  3:104115/52‑1 (MQ=255)
                                                                                                                cAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCtt                                                                                                                                >  4:104115/1‑52 (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                                           >  4:202609/1‑130 (MQ=255)
                                                                                                                       tcAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGc                                           <  3:202609/130‑1 (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)
                                                                                                                                              cTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAGGCGTTGCGCAGCCGCGATcagggcagg  <  7:69421/148‑1 (MQ=255)
                                                                                                                                              cTCAGACTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAGGCGTTGCGCAGCCGCGATcagggcagg  <  1:260979/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCCCCTTCTCATGCAGGCGTTGCGCAGCCGCGATCAGGGCAGG  >  NC_002947/233732‑234021

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