Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,458,940:1 +A 13.8% intergenic (+17/‑45) leuS → / → PP_4795 leucine‑‑tRNA ligase/lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,458,9401.A13.8% 303.5 / 48.8 94intergenic (+17/‑45)leuS/PP_4795leucine‑‑tRNA ligase/lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (41/40);  new base A (5/7);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

CCGCTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAA  >  NC_002947/5458793‑5459082
                                                                                                                                                    |                                                                                                                                              
ccGCTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCaa                                                                                                                                              >  2:1425020/1‑149 (MQ=255)
ccGCTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAt                                                                                                                                              >  2:2256148/1‑150 (MQ=255)
ccGCTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAt                                                                                                                                              >  2:1623501/1‑150 (MQ=255)
   cTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGa                                                                                                                                           <  2:535446/150‑1 (MQ=255)
           cgcgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCAGCCGCGATCGATGCGGGt                                                                                                                                   <  2:1235735/150‑1 (MQ=255)
             cgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTgg                                                                                                                                 <  1:1425020/150‑1 (MQ=255)
             cgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTgg                                                                                                                                 <  2:2297221/150‑1 (MQ=255)
             cgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTgg                                                                                                                                 <  2:1522439/150‑1 (MQ=255)
                agaagTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTggg                                                                                                                                <  1:1597758/148‑1 (MQ=255)
                agaagTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTggg                                                                                                                                >  2:1597758/1‑148 (MQ=255)
                      cGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAg                                                                                                                        >  1:546212/1‑150 (MQ=255)
                      cGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAg                                                                                                                        >  1:1745333/1‑150 (MQ=255)
                       gAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTgg                                                                                                                                 <  1:1465922/140‑1 (MQ=255)
                       gAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTgg                                                                                                                                 >  2:1465922/1‑140 (MQ=255)
                       gAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAgg                                                                                                                       >  2:2296777/1‑150 (MQ=255)
                        aaGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTTATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGc                                                                                                                      <  1:312832/150‑1 (MQ=255)
                           ccGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCa                                                                                                                   <  2:1445832/150‑1 (MQ=255)
                            cGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCaa                                                                                                                  >  1:2125894/1‑150 (MQ=255)
                              cTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCaaaa                                                                                                                <  1:391833/150‑1 (MQ=255)
                                gCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAgg                                                                                                              <  1:2256148/150‑1 (MQ=255)
                                 cccGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAggg                                                                                                             >  1:1528931/1‑150 (MQ=255)
                                  ccGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCTCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTgg                                                                                                                                 >  1:271544/1‑129 (MQ=255)
                                  ccGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCTCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTgg                                                                                                                                 <  2:271544/129‑1 (MQ=255)
                                  ccGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                            >  2:1585361/1‑150 (MQ=255)
                                   cGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAg                                                                                                           <  2:1745333/150‑1 (MQ=255)
                                   cGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAg                                                                                                           <  2:546212/150‑1 (MQ=255)
                                         aaCGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTgg                                                                                                                                 <  1:57517/122‑1 (MQ=255)
                                         aaCGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTgg                                                                                                                                 >  2:57517/1‑122 (MQ=255)
                                          aCGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCaacaac                                                                                                    >  1:225779/1‑150 (MQ=255)
                                           cGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCCCCCATATCCATTTCGGGTTGGCGGAACAGGCCCAAAAGGGAGCAacaaca                                                                                                   <  1:1623501/150‑1 (MQ=255)
                                            gagaATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAgg                                                                                                                       >  2:1245841/1‑129 (MQ=255)
                                            gagaATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAgg                                                                                                                       <  1:1245841/129‑1 (MQ=255)
                                            gagaATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGaacaac                                                                                                       >  1:1697896/1‑145 (MQ=255)
                                            gagaATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGcaac                                                                                                       <  2:1697896/145‑1 (MQ=255)
                                            gagaATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACAt                                                                                                  >  1:87563/1‑150 (MQ=255)
                                                   tCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                            <  2:678695/133‑1 (MQ=255)
                                                   tCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                            >  1:678695/1‑133 (MQ=255)
                                                   tCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCaaa                                                                                           <  2:1992333/150‑1 (MQ=255)
                                                         gcTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGACAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAAt                                                                                     >  2:891088/1‑150 (MQ=255)
                                                           ttgaTCGATGGCCTGACCATCCGCAAGGTTATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAAttt                                                                                   <  2:436919/147‑1 (MQ=255)
                                                           ttCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAAttt                                                                                   >  2:1207350/1‑150 (MQ=255)
                                                            tCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTg                                                                                  >  2:1619262/1‑150 (MQ=255)
                                                            tCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGCGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTg                                                                                  >  1:37140/1‑150 (MQ=255)
                                                             cATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGc                                                                                 <  2:1528931/150‑1 (MQ=255)
                                                              aTCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCt                                                                                <  1:1359630/150‑1 (MQ=255)
                                                               tCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGGAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTg                                                                               <  2:1592855/150‑1 (MQ=255)
                                                               tCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTg                                                                               <  1:1207350/150‑1 (MQ=255)
                                                               tCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTg                                                                               <  1:1890116/150‑1 (MQ=255)
                                                                 gATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCgg                                                                                                                             >  1:301874/1‑102 (MQ=255)
                                                                 gATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCgg                                                                                                                             <  2:301874/102‑1 (MQ=255)
                                                                  aTGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTCGCCGCTGCAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTa                                                                            >  2:2244737/1‑150 (MQ=255)
                                                                        tGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGc                                                                      <  1:1619262/150‑1 (MQ=255)
                                                                         gACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGcaa                                                                                                        >  2:1805701/1‑115 (MQ=255)
                                                                         gACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGcaa                                                                                                        <  1:1805701/115‑1 (MQ=255)
                                                                            cATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAATGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTgg                                                                  >  1:1404423/1‑150 (MQ=255)
                                                                            cATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTgg                                                                  <  2:225779/150‑1 (MQ=255)
                                                                            cATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTgg                                                                  >  1:201107/1‑150 (MQ=255)
                                                                                 gCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCAGTc                                                             >  1:2394310/1‑150 (MQ=255)
                                                                                  cAAGGTCATCGTGGTACCGGGCAAGATGTTAAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCa                                                            <  1:1785734/150‑1 (MQ=255)
                                                                                     ggTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCgg                                                                       <  2:2304166/136‑1 (MQ=255)
                                                                                     ggTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCgg                                                                       >  1:2304166/1‑136 (MQ=255)
                                                                                      gTCATCGTGGTACCAGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCt                                                        >  2:1956719/1‑150 (MQ=255)
                                                                                         aTCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAg                                                     >  1:895010/1‑150 (MQ=255)
                                                                                             tGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGcc                                                 <  2:1676056/150‑1 (MQ=255)
                                                                                               gTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGCGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTg                                               <  2:37140/150‑1 (MQ=255)
                                                                                                  ccGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCgg                                            >  1:1889988/1‑150 (MQ=255)
                                                                                                         aGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAg                                     <  2:1497785/150‑1 (MQ=255)
                                                                                                            tGGTCAACATCGTAGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCgg                                                               >  1:1430285/1‑121 (MQ=255)
                                                                                                            tGGTCAACATCGTAGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCgg                                                               <  2:1430285/121‑1 (MQ=255)
                                                                                                             ggTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGt                                           <  2:838465/140‑1 (MQ=255)
                                                                                                             ggTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGt                                           >  1:838465/1‑140 (MQ=255)
                                                                                                               tCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAg                                                                                                           >  1:1782025/1‑74 (MQ=255)
                                                                                                               tCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAg                                                                                                           <  2:1782025/74‑1 (MQ=255)
                                                                                                               tCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAAt                                                                                     >  2:1936052/1‑96 (MQ=255)
                                                                                                               tCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAAt                                                                                     <  1:1936052/96‑1 (MQ=255)
                                                                                                                 aaCATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCgg                             >  2:2183946/1‑150 (MQ=255)
                                                                                                                 aaCATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCgg                             <  2:201107/150‑1 (MQ=255)
                                                                                                                   cATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCa                           >  2:1661000/1‑150 (MQ=255)
                                                                                                                     tcgtcgCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCAcc                         >  1:2209452/1‑150 (MQ=255)
                                                                                                                      cgtcgCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAATGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCg                        <  2:1404423/150‑1 (MQ=255)
                                                                                                                        tcgCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAgc                                                    >  1:1383958/1‑120 (MQ=255)
                                                                                                                        tcgCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAgc                                                    <  2:1383958/120‑1 (MQ=255)
                                                                                                                                tgatgaAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGc                                                                      <  2:334727/94‑1 (MQ=255)
                                                                                                                                tgatgaAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGc                                                                      >  1:334727/1‑94 (MQ=255)
                                                                                                                                 gatgaAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAg             <  2:895010/150‑1 (MQ=255)
                                                                                                                                     aaaaCGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGa         <  1:1585361/150‑1 (MQ=255)
                                                                                                                                      aaaCGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTg                                               <  1:2202438/111‑1 (MQ=255)
                                                                                                                                      aaaCGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTg                                               >  2:2202438/1‑111 (MQ=255)
                                                                                                                                        aCGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCAAGCTGCGCGGCACGGGCGCCACCGAGCTGAGCg      >  1:200772/1‑150 (MQ=255)
                                                                                                                                          gCCCACCCGCAATCGATGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTgg                                                                              <  1:723454/76‑1 (MQ=37)
                                                                                                                                          gCCCACCCGCAATCGATGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTgg                                                                              >  2:723454/1‑76 (MQ=37)
                                                                                                                                            ccACCCGCAATCGATGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGaa  >  1:1869328/1‑150 (MQ=255)
                                                                                                                                            ccACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGaa  >  2:673128/1‑150 (MQ=255)
                                                                                                                                                ccGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCt                                                                                >  1:1849128/1‑68 (MQ=255)
                                                                                                                                                ccGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCt                                                                                <  2:1849128/68‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CCGCTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAA  >  NC_002947/5458793‑5459082

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