Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 922,760 G→C 5.0% Q91E (CAG→GAG)  PP_0803 ← protein secretion ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947922,7600GC 5.0% 251.6 / 3.4 79Q91E (CAG→GAG) PP_0803protein secretion ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (43/32);  new base C (2/2);  total (45/34)
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.67e-01

GGTCAGGAGAGTCGCGCAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCACCTCGAACCAG  >  NC_002947/922616‑922908
                                                                                                                                                |                                                                                                                                                    
ggtcaggAGAGTCGCGCAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCa                                                                                                                                                 <  1:1915394/150‑1 (MQ=255)
 gtcaggAGAGTCGCGCAGGCTCTCCGGTAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAg                                                                                                                                                >  2:839891/1‑150 (MQ=255)
 gtcaggAGAGTCGCGCAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAg                                                                                                                                                >  2:692390/1‑150 (MQ=255)
     ggAGAGTCGCGCAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCAcc                                                                                                                                            <  1:1822658/150‑1 (MQ=255)
      gagaGTCGCGCAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCt                                                                                                                                           >  2:1242295/1‑150 (MQ=255)
           tCGCGCTGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCcg                                                                                                                                      >  2:1128682/1‑150 (MQ=255)
             gcgcAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCcgcg                                                                                                                                    >  2:1441501/1‑150 (MQ=255)
                    ctctCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCCCGCTGGACGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGAcc                                                                                                                             <  1:1475314/150‑1 (MQ=255)
                     tctcCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCa                                                                                                                            >  1:375275/1‑150 (MQ=255)
                     tctcCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCa                                                                                                                            >  2:489860/1‑150 (MQ=255)
                       tcCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAgg                                                                                                                          >  1:1117528/1‑150 (MQ=255)
                       tcCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAgg                                                                                                                          <  1:1060582/150‑1 (MQ=255)
                       tcCAGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCAGGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAgg                                                                                                                          >  2:1205626/1‑150 (MQ=255)
                           ggAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCAc                                                                                                                      >  2:515147/1‑150 (MQ=255)
                                 tCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCGGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCgg                                                                                                                >  1:106336/1‑150 (MQ=255)
                                 tCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCgg                                                                                                                >  1:105817/1‑150 (MQ=255)
                                       gCTTGCCGATCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTTGCCGCGCTCGACCAGGTCACCTTCGGCCACGc                                                                                                          >  2:22074/1‑150 (MQ=255)
                                           gCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCAt                                                                                                      <  1:692390/150‑1 (MQ=255)
                                           gCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCAt                                                                                                      <  1:1242295/150‑1 (MQ=255)
                                                 cACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTgcgc                                                                                                <  2:443627/150‑1 (MQ=255)
                                                  aCCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGcgcg                                                                                               >  1:1115724/1‑150 (MQ=255)
                                                    ctctGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTGGGCGTGGATTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGcgac                                                                                             >  2:927902/1‑150 (MQ=255)
                                                     tctGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAgg                                                                                                                          >  2:1024969/1‑120 (MQ=255)
                                                     tctGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAgg                                                                                                                          <  1:1024969/120‑1 (MQ=255)
                                                        gcctgcAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATa                                                                                         <  2:12664/150‑1 (MQ=255)
                                                        gcctgcAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATa                                                                                         <  2:645837/150‑1 (MQ=255)
                                                        gcctgcAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATa                                                                                         <  1:839891/150‑1 (MQ=255)
                                                          ctgcAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCa                                                                                                                            <  2:1736823/113‑1 (MQ=255)
                                                          ctgcAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCa                                                                                                                            >  1:1736823/1‑113 (MQ=255)
                                                          ctgcAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATAcc                                                                                       <  2:1117528/150‑1 (MQ=255)
                                                                  gggCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGa                                                                               >  1:618452/1‑150 (MQ=255)
                                                                   ggCCTGGCTGGCCCTGGCGGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGaa                                                                              >  2:25437/1‑150 (MQ=255)
                                                                   ggCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGTTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACAATACCGCCCTCGaa                                                                              >  2:788991/1‑150 (MQ=255)
                                                                    gCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAAc                                                                             <  2:1881929/150‑1 (MQ=255)
                                                                    gCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCAGGTCTCGACCAGGTCACCTTCGGCCACGCGCATCTGCGCGACGATACCGCCCTCGAAc                                                                             >  1:1454145/1‑150 (MQ=255)
                                                                     ccTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGcc                                                                                    <  2:588146/142‑1 (MQ=255)
                                                                     ccTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGcc                                                                                    >  1:588146/1‑142 (MQ=255)
                                                                      ctggctggCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGa                                                                           >  1:1356541/1‑150 (MQ=255)
                                                                       tggctggcCCTGGCGGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGAc                                                                          <  2:106336/150‑1 (MQ=255)
                                                                       tggctggcCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGAc                                                                          <  2:105817/150‑1 (MQ=255)
                                                                       tggctggcCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGAc                                                                          <  1:1441501/150‑1 (MQ=255)
                                                                                cTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCt                                                                                                                    <  2:57599/99‑1 (MQ=255)
                                                                                cTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCt                                                                                                                    >  1:57599/1‑99 (MQ=255)
                                                                                cTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGATTTAACGCGACCAGGTCACCTTCGGCCCCGCCCCCCTGCGCGACAATGCCGCCCCCGAACGACTGCATGAcc                                                                 >  1:798302/1‑150 (MQ=255)
                                                                                  ggCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTg                                                               <  2:605128/150‑1 (MQ=255)
                                                                                   gCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGc                                                              >  1:1762071/1‑150 (MQ=255)
                                                                                        gcgGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACATTAAGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTcg                                                         >  1:1453896/1‑150 (MQ=255)
                                                                                                  gCCTCGCTTTCACCCACGCTGGGAGCCGTCTTGGGCGGGCCGAGCTGGGCCAGCACCTGGCCGCCCTCGCCCAGGTCACCTTCGCCCGCGATGCTCCGACCGAGGATACCGCCCTCCAACGCCTGACGGGCCTGTTCCCGCGCACTTgca                                               >  2:910567/1‑148 (MQ=255)
                                                                                                   ccTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGaa                                              <  1:25437/150‑1 (MQ=255)
                                                                                                    cTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAAt                                             >  2:1291035/1‑150 (MQ=255)
                                                                                                    cTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAAt                                             >  1:56053/1‑150 (MQ=255)
                                                                                                      cGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGa                                           <  1:489860/150‑1 (MQ=255)
                                                                                                         tttCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTc                                                                                                                        <  2:135438/70‑1 (MQ=255)
                                                                                                         tttCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTc                                                                                                                        >  1:135438/1‑70 (MQ=255)
                                                                                                         tttCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCt                                        >  1:1645051/1‑150 (MQ=255)
                                                                                                             aCCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTc                                                                                 <  2:1800872/105‑1 (MQ=255)
                                                                                                             aCCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTc                                                                                 >  1:1800872/1‑105 (MQ=255)
                                                                                                             aCCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGcc                                    >  2:1436922/1‑150 (MQ=255)
                                                                                                              cccACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGa                                           <  2:1447357/142‑1 (MQ=255)
                                                                                                              cccACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGa                                           >  1:1447357/1‑142 (MQ=255)
                                                                                                              cccACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCt                                   >  1:342968/1‑150 (MQ=255)
                                                                                                                  cGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGc                               <  2:1645051/150‑1 (MQ=255)
                                                                                                                       gAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTg                          >  2:42558/1‑150 (MQ=255)
                                                                                                                               cTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGTCCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATCCCACCCTAGCAAGAATGACTGAAGTGCTAATGCGAATTTAGGATTACTGTTTCTTTTCCGGTGGATACCTc                  <  2:1454145/150‑1 (MQ=255)
                                                                                                                               cTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTc                  >  1:1198535/1‑150 (MQ=255)
                                                                                                                                ttGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCCCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCCGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCg                 >  2:272927/1‑150 (MQ=255)
                                                                                                                                    tcgggtcgAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCTGCCACGCTCATCTGCGCGGCGATAACGCCCTCGAACGACTGAAATACCTGCTCTCGCGAAGTTGGAATGACCTTGCCTGTGCCGGGGGGTACCTCGGTcg             >  1:307081/1‑149 (MQ=255)
                                                                                                                                    tcgggtcgAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACGTTGCCTGTGCAGGTGGATAGCTCggcga             >  2:1354872/1‑147 (MQ=255)
                                                                                                                                    tcgggtcgAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCa             >  1:1038146/1‑150 (MQ=255)
                                                                                                                                     ctggtcgAGCGCGGCCAGGTGCTGGCCCAGCTCGACCCGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTg                          <  2:686452/134‑1 (MQ=37)
                                                                                                                                     ctggtcgAGCGCGGCCAGGTGCTGGCCCAGCTCGACCCGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTg                          >  1:686452/3‑136 (MQ=37)
                                                                                                                                     ctggtcgAGCGCGGCCAGGTGCTGGCCCAGCTCGACCATGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTcgc                                                        <  1:55183/104‑1 (MQ=25)
                                                                                                                                     ctggtcgAGCGCGGCCAGGTGCTGGCCCAGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTcgc                                                        >  2:55183/3‑106 (MQ=25)
                                                                                                                                     cgggtcgAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCGCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCAc            <  2:508414/150‑1 (MQ=255)
                                                                                                                                        gtcgAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCg                 >  2:485866/1‑142 (MQ=255)
                                                                                                                                        gtcgAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCg                 <  1:485866/142‑1 (MQ=255)
                                                                                                                                            aGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCACCTCGAAc     <  2:375275/150‑1 (MQ=255)
                                                                                                                                               tGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGCAAGAATGGATGGACTGCTATCACGCTCTTAGTATGCCCTATCCTGTGCCGGTGAAGACCTCGTTCACCTCGAACCAg  <  2:1453896/150‑1 (MQ=255)
                                                                                                                                               tGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCACCTCGAACCAg  <  1:1354872/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGTCAGGAGAGTCGCGCAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGCGCGGTACTTGGCCTCGCTTTCACCCACGCTGGAAGCCGTCTTGGTCGGGTCGAGCTGGGCCAGCACCTGGCCGCGCTCGACCAGGTCACCTTCGGCCACGCTCATCTGCGCGACGATACCGCCCTCGAACGACTGAATGACCTGCTCTCGCGAACTTGGAATGACCTTGCCTGTGCCGGTGGATACCTCGGTCACCTCGAACCAG  >  NC_002947/922616‑922908

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