Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,660,236 A→C 13.5% intergenic (+38/‑191) metK → / → ligB methionine adenosyltransferase/DNA ligase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,660,2360AC13.5% 201.9 / 30.0 89intergenic (+38/‑191)metK/ligBmethionine adenosyltransferase/DNA ligase B
Reads supporting (aligned to +/- strand):  ref base A (37/40);  new base C (7/5);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-02

TTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCG  >  NC_002947/5660091‑5660385
                                                                                                                                                 |                                                                                                                                                     
ttCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCgg                                                                                                                                                   <  1:894415/150‑1 (MQ=255)
  cGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCa                                                                                                                                                 <  1:1655374/150‑1 (MQ=255)
  cGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCa                                                                                                                                                 <  2:561623/150‑1 (MQ=255)
   ggCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAg                                                                                                                                                <  1:397339/150‑1 (MQ=255)
    gCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAgg                                                                                                                                               <  2:1495772/150‑1 (MQ=255)
       gTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGc                                                                                                                                            <  2:197305/150‑1 (MQ=255)
           tCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCtttt                                                                                                                                        >  2:1796830/1‑150 (MQ=255)
           tCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCtttt                                                                                                                                        >  2:725321/1‑150 (MQ=255)
               cagcagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTtgt                                                                                                                                    <  1:1258568/150‑1 (MQ=255)
                       gaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTa                                                                                                                            <  2:1260233/150‑1 (MQ=255)
                               gCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCAttt                                                                                                                    >  2:82942/1‑150 (MQ=255)
                                   cgacACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCATACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCAt                                                                                                                <  1:725321/150‑1 (MQ=255)
                                   cgacACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCAt                                                                                                                <  1:82942/150‑1 (MQ=255)
                                    gacACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATc                                                                                                               >  1:1451766/1‑150 (MQ=255)
                                       accttcaccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGcc                                                                                                            >  1:1092229/1‑150 (MQ=255)
                                       accttcaccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGcc                                                                                                            >  1:1737055/1‑150 (MQ=255)
                                        ccttcaccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTCCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGccc                                                                                                           <  1:407160/150‑1 (MQ=255)
                                         cttcaccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCg                                                                                                          >  2:598287/1‑150 (MQ=255)
                                         cttcaccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCg                                                                                                          >  2:1608288/1‑150 (MQ=255)
                                          ttcaccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGt                                                                                                         <  1:781998/150‑1 (MQ=255)
                                              ccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAg                                                                                                                                                >  2:1085403/1‑107 (MQ=255)
                                              ccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAg                                                                                                                                                <  1:1085403/107‑1 (MQ=255)
                                              ccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGc                                                                                                                                  <  2:26435/121‑1 (MQ=255)
                                              ccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGc                                                                                                                                  >  1:26435/1‑121 (MQ=255)
                                              ccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGa                                                                                                     <  2:441277/150‑1 (MQ=255)
                                              ccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGa                                                                                                     <  1:272339/150‑1 (MQ=255)
                                              ccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACAAAAGCCCCTGCCGAGCGATCGGCAGGGGCttttgtctgg                                                                                                                                  >  2:1064688/1‑115 (MQ=255)
                                              ccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACAAAAGCCCCTGCCGAGCGATCGGCAGGGGCttttgtctgg                                                                                                                                  <  1:1064688/121‑7 (MQ=255)
                                              ccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTACGTGACGCTGCGGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGa                                                                                                     >  1:1892635/1‑150 (MQ=255)
                                               caccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTTCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGAt                                                                                                    >  1:1414885/1‑150 (MQ=255)
                                               caccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGAt                                                                                                    >  1:80795/1‑150 (MQ=255)
                                               caccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGAt                                                                                                    >  2:596411/1‑150 (MQ=255)
                                                   ttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAg                                                                                                <  1:1608288/150‑1 (MQ=255)
                                                      accTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGcc                                                                                                                                 <  2:576102/114‑1 (MQ=255)
                                                      accTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGcc                                                                                                                                 >  1:576102/1‑114 (MQ=255)
                                                        cTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGa                                                                                                 >  2:1083472/1‑144 (MQ=255)
                                                        cTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGa                                                                                                 <  1:1083472/144‑1 (MQ=255)
                                                          gggAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGa                                                                                                     >  1:728146/1‑138 (MQ=255)
                                                          gggAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGa                                                                                                     <  2:728146/138‑1 (MQ=255)
                                                                  accgaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTc                                                                                 >  2:231640/1‑150 (MQ=255)
                                                                  accgaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTc                                                                                 >  2:669315/1‑150 (MQ=255)
                                                                     gaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTg                                                                              <  2:99231/150‑1 (MQ=255)
                                                                       ccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCt                                                                            <  2:1158587/150‑1 (MQ=255)
                                                                             ccAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTGCGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGtct                                                                      >  2:1733199/1‑150 (MQ=255)
                                                                             ccAGTCGTTGCGTGACGCTGCCGCCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCCGCAGGGGCTTTTTTTTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGACCCCCGCGCTCGGTTCTTGATcccgcc                                                                      >  2:111488/1‑146 (MQ=255)
                                                                              cAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGtctc                                                                     <  2:1815074/150‑1 (MQ=255)
                                                                                gTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCgg                                                                   <  2:1737055/150‑1 (MQ=255)
                                                                                 tCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCggg                                                                  >  1:1628863/1‑150 (MQ=255)
                                                                                    ttGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCtt                                                               >  1:1613113/1‑150 (MQ=255)
                                                                                     tGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTg                                                              <  1:231640/150‑1 (MQ=255)
                                                                                       cGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCt                                                                                >  1:566933/1‑130 (MQ=255)
                                                                                       cGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCt                                                                                <  2:566933/130‑1 (MQ=255)
                                                                                         tGACGCTGCCGGCCTGTAAGCTTCCTGACGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGCGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGgaga                                                          <  1:493802/150‑1 (MQ=255)
                                                                                               tGCCGGCCTATAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGAGCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCCGCTCGAGCGCACCGCTAGTTCCCTGCCCTTGCTCGGGCGTTGAGGCGAGCa                                                    >  2:1181777/1‑150 (MQ=255)
                                                                                                  cGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTtgt                                                                                                                                    <  2:169129/67‑1 (MQ=255)
                                                                                                  cGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTtgt                                                                                                                                    >  1:169129/1‑67 (MQ=255)
                                                                                                  cGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCa                                                 <  1:1796830/150‑1 (MQ=255)
                                                                                                   ggCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAGGGGCtttcgt                                                                                                                                      <  1:1611552/64‑4 (MQ=21)
                                                                                                   ggCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAGGGGCtttcgt                                                                                                                                      >  2:1611552/1‑61 (MQ=18)
                                                                                                    gCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGc                                               >  1:1433996/1‑150 (MQ=255)
                                                                                                       tGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGc                                                                                                                                                  >  2:1300646/1‑48 (MQ=255)
                                                                                                       tGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGc                                                                                                                                                  <  1:1300646/48‑1 (MQ=255)
                                                                                                          aaGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAata                                         <  2:1613113/150‑1 (MQ=255)
                                                                                                             cTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGCACCCGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCGCCAGCGGAatagcc                                      >  2:444882/1‑147 (MQ=255)
                                                                                                             cTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCa                                      <  2:1892635/150‑1 (MQ=255)
                                                                                                             cTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCa                                      >  2:1024817/1‑150 (MQ=255)
                                                                                                              ttCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGgaga                                                          >  1:907542/1‑129 (MQ=255)
                                                                                                              ttCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGgaga                                                          <  2:907542/129‑1 (MQ=255)
                                                                                                              ttCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCAc                                     <  1:1024817/150‑1 (MQ=255)
                                                                                                              ttCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCAc                                     <  2:1685709/150‑1 (MQ=255)
                                                                                                                          cAGACGAAAGCCCCTGCCGAGCGATCGGCATGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACaa                              >  2:878736/1‑145 (MQ=255)
                                                                                                                          cAGACGAAAGCCCCTGCCGAGCGATCGGCATGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACaa                              <  1:878736/145‑1 (MQ=255)
                                                                                                                           acaaaaaaGCCCCTGCCGATCGCTCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGCCACCGCCTTGGCGACGAGCATCAGCGAAATATCACCTGACAACTgcac                        >  2:212584/6‑150 (MQ=38)
                                                                                                                             aaaaaaGCCCCTGCCGATCGCTCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTgcacgc                      >  1:333568/4‑150 (MQ=255)
                                                                                                                              cGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTgcacgca                     <  1:210751/150‑1 (MQ=255)
                                                                                                                                   gCCCCTGCCGATCGCTCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCa                                      >  2:1880792/1‑128 (MQ=255)
                                                                                                                                   gCCCCTGCCGATCGCTCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCa                                      <  1:1880792/128‑1 (MQ=255)
                                                                                                                                     cccTGCCGATCGCTCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGga                                                            >  2:473944/1‑104 (MQ=255)
                                                                                                                                     cccTGCCGATCGCTCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGga                                                            <  1:473944/104‑1 (MQ=255)
                                                                                                                                     cccTGCCGATCGCTCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCAccgcccg              >  1:540006/1‑150 (MQ=255)
                                                                                                                                     cccTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCAccgcccg              >  2:1653572/1‑150 (MQ=255)
                                                                                                                                      ccTGCCGATCGCTCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCAccgcccgc             >  1:1092400/1‑150 (MQ=255)
                                                                                                                                      ccTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTgcacgca                     <  2:884853/142‑1 (MQ=255)
                                                                                                                                      ccTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTgcacgca                     >  1:884853/1‑142 (MQ=255)
                                                                                                                                        tGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCtt                                                                               >  2:726058/1‑82 (MQ=255)
                                                                                                                                        tGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCtt                                                                               <  1:726058/82‑1 (MQ=255)
                                                                                                                                           cGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGAc                                >  1:1410001/1‑126 (MQ=255)
                                                                                                                                           cGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGAc                                <  2:1410001/126‑1 (MQ=255)
                                                                                                                                           cGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGcccgcccagc        <  1:1653572/150‑1 (MQ=255)
                                                                                                                                                 aTCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGc                                               <  1:805735/105‑1 (MQ=255)
                                                                                                                                                 aTCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGc                                               >  2:805735/1‑105 (MQ=255)
                                                                                                                                                 aTCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCg  >  1:203622/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCCGCCCCGCCACCCG  >  NC_002947/5660091‑5660385

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: