Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 589,856 A→G 5.4% intergenic (‑1208/‑24) PP_0498 ← / → PP_0500 hypothetical protein/dTDP‑4‑rhamnose reductase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947589,8560AG 5.4% 227.8 / 2.0 75intergenic (‑1208/‑24)PP_0498/PP_0500hypothetical protein/dTDP‑4‑rhamnose reductase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (35/35);  new base G (2/2);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

GTCAGCCCATTCCTCGTATGTCGTAGGAAATTTCTGTCTGTGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTGGGGCCGAGGAGGACATCGCTTTCCTGGCACCGCGCCCGCCCGAGAACGGCTGGGCGCCC  >  NC_002947/589709‑590002
                                                                                                                                                   |                                                                                                                                                  
gTCAGCCCATTCCTCGTATGTCGTAGGAAATTTCTGTCTGTGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAAc                                                                                                                                                  >  2:1389683/1‑150 (MQ=255)
   aGCCCATTCCTCGTATGTCGTAGGAAATTTCTGTCTGTGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACctc                                                                                                                                               >  1:305943/1‑150 (MQ=255)
    gCCCATTCCTCGTATGTCGTAGGAAATTTCTGTCTGTGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACctct                                                                                                                                              >  2:1011887/1‑150 (MQ=255)
         ttCCTCGTATGTCGTAGGAAATTTCTGTCTGTGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGa                                                                                                                                         >  1:233451/1‑150 (MQ=255)
           ccTCGTATGTCGTAGGAAATTTCTGTCTGTGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATc                                                                                                                                       <  2:1294570/150‑1 (MQ=255)
            cTCGTATGTCGTAGGAAATTTCTGTCTGTGGCGGTATTGAGCGGTGAAACTGTGGAAGCACGCCCTGTGCCATCGTACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATct                                                                                                                                      <  2:284211/150‑1 (MQ=255)
             tCGTATGTCGTAGGAAATTTCTGTCTGTGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATctg                                                                                                                                     >  2:317211/1‑150 (MQ=255)
              cGTATGTCGTAGGAAATTTCTGTCTGTGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATctgc                                                                                                                                    <  1:1140696/150‑1 (MQ=255)
               gTATGTCGTAGGAAATTTCTGTCTGTGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATctgct                                                                                                                                   <  2:233451/150‑1 (MQ=255)
                       tAGGAAATTTCTGTCTGTGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTCTCTCATCAACCTCTCCGGATCTGCTGACTCtat                                                                                                                           <  1:1365648/150‑1 (MQ=255)
                             aTTTCTGTCTGTGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCtatgcgtat                                                                                                                     <  2:1231655/150‑1 (MQ=255)
                             aTTTCTGTCTGTGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCtatgcgtat                                                                                                                     <  1:297039/150‑1 (MQ=255)
                              tttCTGTCTGTGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCtatgcgtatg                                                                                                                    >  1:482463/1‑150 (MQ=255)
                                  tgtctgtGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCaa                                                                                                                                                   >  1:982570/1‑115 (MQ=255)
                                  tgtctgtGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCaa                                                                                                                                                   <  2:982570/115‑1 (MQ=255)
                                    tctgtGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTa                                                                                                              >  1:880900/1‑150 (MQ=255)
                                          gCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAAc                                                                                                                                                  <  1:182751/108‑1 (MQ=255)
                                          gCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAAc                                                                                                                                                  >  2:182751/1‑108 (MQ=255)
                                          gCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGa                                                                                                                                 <  1:1444649/125‑1 (MQ=255)
                                          gCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGa                                                                                                                                 >  2:1444649/1‑125 (MQ=255)
                                              tATTGCGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTgg                                                                                                    <  2:189234/150‑1 (MQ=255)
                                               aTTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGGAGGAAACCGCAAAATGACTCTTTTTTCTCATCATCCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTggg                                                                                                   <  2:895940/150‑1 (MQ=255)
                                                  gAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGggtgg                                                                                                >  2:704622/1‑150 (MQ=255)
                                                   aGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGggtggt                                                                                               >  1:1073390/1‑150 (MQ=255)
                                                       gTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCtat                                                                                                                           >  2:720870/1‑118 (MQ=255)
                                                       gTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCtat                                                                                                                           <  1:720870/118‑1 (MQ=255)
                                                         gAAACTGTGGAATCACGCCCTGTGGCATCGAACTGGTCCGGTGTTGGACTTCCCGGGCAGGAAAACGCAAATGGACTCTTTTTTCTCACCCACCTCTCCGCCTCTGCTCACTCTATGCGTATGCGCCTAATGCTGGTGGTtggtgggcac                                                                                         >  1:708488/1‑146 (MQ=255)
                                                            attgtgGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGcc                                                                                      <  2:1016368/148‑1 (MQ=255)
                                                              tgtgGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCtatgcgt                                                                                                                       <  1:405209/115‑1 (MQ=255)
                                                              tgtgGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCtatgcgt                                                                                                                       >  2:405209/1‑115 (MQ=255)
                                                                tgGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCg                                                                                  >  2:935294/1‑150 (MQ=255)
                                                                 ggAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCgg                                                                                 <  1:1389683/150‑1 (MQ=255)
                                                                  gAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCggg                                                                                <  2:288921/150‑1 (MQ=255)
                                                                    aTCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCa                                                                              <  1:704622/150‑1 (MQ=255)
                                                                       aCGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAgc                                                                           >  2:1049799/1‑150 (MQ=255)
                                                                       aCGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAgc                                                                           >  1:321112/1‑150 (MQ=255)
                                                                       aCGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAgc                                                                           >  1:393930/1‑150 (MQ=255)
                                                                         gCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAgcgc                                                                         <  2:393930/150‑1 (MQ=255)
                                                                              gtgGCATCGAACGGGGCCGGGGTTGGACTTCCCGGGCATGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGAtt                                                                    <  2:532354/150‑1 (MQ=255)
                                                                               tgGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTc                                                                   >  1:1120977/1‑150 (MQ=255)
                                                                               tgGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTc                                                                   <  2:264701/150‑1 (MQ=255)
                                                                                  cATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTACTGACTCTATGCGTATGCGCCTAATGCTGTTGGgtggtg                                                                                              <  1:824959/120‑1 (MQ=255)
                                                                                  cATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTACTGACTCTATGCGTATGCGCCTAATGCTGTTGGgtggtg                                                                                              >  2:824959/1‑120 (MQ=255)
                                                                                   aTCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGggtggt                                                                                               >  1:987146/1‑118 (MQ=255)
                                                                                   aTCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGggtggt                                                                                               <  2:987146/118‑1 (MQ=255)
                                                                                   aTCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCt                                                               <  2:1073390/150‑1 (MQ=255)
                                                                                    tCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCCTCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCtt                                                              <  2:995045/150‑1 (MQ=255)
                                                                                    tCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCtt                                                              >  2:1018343/1‑150 (MQ=255)
                                                                                    tCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCtt                                                              <  2:1120977/150‑1 (MQ=255)
                                                                                      gAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTgg                                                            <  2:321112/150‑1 (MQ=255)
                                                                                      gAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTgg                                                            >  2:1245907/1‑150 (MQ=255)
                                                                                           ggTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGtt                                                                                                      <  2:1558711/103‑1 (MQ=255)
                                                                                           ggTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGtt                                                                                                      >  1:1558711/1‑103 (MQ=255)
                                                                                            gTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTGGGGCCga                                                      <  1:317211/150‑1 (MQ=255)
                                                                                                ggTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTGGGGCCgaggag                                                  >  1:180669/1‑150 (MQ=255)
                                                                                                 gtgtTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTGGGGCCGaggagg                                                 <  2:305943/150‑1 (MQ=255)
                                                                                                      ggACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTGGGGCCGAGGAGGACATc                                            >  1:1221343/1‑150 (MQ=255)
                                                                                                       gACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCtat                                                                                                                           >  1:111909/1‑70 (MQ=255)
                                                                                                       gACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCtat                                                                                                                           <  2:111909/70‑1 (MQ=255)
                                                                                                       gACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTg                                                             <  1:475867/132‑1 (MQ=255)
                                                                                                       gACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTg                                                             >  2:475867/1‑132 (MQ=255)
                                                                                                       gACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTGGGGCCGAGGAGGACATCg                                           >  2:1144335/1‑150 (MQ=255)
                                                                                                       gACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTGGGGCCGAGGAGGACATCg                                           >  2:604705/1‑150 (MQ=255)
                                                                                                        aCTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTGGGGCCGAGGAGGACATCGc                                          >  2:849206/1‑150 (MQ=255)
                                                                                                        aCTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTGGGGCCGAGGAGGACATCGc                                          <  1:604705/150‑1 (MQ=255)
                                                                                                          ttCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCGACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGcccgggaa                                                                                <  2:1203441/110‑6 (MQ=255)
                                                                                                          ttCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCGACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGcccgggaa                                                                                >  1:1203887/1‑105 (MQ=255)
                                                                                                          ttCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCGACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGcccgggaa                                                                                >  1:1203441/1‑105 (MQ=255)
                                                                                                          ttCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCGACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGcccgggaa                                                                                <  2:1203887/110‑6 (MQ=255)
                                                                                                             ccGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTGGGGCCGAGGAGGACATCGCTTTcc                                     >  1:1389381/1‑150 (MQ=255)
                                                                                                             ccGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTGGGGCCGAGGAGGACATCGCTTTcc                                     >  2:454849/1‑150 (MQ=255)
                                                                                                               gggCAGGAAACCGCACAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGGGTATGCGCCTAATGCTGTTGGTTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTGGGGCCGAGGAGGACATCGCTTTCCTg                                   <  2:708488/150‑1 (MQ=255)
                                                                                                                ggCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGtggtgg                                                                                             >  2:235743/1‑91 (MQ=255)
                                                                                                                ggCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGtggtgg                                                                                             <  1:235743/91‑1 (MQ=255)
                                                                                                                    ggAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCATATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTTGTCTTGGGGCCGAGGAGGACATCGCTTTCATGGCAcc                              >  1:264873/1‑150 (MQ=255)
                                                                                                                                          tttCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTGGGGCCGAGGAGGACAt                                             >  1:1015086/1‑113 (MQ=255)
                                                                                                                                          tttCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTGGGGCCGAGGAGGACAt                                             <  2:1015086/113‑1 (MQ=255)
                                                                                                                                                atcaACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTGGGGCCGAGGAGGACATCGCTTTCCTGGCACCGCGCCCGCCCGAGAACGGCTGGGCGccc  >  1:178415/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTCAGCCCATTCCTCGTATGTCGTAGGAAATTTCTGTCTGTGGCGGTATTGAGCGGTGAAACTGTGGAATCACGCCCTGTGGCATCGAACGGGTCCGGTGTTGGACTTCCCGGGCAGGAAACCGCAAAATGACTCTTTTTTCTCATCAACCTCTCCGGATCTGCTGACTCTATGCGTATGCGCCTAATGCTGTTGGGTGGTGGCAATGCCCTCGGGCAAGCGCTGATTCGTCTTGGGGCCGAGGAGGACATCGCTTTCCTGGCACCGCGCCCGCCCGAGAACGGCTGGGCGCCC  >  NC_002947/589709‑590002

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