Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,517,189 A→C 6.3% intergenic (+26/+60) 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,1890AC 6.3% 248.7 / 4.6 79intergenic (+26/+60)infA/clpAtranslation initiation factor IF‑1/ATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base A (36/38);  new base C (2/3);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CGCACATCTCCGGAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGC  >  NC_002947/4517045‑4517336
                                                                                                                                                |                                                                                                                                                    
cgcACATCTCCGGAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTg                                                                                                                                                 <  3:32312/149‑1 (MQ=255)
   acaTCTCCGGAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCg                                                                                                                                              >  4:113313/1‑149 (MQ=255)
         ccGGAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCttt                                                                                                                                        >  4:352958/1‑149 (MQ=255)
          cGGAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCttt                                                                                                                                        <  5:69846/148‑1 (MQ=255)
          cGGAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCttt                                                                                                                                        <  6:163522/148‑1 (MQ=255)
          cGGAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCtt                                                                                                                                         <  7:287319/147‑1 (MQ=255)
          cGGAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCACCTGGCCGGGCtttt                                                                                                                                       <  3:149680/149‑1 (MQ=255)
           ggAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCttttt                                                                                                                                      >  8:310784/1‑149 (MQ=255)
            gAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCtttttt                                                                                                                                     <  4:241061/149‑1 (MQ=255)
             aaaGATTCGCAAGAACTACAACCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCTTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCACCTGGCCGGGCtttttca                                                                                                                                    <  1:274889/149‑3 (MQ=255)
              aaGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTtgt                                                                                                                                   <  7:146486/149‑1 (MQ=255)
                  tGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCt                                                                                                                               <  3:352958/149‑1 (MQ=255)
                         gAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGaaa                                                                                                                        >  4:265067/1‑149 (MQ=255)
                         gAACTACAACCGCATTCTCACTGGCGACAAGGTCCGCGTCGCACTGACGCCGTACGCCCTGAGCAAGGGCCCCATCACCTACCCTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCCCCCAGGTGCCCGGCTTTTTTGTGCCTCAAGaaa                                                                                                                        >  5:286872/1‑149 (MQ=255)
                           aCTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAAtt                                                                                                                      <  4:310671/149‑1 (MQ=255)
                                  ccGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCgg                                                                                                                                             >  1:52098/1‑119 (MQ=255)
                                  ccGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCgg                                                                                                                                             <  2:52098/119‑1 (MQ=255)
                                   cGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGc                                                                                                              >  1:360026/1‑149 (MQ=255)
                                   cGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGc                                                                                                              >  4:242646/1‑149 (MQ=255)
                                    gCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCt                                                                                                             >  3:265981/1‑149 (MQ=255)
                                            aCTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCg                                                                                                                                              >  5:94446/1‑108 (MQ=255)
                                            aCTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCg                                                                                                                                              <  6:94446/108‑1 (MQ=255)
                                            aCTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCa                                                                                                     <  2:305112/149‑1 (MQ=255)
                                              tGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGc                                                                                                   >  1:229490/1‑149 (MQ=255)
                                                gCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGc                                                                                                   <  2:209058/147‑1 (MQ=255)
                                                   aCAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTg                                                                                              <  7:155350/149‑1 (MQ=255)
                                                   aCAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTg                                                                                              <  4:265981/149‑1 (MQ=255)
                                                        gTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGc                                                                                                   <  8:25483/139‑1 (MQ=255)
                                                        gTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGc                                                                                                   >  7:25483/1‑139 (MQ=255)
                                                          ccGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAGAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAgg                                                                                                                 >  7:242940/1‑123 (MQ=255)
                                                          ccGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAACAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAgg                                                                                                                 <  8:242940/123‑1 (MQ=255)
                                                          ccGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCa                                                                                       >  8:151489/1‑149 (MQ=255)
                                                             cgTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCAcgc                                                                                    >  8:47620/1‑149 (MQ=255)
                                                             cgTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCAcgc                                                                                    <  5:188229/149‑1 (MQ=255)
                                                              gTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAgggg                                                                                                               >  8:81068/1‑121 (MQ=255)
                                                              gTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAgggg                                                                                                               <  7:81068/121‑1 (MQ=255)
                                                                 gAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGAcc                                                                                <  3:113313/149‑1 (MQ=255)
                                                                 gAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGAcc                                                                                >  8:195924/1‑149 (MQ=255)
                                                                  aaCTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCACCTGGCCGGGCtttttcat                                                                                                                                   <  4:187114/97‑4 (MQ=37)
                                                                  aaCTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCACCTGGCCGGGCtttttcat                                                                                                                                   >  3:187114/1‑94 (MQ=37)
                                                                             tACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTCGCCGGGCTTTTTTGTGCCTCAAGAAAtt                                                                                                                      <  1:203539/100‑1 (MQ=255)
                                                                             tACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAAtt                                                                                                                      >  2:203539/1‑99 (MQ=255)
                                                                              aCGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGttt                                                                   >  5:129602/1‑149 (MQ=255)
                                                                               cGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTc                                                                  <  2:229490/149‑1 (MQ=255)
                                                                                gACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTc                                                                  >  1:321799/1‑148 (MQ=255)
                                                                                  ccTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGc                                                                                                   <  7:249625/113‑1 (MQ=255)
                                                                                  ccTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGc                                                                                                   >  8:249625/1‑113 (MQ=255)
                                                                                  ccTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCg                                                                 >  5:33727/1‑147 (MQ=255)
                                                                                   cTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAg                                                              <  7:195924/149‑1 (MQ=255)
                                                                                   cTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAg                                                              <  6:129602/149‑1 (MQ=255)
                                                                                   cTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAg                                                              >  4:192401/1‑149 (MQ=255)
                                                                                     gAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTc                                                            >  3:137365/1‑149 (MQ=255)
                                                                                     gAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTc                                                            >  2:60219/1‑149 (MQ=255)
                                                                                      aGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGt                                                             <  3:242646/147‑1 (MQ=255)
                                                                                      aGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCg                                                           <  2:360026/149‑1 (MQ=255)
                                                                                       gCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGa                                                          <  4:137365/149‑1 (MQ=255)
                                                                                           gggCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAg                                                              <  5:127109/141‑1 (MQ=255)
                                                                                           gggCCGCATCACCTACCGTGAGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAg                                                              >  6:127109/1‑141 (MQ=255)
                                                                                             gCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCa                                                    <  3:265067/149‑1 (MQ=255)
                                                                                                gCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGtt                                                 <  7:151489/149‑1 (MQ=255)
                                                                                                 cATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTc                                                >  4:60983/1‑149 (MQ=255)
                                                                                                     aCCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTc                                                                              <  5:50402/115‑1 (MQ=255)
                                                                                                     aCCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTc                                                                              >  6:50402/1‑115 (MQ=255)
                                                                                                     aCCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGcc                                             >  6:331904/1‑148 (MQ=255)
                                                                                                           cGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCgcgc                                      <  2:321799/149‑1 (MQ=255)
                                                                                                             tGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCAGATGACGCGACCTCAGCCGTGGTTTCGAAGTCGTACACCAGTTCGCCATCGCGCAg                                    <  6:286872/149‑1 (MQ=255)
                                                                                                               cgcgcTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAgg                                   <  7:47620/148‑1 (MQ=255)
                                                                                                                    tAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACa                             >  2:351368/1‑149 (MQ=255)
                                                                                                                      aGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATg                           >  5:112197/1‑149 (MQ=255)
                                                                                                                       gCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGc                          <  7:310784/149‑1 (MQ=255)
                                                                                                                        cTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGca                         <  1:351368/149‑1 (MQ=255)
                                                                                                                          cAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGcacc                       >  1:337940/1‑149 (MQ=255)
                                                                                                                               cATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGcc                  >  1:235694/1‑149 (MQ=255)
                                                                                                                                 tGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGccac                >  7:285923/1‑149 (MQ=255)
                                                                                                                                  aaaaaaGCCCGGCACCTGGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGccac                >  8:250626/2‑148 (MQ=255)
                                                                                                                                           cGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGc      <  8:285923/149‑1 (MQ=255)
                                                                                                                                              ccAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCgcgc                                      >  4:196285/1‑114 (MQ=255)
                                                                                                                                              ccAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCgcgc                                      <  3:196285/114‑1 (MQ=255)
                                                                                                                                               cAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGc  >  3:168160/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CGCACATCTCCGGAAAGATGCGCAAGAACTACATCCGCATTCTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGC  >  NC_002947/4517045‑4517336

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: