Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,122,760 A→G 6.3% intergenic (+31/+116) PP_0983 → / ← PP_0984 membrane protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,122,7600AG 6.3% 181.9 / 2.0 63intergenic (+31/+116)PP_0983/PP_0984membrane protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base G (2/2);  total (32/31)
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.66e-01

CTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCACC  >  NC_002947/1122619‑1122897
                                                                                                                                             |                                                                                                                                            
cTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCa                                                                                                                                              <  4:320523/142‑1 (MQ=255)
cTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCa                                                                                                                                              >  3:320523/1‑142 (MQ=255)
cTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCgg                                                                                                                                       <  3:226817/149‑1 (MQ=255)
  gCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCt                                                                                                                                     >  2:29137/1‑149 (MQ=255)
  gCTTGGGCAGTCGAGACAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATATCCTGAGCTGCACAAGAAAGCCGACCTCACGGTCGGCt                                                                                                                                     >  6:18703/1‑149 (MQ=255)
      gggCCGTCGAGCCAGCTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCGGAGCTGCACAAGAAAGCCGACCTCAGGGGCGGCtttc                                                                                                                                  <  1:42139/148‑2 (MQ=255)
       ggCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCtt                                                                                                                                    >  1:141209/1‑145 (MQ=255)
       ggCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCtt                                                                                                                                    <  2:141209/145‑1 (MQ=255)
         ccGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACGAAAAAGCCGACCCTGAGGTCGGCtttcttgt                                                                                                                              >  7:236776/1‑144 (MQ=38)
          cGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACGAAAAAGCCGACCCTGAGGTCGGCTTTCTtgtg                                                                                                                             <  8:236776/149‑1 (MQ=38)
          cGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTg                                                                                                                             <  4:171653/149‑1 (MQ=255)
                 ccAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGtt                                                                                                                       >  5:23475/1‑148 (MQ=255)
                 ccAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGtt                                                                                                                       >  7:182872/1‑148 (MQ=255)
                        ttCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCt                                                                                                                                     >  4:316694/1‑127 (MQ=255)
                        ttCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCt                                                                                                                                     <  3:316694/127‑1 (MQ=255)
                        ttCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTc                                                                                                                                <  4:36830/132‑1 (MQ=255)
                        ttCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTc                                                                                                                                >  3:36830/1‑132 (MQ=255)
                        ttCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGc                                                                                                               >  4:158303/1‑149 (MQ=255)
                          cGGCTTCCCGCGACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGTGTAATAGCCTGAGCTTCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGcc                                                                                                              <  5:18703/148‑1 (MQ=255)
                           ggCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCg                                                                                                            <  6:23475/149‑1 (MQ=255)
                           ggCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCg                                                                                                            <  8:182872/149‑1 (MQ=255)
                             cTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGc                                                                                                          <  1:33681/149‑1 (MQ=255)
                              ttCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCAGAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTg                                                                                                                 >  1:285003/1‑141 (MQ=255)
                              ttCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTg                                                                                                                 <  2:285003/141‑1 (MQ=255)
                                                gtggtgATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGc                                                                                                                                      >  4:18285/1‑102 (MQ=255)
                                                gtggtgATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGc                                                                                                                                      <  3:18285/102‑1 (MQ=255)
                                                   gtgATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCgcgc                                                                                    >  1:344054/1‑149 (MQ=255)
                                                          cAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAgg                                                                          >  8:23152/1‑149 (MQ=255)
                                                           aGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTc                                                                                                    >  7:250180/1‑125 (MQ=255)
                                                           aGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTc                                                                                                    <  8:250180/125‑1 (MQ=255)
                                                             cGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTAc                                                                       >  1:137554/1‑149 (MQ=255)
                                                                         ccTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATc                                                           <  3:158303/149‑1 (MQ=255)
                                                                          cTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCa                                                          <  5:329406/149‑1 (MQ=255)
                                                                                 ggTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAgg                                                   <  7:23152/149‑1 (MQ=255)
                                                                                     tgtgGTTGTTGCGCCGGGCACGGTAATAGCCTGAGCTGCACAAGAAACCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGAAAGCCCGCTCC‑‑‑TACAGGTACCGCGCAGGCATAAGGAAAGGGccc                                               >  8:182306/1‑149 (MQ=255)
                                                                                      gtgGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGccc                                               >  4:235039/1‑148 (MQ=255)
                                                                                       tgGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCg                                              >  5:31690/1‑148 (MQ=255)
                                                                                         gttgttGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAg                                           >  2:14130/1‑149 (MQ=255)
                                                                                            gttgCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGc                                         >  2:282011/1‑148 (MQ=255)
                                                                                               gcgcCGGGCAGTGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGGCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCTCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAAc                                      <  7:182306/148‑1 (MQ=255)
                                                                                                cgcCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACt                                     <  1:29137/148‑1 (MQ=255)
                                                                                                     ggCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACAGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATc                               >  1:267654/1‑149 (MQ=255)
                                                                                                         gggTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCctct                           <  2:344054/149‑1 (MQ=255)
                                                                                                              aTAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCgg                                                                                                >  1:120085/1‑78 (MQ=255)
                                                                                                              aTAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCgg                                                                                                <  2:120085/78‑1 (MQ=255)
                                                                                                                  ccTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGaaa                                                     <  5:46731/114‑1 (MQ=255)
                                                                                                                  ccTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGaaa                                                     >  6:46731/1‑114 (MQ=255)
                                                                                                                   cTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttc                 <  6:31690/149‑1 (MQ=255)
                                                                                                                      aGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttcttc              <  3:235039/149‑1 (MQ=255)
                                                                                                                           cacaAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCACCTCCAATCCTCTGGCTTATTTCTTCTGCtt         >  7:57789/1‑149 (MQ=255)
                                                                                                                           cacaAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTa                     <  4:256342/137‑1 (MQ=255)
                                                                                                                           cacaAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTa                     >  3:256342/1‑137 (MQ=255)
                                                                                                                            acaAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TCCAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTc        >  5:47290/1‑149 (MQ=255)
                                                                                                                             cgaaaaaGCCGACCCTGAGGTCGGCTTTCTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCg       <  3:80460/145‑1 (MQ=255)
                                                                                                                                gAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTtcttct             <  2:292357/140‑1 (MQ=255)
                                                                                                                                gAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTtcttct             >  1:292357/1‑140 (MQ=255)
                                                                                                                                gAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCa    >  7:29231/1‑149 (MQ=255)
                                                                                                                                  aaGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGTTAAGCCCGCGCCCTTTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCgg      <  1:282011/148‑1 (MQ=255)
                                                                                                                                  aaGCCGACCCTGAGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAACGGCCCGAAGGCAACTCCAATCATCTGGCTTATTTCTTATGCTTCGCCAcc  >  4:239944/1‑149 (MQ=255)
                                                                                                                                      cGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttcttc              <  7:314757/133‑1 (MQ=255)
                                                                                                                                      cGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttcttc              >  8:314757/1‑133 (MQ=255)
                                                                                                                                       gACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCgg      >  3:297624/1‑140 (MQ=255)
                                                                                                                                       gACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCgg      <  4:297624/140‑1 (MQ=255)
                                                                                                                                           tCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCctct                           <  6:323276/115‑1 (MQ=255)
                                                                                                                                           tCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCctct                           >  5:323276/1‑115 (MQ=255)
                                                                                                                                             |                                                                                                                                            
CTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCC‑‑‑TACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCACC  >  NC_002947/1122619‑1122897

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: