Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 444,580 G→C 5.9% intergenic (+59/‑115) bioD → / → PP_0367 dethiobiotin synthetase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947444,5800GC 5.9% 540.7 / 12.4 169intergenic (+59/‑115)bioD/PP_0367dethiobiotin synthetase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (66/92);  new base C (7/4);  total (73/96)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 170 total aligned reads displayed.

CGGCGCCTTGCCTGGGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTC  >  NC_002947/444433‑444728
                                                                                                                                                   |                                                                                                                                                    
cGGCGCCTTGCCTGGGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGtt                                                                                                                                                    <  1:4155001/150‑1 (MQ=255)
 ggCGCCTTGCCTGGGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTtg                                                                                                                                                   >  2:2855630/1‑150 (MQ=255)
  gcgcCTTGCCTGGGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTtgt                                                                                                                                                  <  2:1498657/150‑1 (MQ=255)
   cgcCTTGCCTGGGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTtgtg                                                                                                                                                 <  2:1766023/150‑1 (MQ=255)
    gcCTTGCCTGGGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTgg                                                                                                                                                <  2:4956669/150‑1 (MQ=255)
    gcCTTGCCTGGGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTgg                                                                                                                                                <  1:966720/150‑1 (MQ=255)
      cTTGCCTGGGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGgc                                                                                                                                              >  2:1229740/1‑150 (MQ=255)
          ccTGGGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCgg                                                                                                                                          <  2:3462689/150‑1 (MQ=255)
          ccTGGGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCgg                                                                                                                                          <  1:650457/150‑1 (MQ=255)
          ccTGGGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCgg                                                                                                                                          <  1:4694641/150‑1 (MQ=255)
          ccTGGGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCgg                                                                                                                                          <  1:650473/150‑1 (MQ=255)
          ccTGGGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCgg                                                                                                                                          <  2:3462783/150‑1 (MQ=255)
           cTGGGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCggg                                                                                                                                         <  1:2808808/150‑1 (MQ=255)
           cTGGGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCggg                                                                                                                                         <  1:4344469/150‑1 (MQ=255)
             ctcctgggTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGtt                                                                                                                                       <  1:3516896/145‑1 (MQ=255)
                  gggTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTAccc                                                                                                                                  <  2:1525116/150‑1 (MQ=255)
                   ggTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCcg                                                                                                                                 >  2:3120559/1‑150 (MQ=255)
                       cccAAGCTCAAGCAGGCCAGTGCCGTCATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCGCCCACAACGGTATGCGTTCCCTCGTAGGAGCGGGTTTACCCGCgaa                                                                                                                             <  1:554039/150‑1 (MQ=18)
                       cccAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCGCCCACAACGGTATGCGTTCCCTCGTAGGAGCGGGTTTACCCGCgaa                                                                                                                             <  1:554037/150‑1 (MQ=21)
                        ccAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCgaag                                                                                                                            <  2:4940766/150‑1 (MQ=255)
                        ccAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCgaag                                                                                                                            >  2:524997/1‑150 (MQ=255)
                        ccAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCgaag                                                                                                                            >  2:525099/1‑150 (MQ=255)
                        ccAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGGTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCgaag                                                                                                                            >  2:4555320/1‑150 (MQ=255)
                         cAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGgc                                                                                                                                              >  1:5094723/1‑131 (MQ=255)
                         cAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGgc                                                                                                                                              <  2:5094723/131‑1 (MQ=255)
                         cAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCgaaga                                                                                                                           >  1:2701678/1‑150 (MQ=255)
                         cAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCgaaga                                                                                                                           >  2:395361/1‑150 (MQ=255)
                         cAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCGCCCACAACGGTATGCGTTCCCTCGTAGGAGCGGGttt                                                                                                                                      <  1:1004062/139‑1 (MQ=21)
                          aaGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCgaagaa                                                                                                                          <  1:3673760/150‑1 (MQ=255)
                            gCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGc                                                                                                                        <  2:1018017/150‑1 (MQ=255)
                              tCAAGCTGGCCAGTGCCGACTTGGTGGCTGAGCATCTGAAAGTGGATCTGGTGGATTTGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGGGGGTTTACCCGCGAAGAAGCCa                                                                                                                      <  1:5156194/150‑1 (MQ=38)
                              tCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTAccc                                                                                                                                  <  1:1660656/138‑1 (MQ=255)
                                aaGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGt                                                                                                                    >  2:4280383/1‑150 (MQ=255)
                                 tgccggCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTg                                                                                                                   >  2:941198/5‑150 (MQ=255)
                                 aGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTg                                                                                                                   >  1:416831/1‑150 (MQ=255)
                                    aGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTGGGTTGCAGGCCCGGCCTATTCGGGGGTAAACCCGCTCATACGTGGGAACGCTTACCGTTGTGGGCGCGGGTTTACCCGCTAAGAAGCCAGTGCAg                                                                                                                <  1:3023457/150‑1 (MQ=255)
                                      gCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCcgcg                                                                                                                               <  2:1459112/133‑1 (MQ=255)
                                      gCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGa                                                                                                              >  2:4596207/1‑150 (MQ=255)
                                      gCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGa                                                                                                              >  2:4596357/1‑150 (MQ=255)
                                      gCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCGCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCcgcg                                                                                                                               >  1:1459010/1‑133 (MQ=255)
                                       ccAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGaa                                                                                                             <  1:1069842/150‑1 (MQ=255)
                                              cGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCaa                                                                                                      >  2:4176180/1‑150 (MQ=255)
                                              cGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCaa                                                                                                      <  2:1932597/150‑1 (MQ=255)
                                              cGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCaa                                                                                                      <  2:2999735/150‑1 (MQ=255)
                                                aCATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAgg                                                                                                    <  1:1703882/150‑1 (MQ=255)
                                                        gCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGCAGAAGCCAGTGCAGGAATACCCAAGGCCTATCaa                                                                                            >  2:2204218/1‑150 (MQ=255)
                                                          tGAGCATCTGCAACTGGATCTGCTGGTTTTGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACt                                                                                          <  1:36224/150‑1 (MQ=255)
                                                          tGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGTTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACt                                                                                          <  1:36417/150‑1 (MQ=255)
                                                          tGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACt                                                                                          <  1:4802161/150‑1 (MQ=255)
                                                          tGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACt                                                                                          >  1:2201396/1‑150 (MQ=255)
                                                          tGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCACCTACGAGGGAACGCATAACGTTGTGGGCGCGGGGTTACCCGCGAAGAAGCCAGTGCAGGAATACCAAAGGCCTATCaaat                                                                                          >  1:1158355/1‑148 (MQ=255)
                                                           gAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGCACGCATACCGTGGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTg                                                                                         <  2:1682193/150‑1 (MQ=255)
                                                           gAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAAACCGCGCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGGTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCGaacaacag                                                                                         >  1:557898/1‑148 (MQ=255)
                                                                  tGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAg                                                                                                                >  2:1238387/1‑120 (MQ=255)
                                                                  tGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAg                                                                                                                <  1:1238387/120‑1 (MQ=255)
                                                                  tGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGa                                                                                  <  2:3386370/150‑1 (MQ=255)
                                                                  tGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGa                                                                                  <  1:4763379/150‑1 (MQ=255)
                                                                  tGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATAACGGTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGa                                                                                  >  1:474226/1‑150 (MQ=255)
                                                                    cAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTAGCGGGTAAACCCCCTCCTACGAGGGAACGCAAACCCTTGTGGGCGCGGGGTTTCCCCCGAAGAAGCCAGTGCAGGAAAACCCAACGCCTATCAACTGGCAATGatc                                                                                >  1:1896498/1‑148 (MQ=255)
                                                                     aaCTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAACCCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGcc                                                                               <  1:1486264/150‑1 (MQ=255)
                                                                      aCTGGATCTGCTGGTTTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGTTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCa                                                                              <  1:1834778/150‑1 (MQ=255)
                                                                          gATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAg                                                                          <  1:2740713/150‑1 (MQ=255)
                                                                          gATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATAAACAGGCATTGAGGCATTAg                                                                          >  1:3397562/1‑150 (MQ=255)
                                                                              tgctgGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGAt                                                                      >  2:2835575/1‑150 (MQ=255)
                                                                                ctgGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGAttt                                                                    >  2:3972652/1‑150 (MQ=255)
                                                                                    aTTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGATGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAAc                                                                <  1:4921682/150‑1 (MQ=255)
                                                                                      tAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACgg                                                              <  2:4796415/150‑1 (MQ=255)
                                                                                         tCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTTGGGATTTAAACggggg                                                           >  2:1535360/1‑148 (MQ=255)
                                                                                         tCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGcc                                                           <  1:3883834/150‑1 (MQ=255)
                                                                                           tGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGCAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCCTTGAGCCATTAGGGATTTAAACGGGCCtt                                                         >  2:2544490/1‑150 (MQ=255)
                                                                                           tGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTa                                                                           >  2:329554/1‑132 (MQ=255)
                                                                                              tGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTc                                                      <  2:69305/150‑1 (MQ=255)
                                                                                                  ggCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGccc                                                  <  1:1748657/150‑1 (MQ=255)
                                                                                                  ggCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGccc                                                  <  1:2855630/150‑1 (MQ=255)
                                                                                                     ccGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTa                                               <  1:323775/150‑1 (MQ=255)
                                                                                                     ccGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTa                                               >  2:4696903/1‑150 (MQ=255)
                                                                                                          ctctTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTg                                          <  1:3166081/150‑1 (MQ=255)
                                                                                                             ttCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGttt                                       >  1:1311129/1‑150 (MQ=255)
                                                                                                                gcgGGTAAACCCGCTCCTACGAGGGACCGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCAGTAGGGATTTAAACGGGCCTTTTCGCCCCTCGCCTGTTGTaa                                    >  2:830591/1‑150 (MQ=255)
                                                                                                                  gggTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAAtt                                  <  2:3267216/150‑1 (MQ=255)
                                                                                                                   ggTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTg                                 >  2:1835274/1‑150 (MQ=255)
                                                                                                                    gTAACCCCGCTCCTACGCGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTgg                                <  1:941198/150‑1 (MQ=255)
                                                                                                                    gTAACCCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGTTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTgg                                <  1:4581076/150‑1 (MQ=255)
                                                                                                                    gTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTgg                                <  1:4581047/150‑1 (MQ=255)
                                                                                                                    gTAAACCCGCGCCCACAACGGTATGCGTTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGc                                                            >  2:166451/1‑122 (MQ=17)
                                                                                                                     tAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTggg                               <  1:1486163/150‑1 (MQ=255)
                                                                                                                     tAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTggg                               <  1:4999072/150‑1 (MQ=255)
                                                                                                                     tAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTggg                               <  2:1311129/150‑1 (MQ=255)
                                                                                                                     tAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTGTACCCGCGACGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTggg                               <  1:4103339/150‑1 (MQ=255)
                                                                                                                     tAAACCCGCGCCCACAACGGTATGCGTTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTggg                               >  2:728264/1‑150 (MQ=21)
                                                                                                                       aaCCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGc                             >  2:733812/1‑150 (MQ=255)
                                                                                                                        aCCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTaaa                                                                 <  2:214772/113‑1 (MQ=255)
                                                                                                                         cccGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTg                           >  2:3840930/1‑150 (MQ=255)
                                                                                                                                  aCGAGGGAACGCATACCGTTGTGGGCGCGGGATTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATGAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTATGTGGGGCTGTGCGTCATc                  >  1:894109/1‑150 (MQ=255)
                                                                                                                                            gCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAAt        <  2:4102705/150‑1 (MQ=255)
                                                                                                                                            gCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAAt        <  2:3678247/150‑1 (MQ=255)
                                                                                                                                            gCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAAt        >  1:3695269/1‑150 (MQ=255)
                                                                                                                                               tACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGa     <  1:851548/150‑1 (MQ=255)
                                                                                                                                                 ccGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATgggta   >  1:2622277/1‑147 (MQ=255)
                                                                                                                                                  cGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTc  >  2:1884004/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGGCGCCTTGCCTGGGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTC  >  NC_002947/444433‑444728

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: