Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,459,609 C→A 5.5% intergenic (+19/‑21) PP_4795 → / → holA lipoprotein/DNA polymerase III subunit delta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,459,6090CA 5.5% 290.8 / 3.2 90intergenic (+19/‑21)PP_4795/holAlipoprotein/DNA polymerase III subunit delta
Reads supporting (aligned to +/- strand):  ref base C (42/43);  new base A (2/3);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

GCCCCAACTGGATGAGCTGCAGCGTAAAGCAGACGAGCGCGCCCAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAGCCGCCGACGCCATCCGCAATGCCGCGCGTCA  >  NC_002947/5459461‑5459757
                                                                                                                                                    |                                                                                                                                                    
gCCCCAACTGGATGAGCTGCAGCGTAAAGCAGACGAGCGCGCCCAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCg                                                                                                                                                     <  1:724368/150‑1 (MQ=255)
 ccccAACTGGATGAGCTGCAGCGTAAAGCAGACGAGCGCGCCCAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCgg                                                                                                                                                    <  1:1635122/150‑1 (MQ=255)
  cccAACTGGATGAGCTGCAGCGTAAAGCAGACGAGCGCGCCCAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCggg                                                                                                                                                   <  2:2191371/150‑1 (MQ=255)
     aaCTGGATGAGCTGCAGCGTAAAGCAGACGAGCGCGCCCAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGgcgc                                                                                                                                                >  2:1352960/1‑150 (MQ=255)
          gATGAGCTGCAGCGTAAAGCAGACGAGCGCGCCCAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGCGGGGCGCCCGAGGGCGcccctc                                                                                                                                           >  2:897808/1‑148 (MQ=255)
          gATGAGCTGCAGCGTAAAGCAGACGAGCGCGCCCAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGc                                                                                                                                           <  1:106510/150‑1 (MQ=255)
            tGAGCTGCAGCGTAAAGCAGACGAGCGCGCCCAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGCGGGGCGCCCGAGGGCGctcctcgc                                                                                                                                         <  1:780017/150‑5 (MQ=255)
               gCTGCAGCGTAAAGCAGACGAGCGCGCCCAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGCGGGGCGCCCGAGGGCGcccctcgcggg                                                                                                                                      <  1:897808/150‑8 (MQ=255)
               gCTGCAGCGTAAAGCAGACGAGCGCGCCCAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGtt                                                                                                                                      >  2:2287488/1‑150 (MQ=255)
                 tGCAGCGTAAAGCAGACGAGCGCGCCCAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTcc                                                                                                                                    >  2:2071615/1‑150 (MQ=255)
                    aGCGTAAAGCAGACGAGCGCGCCCAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCa                                                                                                                                 <  1:2071615/150‑1 (MQ=255)
                             cAGACGAGCGCGCCCAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCg                                                                                                                        <  2:288352/150‑1 (MQ=255)
                                  gAGCGCGCCCAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCTGAGACGCCTCAGGAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGccccc                                                                                                                   <  2:1636733/150‑1 (MQ=255)
                                          ccAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACt                                                                                                           >  1:1847325/1‑150 (MQ=255)
                                          ccAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACt                                                                                                           >  2:1240040/1‑150 (MQ=255)
                                          ccAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGAGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACt                                                                                                           >  2:555980/1‑150 (MQ=255)
                                             aGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTcaa                                                                                                        >  2:1678031/1‑150 (MQ=255)
                                                cGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCACCGCCCAACTcaacaa                                                                                                     >  2:2295786/1‑150 (MQ=255)
                                                 gAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGa                                                                                                                              <  2:415591/124‑1 (MQ=255)
                                                 gAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGa                                                                                                                              >  1:415591/1‑124 (MQ=255)
                                                 gAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAg                                                                                                    >  1:769244/1‑150 (MQ=255)
                                                  aaGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGc                                                                                                   >  1:1808356/1‑150 (MQ=255)
                                                     cTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCAcc                                                                                                >  2:2116175/1‑150 (MQ=255)
                                                         cgcgCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCa                                                                                            >  1:91923/1‑150 (MQ=255)
                                                          ttgcTGGCGGCGGGCAGCCTCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCaa                                                                                           <  1:1745512/148‑1 (MQ=255)
                                                          gcgcTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGGCCCCGCCCAACTCAACAAGCACCTGCaa                                                                                           >  1:1072363/1‑150 (MQ=255)
                                                             cTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGc                                                                                                                       <  1:200299/119‑1 (MQ=255)
                                                             cTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGc                                                                                                                       >  2:200299/1‑119 (MQ=255)
                                                              tGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGgcg                                                                                       >  2:1898452/1‑150 (MQ=255)
                                                                 aGGCTGCCCGCCGCCATCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCTCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGccc                                                                                    <  1:209188/150‑1 (MQ=255)
                                                                    cTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTgg                                                                                 <  1:1898452/150‑1 (MQ=255)
                                                                      gCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCa                                                                               >  1:65741/1‑150 (MQ=255)
                                                                        ccgccgCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCAcc                                                                             >  1:821739/1‑150 (MQ=255)
                                                                        ccgccgCCAGCAAGCCGAGACGCATCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCCATGGGCGCCCTCGGGCGCCCCGCCTGTCCCCCATGACACTACCGCTCCTGCGGGCCAACAAGCACCTGCCAGGCGCCCTCGCAcc                                                                             >  2:2301729/1‑150 (MQ=255)
                                                                         cgccgcCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCg                                                                            <  1:1352960/150‑1 (MQ=255)
                                                                         cgccgcCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAACTACCGGGCAACTACGCCCGCCAGGAGCGCCCTCGGACGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACCAGCGCCTGCAAGGCTCCCTGGCACCg                                                                            <  2:1412913/150‑1 (MQ=255)
                                                                              ccAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTa                                                                       >  2:915544/1‑150 (MQ=255)
                                                                               cAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGcc                                                                                                                                          <  2:1352620/82‑1 (MQ=255)
                                                                               cAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGcc                                                                                                                                          >  1:1352620/1‑82 (MQ=255)
                                                                                   aaGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTgg                                                                  <  2:91923/150‑1 (MQ=255)
                                                                                   aaGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCCGTCTACgggg                                                                  >  1:1054866/1‑150 (MQ=255)
                                                                                             gCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGcga                                                        <  1:555980/150‑1 (MQ=255)
                                                                                             gCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGcga                                                        >  2:663148/1‑150 (MQ=255)
                                                                                               cTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGc                                                                                <  1:2054092/124‑1 (MQ=255)
                                                                                               cTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGc                                                                                >  2:2054092/1‑124 (MQ=255)
                                                                                                tCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGcgacga                                                     <  2:1442241/150‑1 (MQ=255)
                                                                                                   gCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGAccc                                                  >  2:2032655/1‑150 (MQ=255)
                                                                                                    cAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCg                                                 >  1:365351/1‑150 (MQ=255)
                                                                                                    cAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCg                                                 <  1:2287488/150‑1 (MQ=255)
                                                                                                           cACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCt                                                                        <  1:384725/120‑1 (MQ=255)
                                                                                                           cACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCt                                                                        >  2:384725/1‑120 (MQ=255)
                                                                                                           cACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCtgt                                          <  1:204226/150‑1 (MQ=255)
                                                                                                           cACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCtgt                                          <  2:65741/150‑1 (MQ=255)
                                                                                                            aCTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCACCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCtgtg                                         <  1:2295786/150‑1 (MQ=255)
                                                                                                             cTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGc                                        <  1:663148/150‑1 (MQ=255)
                                                                                                              tGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGcc                                       <  1:2238589/150‑1 (MQ=255)
                                                                                                                gAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAg                                     <  2:821739/150‑1 (MQ=255)
                                                                                                                   gTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGaa                                  <  1:915544/150‑1 (MQ=255)
                                                                                                                     aCCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAgc                                <  1:1678031/150‑1 (MQ=255)
                                                                                                                      ccGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACCCGACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAgcc                               <  1:155876/150‑1 (MQ=255)
                                                                                                                       cGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAgccg                              >  1:1338083/1‑150 (MQ=255)
                                                                                                                       cGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAgccg                              >  1:2150459/1‑150 (MQ=255)
                                                                                                                        gggCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAgccgc                             <  1:1440061/150‑1 (MQ=255)
                                                                                                                        gggCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAgccgc                             >  2:212504/1‑150 (MQ=255)
                                                                                                                           cAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTATACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAGccgccgc                          >  1:1681668/1‑149 (MQ=255)
                                                                                                                            aaCTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTAc                                                                      <  2:2133617/105‑1 (MQ=255)
                                                                                                                            aaCTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTAc                                                                      >  1:2133617/1‑105 (MQ=255)
                                                                                                                            aaCTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAGCCGCCGAc                         <  1:2326771/150‑1 (MQ=255)
                                                                                                                             aCTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGTGTCGACCCGCTGCTGTGCCAGGAAGCCGCCGACg                        <  2:115398/150‑1 (MQ=255)
                                                                                                                              cTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTcaacaa                                                                                                     >  1:1478094/1‑72 (MQ=255)
                                                                                                                              cTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTcaacaa                                                                                                     <  2:1478094/72‑1 (MQ=255)
                                                                                                                              cTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCgacgac                                                    >  1:702163/1‑121 (MQ=255)
                                                                                                                              cTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCgacgac                                                    <  2:702163/121‑1 (MQ=255)
                                                                                                                              cTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAGCCGCCGACGc                       >  1:2215406/1‑150 (MQ=255)
                                                                                                                              cTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAGCCGCCGACGc                       >  1:1845789/1‑150 (MQ=255)
                                                                                                                              cTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAGCCGCCGACGc                       >  2:1547026/1‑150 (MQ=255)
                                                                                                                               tAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAGCCGCCGACGcc                      >  1:1662501/1‑150 (MQ=255)
                                                                                                                                 aGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAgcc                               >  2:649230/1‑139 (MQ=255)
                                                                                                                                 aGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAgcc                               <  1:649230/139‑1 (MQ=255)
                                                                                                                                    ccGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCCGTGCCAGGAAGCCGCCGACGCCATCCg                 >  2:1373130/1‑150 (MQ=255)
                                                                                                                                       cgAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAGCCGCCGACGCCATCCGCaa              <  1:2116175/150‑1 (MQ=255)
                                                                                                                                           gggCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAg                                     <  2:1884988/123‑1 (MQ=255)
                                                                                                                                           gggCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAg                                     >  1:1884988/1‑123 (MQ=255)
                                                                                                                                           gggCGCCCGAGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCa                                      <  2:1713643/122‑1 (MQ=255)
                                                                                                                                           gggCGCCCGAGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCa                                      >  1:1713643/1‑122 (MQ=255)
                                                                                                                                             gcgcCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAGCCGCCGACGCCATCCGCAATGCcgc        <  2:1808356/150‑1 (MQ=255)
                                                                                                                                              cgcCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAGCCGCCGACGCCATCCGCAATGCcgcg       >  1:1703059/1‑150 (MQ=255)
                                                                                                                                                 ccTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAGCCGCCGACGCCATCCGCAATGCCGCGCGt    <  2:2215406/150‑1 (MQ=255)
                                                                                                                                                  cTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAGCCGCCGACGCCATCCGCAATGCCGCGCGTc   <  1:1240040/150‑1 (MQ=255)
                                                                                                                                                   tCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAGCCGCCGACGCCATCCGCAATGCCGCGCGTCa  <  1:1373130/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCCCAACTGGATGAGCTGCAGCGTAAAGCAGACGAGCGCGCCCAAGCCGAAGCTGCCGCGCTGGAGGCTGCCCGCCGCCAGCAAGCCGAGACGCCTCAGCAATCGCCACTGGAAGTACCGGGCAACTAAGCCCGCGAGGGGCGCCCTCGGGCGCCCCGCCTGTTCCCCATGAAACTCGCCCCCGCCCAACTCAACAAGCACCTGCAAGGCGCCCTGGCACCGGTCTACGTGGTCAGTGGCGACGACCCGCTGCTGTGCCAGGAAGCCGCCGACGCCATCCGCAATGCCGCGCGTCA  >  NC_002947/5459461‑5459757

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