Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,517,191 G→T 7.9% intergenic (+28/+58) infA → / ← clpA translation initiation factor IF‑1/ATP‑dependent serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,517,1910GT 7.9% 285.2 / 11.4 88intergenic (+28/+58)infA/clpAtranslation initiation factor IF‑1/ATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base G (38/43);  new base T (3/4);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-01

GCGCACATCTCCGGAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCT  >  NC_002947/4517044‑4517337
                                                                                                                                                   |                                                                                                                                                  
gcgcACATCTCCGGAAAGATGCGCAAGAACTACATCCGCATTCTCGCTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGt                                                                                                                                                   <  6:84585/149‑1 (MQ=255)
    acaTCTCCGGAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCg                                                                                                                                               >  8:13329/1‑149 (MQ=255)
     caTCTCCGGAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCgg                                                                                                                                              <  4:170631/149‑1 (MQ=255)
     caTCTCCGGAAAGATGCGCAAGAACTACAACCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCgg                                                                                                                                              <  3:383440/149‑1 (MQ=255)
      aTCTCCGGAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCggg                                                                                                                                             >  5:235734/1‑149 (MQ=255)
            ggAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCttttt                                                                                                                                       <  3:349124/149‑1 (MQ=255)
            ggAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCttttt                                                                                                                                       >  2:168119/1‑149 (MQ=255)
             gAAAGATGTGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCACCTGGCCGGGCtttttc                                                                                                                                      <  1:142823/149‑2 (MQ=255)
               aaGATGCTCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCTTCGAACTGACTCCGTACGACCTTAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTtgt                                                                                                                                    <  3:45861/149‑1 (MQ=255)
               aaGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTtgt                                                                                                                                    <  1:209369/149‑1 (MQ=255)
               aaGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTtgt                                                                                                                                    <  2:237961/149‑1 (MQ=255)
                aGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTtgtg                                                                                                                                   <  5:360400/149‑1 (MQ=255)
                aGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTtgtg                                                                                                                                   <  3:66757/149‑1 (MQ=255)
                     cgcAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCa                                                                                                                              <  3:315238/149‑1 (MQ=255)
                        aaGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGCGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAg                                                                                                                            >  6:10630/1‑148 (MQ=255)
                             cTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTg                                                                                                                      <  7:13329/149‑1 (MQ=255)
                              tACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTg                                                                                                                      >  4:298408/1‑148 (MQ=255)
                                 aTCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGc                                                                                                                                  >  2:357072/1‑133 (MQ=255)
                                 aTCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGc                                                                                                                                  <  1:357072/133‑1 (MQ=255)
                                 aTCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAgg                                                                                                                  <  5:29877/149‑1 (MQ=255)
                                 aTCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAgg                                                                                                                  <  7:174175/149‑1 (MQ=255)
                                   ccGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGCCGCAGTAAGACCTGAGCAAGGGCCTCATCACCTACCGAGCTCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCATCAAGAAATTGCAgggg                                                                                                                >  6:252840/1‑149 (MQ=255)
                                     gCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGt                                                                                                                                                   <  5:225763/112‑1 (MQ=255)
                                     gCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGt                                                                                                                                                   >  6:225763/1‑112 (MQ=255)
                                     gCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCt                                                                                                              >  5:305216/1‑149 (MQ=255)
                                      cATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTg                                                                                                             <  5:166790/149‑1 (MQ=255)
                                       aTTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGc                                                                                                            <  3:153227/149‑1 (MQ=255)
                                       aTTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGc                                                                                                            <  3:147306/149‑1 (MQ=255)
                                          ctcACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCaa                                                                                                          <  6:235734/148‑1 (MQ=255)
                                             aCTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGGCTCAAGAAATTGCAGGGGCTGCAAGGCa                                                                                                      <  4:352517/149‑1 (MQ=255)
                                                 gCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCGACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGccc                                                                                                  >  1:246947/1‑149 (MQ=255)
                                                     cAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGc                                                                                              <  2:246947/149‑1 (MQ=255)
                                                         gTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACAt                                                                                          <  7:202848/149‑1 (MQ=255)
                                                          tCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAgg                                                                                                                  >  4:233895/1‑124 (MQ=255)
                                                          tCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAgg                                                                                                                  <  3:233895/124‑1 (MQ=255)
                                                          tCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCa                                                                                                      >  1:163912/1‑136 (MQ=255)
                                                                 cGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGaaa                                                                                                                         >  3:57212/1‑110 (MQ=255)
                                                                 cGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGaaa                                                                                                                         <  4:57212/110‑1 (MQ=255)
                                                                 cGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAgggg                                                                                                                <  5:214435/119‑1 (MQ=255)
                                                                 cGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAgggg                                                                                                                >  6:214435/1‑119 (MQ=255)
                                                                   aaCTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCAGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCCCGCGGCCt                                                                                >  2:413352/1‑149 (MQ=255)
                                                                        aCGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGcc                                                                           >  7:94906/1‑149 (MQ=255)
                                                                               aCGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGtt                                                                     >  8:149036/1‑148 (MQ=255)
                                                                                cGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTc                                                                   <  8:94906/149‑1 (MQ=255)
                                                                                 gACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCaa                                                                                                                             >  7:384592/1‑90 (MQ=255)
                                                                                 gACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCaa                                                                                                                             <  8:384592/90‑1 (MQ=255)
                                                                                  aCCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCg                                                                  >  7:278314/1‑148 (MQ=255)
                                                                                  aCCTGAGCAAGGGCCGCATCACCTACCGTGCGCACTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCg                                                                  >  4:315669/1‑148 (MQ=255)
                                                                                   ccTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCg                                                                  >  8:216825/1‑147 (MQ=255)
                                                                                      gAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTc                                                             >  6:44683/1‑149 (MQ=255)
                                                                                       aGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTc                                                             >  4:94560/1‑148 (MQ=255)
                                                                                          aaGGGCCGCATCTCATACCGTGCGCGCTTAGCTCAAGCCATGAAAAAGCCCGGCCAGTTGCCGGGCTTTTTTGTTCCTCAAGAAATTGCAGGGGCTTGAAGGCATCCCCTTCACATCACTCGACCTCAGCCGTGGTTTCGAAGTCGAac                                                         <  5:252840/149‑1 (MQ=255)
                                                                                               ccGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTgg                                                                       <  4:113971/130‑1 (MQ=255)
                                                                                               ccGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTgg                                                                       >  3:113971/1‑130 (MQ=255)
                                                                                                  cATCACCTACCGCGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTc                                                 <  5:10630/149‑1 (MQ=255)
                                                                                                      aCCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCa                                             <  3:298408/149‑1 (MQ=255)
                                                                                                         tACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTc                                                             >  8:109142/1‑130 (MQ=255)
                                                                                                         tACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTc                                                             <  7:109142/130‑1 (MQ=255)
                                                                                                          aCCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGaa                                                                >  4:99729/1‑126 (MQ=255)
                                                                                                          aCCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGaa                                                                <  3:99729/126‑1 (MQ=255)
                                                                                                                   gcTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGa                                <  6:305216/149‑1 (MQ=255)
                                                                                                                     tAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGaa                                                                <  7:156921/115‑1 (MQ=255)
                                                                                                                     tAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGaa                                                                >  8:156921/1‑115 (MQ=255)
                                                                                                                       aGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACAt                             <  3:94560/148‑1 (MQ=255)
                                                                                                                       aGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACAt                             <  8:278314/148‑1 (MQ=255)
                                                                                                                          tCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAgg                                                                                                                  >  1:140558/1‑60 (MQ=255)
                                                                                                                          tCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAgg                                                                                                                  <  2:140558/60‑1 (MQ=255)
                                                                                                                          tCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGTCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGcac                         >  6:289953/1‑149 (MQ=255)
                                                                                                                          tCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATg                            <  8:255342/146‑1 (MQ=255)
                                                                                                                          tCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATg                            >  7:255342/1‑146 (MQ=255)
                                                                                                                               ccATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGc                    >  4:308628/1‑149 (MQ=255)
                                                                                                                                cacaaaaaaGCCCGGCACCTGGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCg                                                                  <  5:167317/98‑1 (MQ=37)
                                                                                                                                cacaaaaaaGCCCGGCACCTGGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGAGACCTCAGCCGTGGTTTCg                                                                  >  6:167317/5‑102 (MQ=37)
                                                                                                                                cATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGCCATGCACCACGcc                   >  3:230166/1‑149 (MQ=255)
                                                                                                                                 aTGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATGACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGcca                  <  1:127697/149‑1 (MQ=255)
                                                                                                                                 aTGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGcca                  <  5:44683/149‑1 (MQ=255)
                                                                                                                                  tGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGc                           >  4:112568/1‑139 (MQ=255)
                                                                                                                                  tGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGc                           <  3:112568/139‑1 (MQ=255)
                                                                                                                                   gAAAAAGCCCGGCCAGGAGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGACCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCATGTCGACATGCACCACGccac                 >  1:404424/1‑148 (MQ=255)
                                                                                                                                   aaaaaaGCCCGGCACCTGGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGccac                 >  6:34046/2‑148 (MQ=255)
                                                                                                                                   aaaaaaGCCCGGCACCTGGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGcca                  <  7:96351/146‑1 (MQ=255)
                                                                                                                                   aaaaaaGCCCGGCACCTGGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGcca                  >  8:96351/2‑147 (MQ=255)
                                                                                                                                    aaaaaGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGTAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGccacc                <  8:254542/148‑1 (MQ=255)
                                                                                                                                     aaaaGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCAt              >  4:341736/1‑149 (MQ=255)
                                                                                                                                     aaaaGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCAt              >  6:157768/1‑149 (MQ=255)
                                                                                                                                         gCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGACCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTc          >  1:202010/1‑149 (MQ=255)
                                                                                                                                                cAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGc   <  1:168119/149‑1 (MQ=255)
                                                                                                                                                 aGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCt  >  1:289551/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCGCACATCTCCGGAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCT  >  NC_002947/4517044‑4517337

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: