Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,875,399 G→A 5.0% V27V (GTG→GTA) ‡ cobS → adenosylcobinamide‑GDP ribazoletransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,875,3990GA 5.0% 304.4 / 3.6 99V27V (GTG→GTA) ‡cobSadenosylcobinamide‑GDP ribazoletransferase
Reads supporting (aligned to +/- strand):  ref base G (43/51);  new base A (3/2);  total (46/53)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

TGGCTTGCCCAGGGAGCAGTTGCTGCAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAGGGTACGCCGTCGCCGCTGCATGCGGCTTTGCTGCT  >  NC_002947/1875251‑1875545
                                                                                                                                                    |                                                                                                                                                  
tgGCTTGCCCAGGGAGCAGTTGCTGAAGCTCTAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGa                                                                                                                                                   <  1:494032/150‑1 (MQ=255)
  gCTTGCCCAGGGAGCAGTTGCTGCAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGTCGGTGAGc                                                                                                                                                 <  1:206450/150‑1 (MQ=255)
  gCTTGCCCAGGGAGCAGTTGCTGCAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGc                                                                                                                                                 >  2:1728227/1‑150 (MQ=255)
   cTTGCCCAGGGAGCAGTTGCTGCAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGcc                                                                                                                                                >  1:930185/1‑150 (MQ=255)
        ccAGGGAGCAGTTGCTGCAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCg                                                                                                                                           >  2:1026419/1‑150 (MQ=255)
          aGGTAGCAGTTGCTGCAGGGCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGACCGAATCGCACTGCGAGGTCGCCTTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGGCggg                                                                                                                                         >  1:718666/1‑150 (MQ=255)
               gCAGTTGCTGCAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGc                                                                                                                                    <  2:1799329‑M2/150‑2 (MQ=255)
                  gTTGCTGCAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGccgg                                                                                                                                 <  1:1291472‑M2/150‑5 (MQ=255)
                  gTTGCTGCAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGccgg                                                                                                                                 >  1:942025‑M2/1‑146 (MQ=255)
                   ttGCTGCAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGccggc                                                                                                                                <  1:1549291‑M2/150‑6 (MQ=255)
                      ctgcAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCAcc                                                                                                                             >  2:1820989/1‑150 (MQ=255)
                       tgcAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACcg                                                                                                                            >  2:1078859/1‑150 (MQ=255)
                            gTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGa                                                                                                                       <  1:508529/150‑1 (MQ=255)
                            gTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGa                                                                                                                       <  1:292423/150‑1 (MQ=255)
                             tCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAg                                                                                                                      >  1:520130/1‑150 (MQ=255)
                                aGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTa                                                                                                                   >  1:1394211/1‑150 (MQ=255)
                                aGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTa                                                                                                                   <  1:1202961/150‑1 (MQ=255)
                                aGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTa                                                                                                                   <  2:1394211/150‑1 (MQ=255)
                                  gTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAgg                                                                                                                 <  1:1728227/150‑1 (MQ=255)
                                   ttGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGa                                                                                                                                                   <  2:1167221/115‑1 (MQ=255)
                                   ttGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGa                                                                                                                                                   >  1:1167221/1‑115 (MQ=255)
                                       ccATGGCGCGTTGATGCGGCTGGTGCCGGGTCGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGtt                                                                                                            >  2:1196295/1‑150 (MQ=255)
                                       ccATGGCGCGTTGATGCGGATGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGGATGCCGGCACCGCGTGAGGTAGGGCGtt                                                                                                            >  1:1440722/1‑150 (MQ=255)
                                            gcgcGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCgctg                                                                                                       <  1:1196295/150‑1 (MQ=255)
                                              gcgTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGCATGCCCGGCAGGCTCAccggcaggctgctgagaaactgcaaggcgatc                                                                                                     >  1:1259414‑M2/1‑118 (MQ=255)
                                                         gCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGc                                                                                          >  2:590067/1‑150 (MQ=255)
                                                         gCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGCATGCCCGGCAGGCTCAccggcaggctgctgagaaactgcaaggcgatccagaa                                                                                                >  2:1103916‑M2/1‑107 (MQ=255)
                                                         gCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGCATGCCCGGCAGGCTCAccggcaggctgctgagaaactgcaaggcgatccagaa                                                                                                <  1:134779‑M2/144‑38 (MQ=255)
                                                         gCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGCATGCCCGGCAGGCTCAccggcaggctgctgagaaactgcaaggcgatccagaa                                                                                                >  2:134779‑M2/1‑107 (MQ=255)
                                                         gCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGCATGCCCGGCAGGCTCAccggcaggctgctgagaaactgcaaggcgatccagaa                                                                                                <  1:1103916‑M2/144‑38 (MQ=255)
                                                            ggTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTgg                                                                                       <  1:1026419/150‑1 (MQ=255)
                                                               gCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCg                                                                                    <  1:1139087/150‑1 (MQ=255)
                                                                  gggTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGgc                                                                                 <  2:1291472/150‑1 (MQ=255)
                                                                   ggTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGgct                                                                                <  1:95794/150‑1 (MQ=255)
                                                                    gTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGAAACCACGTGAGGTATTGCTTGCGCTGCTGTACTACCCGCTGGTCGGgctg                                                                               >  2:1287994/1‑150 (MQ=255)
                                                                     tGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGgctgc                                                                              <  2:442283/150‑1 (MQ=255)
                                                                     tGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGgctgc                                                                              >  1:1356928/1‑150 (MQ=255)
                                                                         gATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGt                                                                                                             >  2:615462/1‑115 (MQ=255)
                                                                         gATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGt                                                                                                             <  1:615462/115‑1 (MQ=255)
                                                                         gATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGtt                                                                          >  1:1278513/1‑150 (MQ=255)
                                                                          aTGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTATGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTc                                                                         <  2:1258686/150‑1 (MQ=255)
                                                                           tGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCg                                                                        <  1:1130987/150‑1 (MQ=255)
                                                                           tGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCg                                                                        >  1:74185/1‑150 (MQ=255)
                                                                            gggCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCgg                                                                       <  1:1078859/150‑1 (MQ=255)
                                                                             ggCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGc                                                                      <  1:1820989/150‑1 (MQ=255)
                                                                             ggCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGc                                                                      <  2:1356928/150‑1 (MQ=255)
                                                                              gCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACTCGCTGGTCGGGCTGCTGTTCGGCt                                                                     <  2:560691/150‑1 (MQ=255)
                                                                               cAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGa                                                                                                                       >  1:496740/1‑99 (MQ=255)
                                                                               cAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGa                                                                                                                       <  2:496740/99‑1 (MQ=255)
                                                                               cAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGACGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCtt                                                                    <  1:1319510/150‑1 (MQ=255)
                                                                                aaCTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCAc                                                                                                                              >  1:226197/1‑91 (MQ=255)
                                                                                aaCTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCAc                                                                                                                              <  2:226197/91‑1 (MQ=255)
                                                                                aaCTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGCGCATGCCGGCACCGCGTGAGGTAGGGCGTTGACTGCTGTATTACCCGCTGGTCGGGCTGCTGTTCGGATtg                                                                   >  1:916279/1‑150 (MQ=255)
                                                                                 aCTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTtgc                                                                  >  1:1341794/1‑150 (MQ=255)
                                                                                   tGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTtgctg                                                                >  1:1524630/1‑150 (MQ=255)
                                                                                     gTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTgctgct                                                              >  2:1757156/1‑150 (MQ=255)
                                                                                     gTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTgctgct                                                              >  2:862820/1‑150 (MQ=255)
                                                                                          ggCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGt                                                         <  1:185015/150‑1 (MQ=255)
                                                                                           gcgcCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGtt                                                        <  2:1524630/150‑1 (MQ=255)
                                                                                            cgcCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCCGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTg                                                       <  1:1696054/150‑1 (MQ=255)
                                                                                             gcgcTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTTCCCGCGGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTgg                                                      <  2:718666/147‑1 (MQ=255)
                                                                                             gcCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTgg                                                      <  1:1201649/150‑1 (MQ=255)
                                                                                               ccTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGcc                                                    >  1:1622195/1‑150 (MQ=255)
                                                                                                cTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTGGTTGTTGGACTACTCGCTGGTCGCGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCa                                                   <  2:520130/150‑1 (MQ=255)
                                                                                                cTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCa                                                   >  1:1000865/1‑150 (MQ=255)
                                                                                                 tGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAg                                                  >  1:1269040/1‑150 (MQ=255)
                                                                                                   agaTGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCtgt                                                            >  1:1924013/1‑138 (MQ=255)
                                                                                                   agaTGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCtgt                                                            <  2:1924013/138‑1 (MQ=255)
                                                                                                    gaTGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTtgc                                                                  >  2:813857/1‑131 (MQ=255)
                                                                                                    gaTGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTtgc                                                                  <  1:813857/131‑1 (MQ=255)
                                                                                                        ttGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTg                                           >  1:1262795/1‑150 (MQ=255)
                                                                                                         tGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGc                                          <  1:1130186/150‑1 (MQ=255)
                                                                                                           ccGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCtgt                                                            <  1:1809710/130‑1 (MQ=255)
                                                                                                           ccGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCtgt                                                            >  2:1809710/1‑130 (MQ=255)
                                                                                                               tCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAgg                                    <  2:1262795/150‑1 (MQ=255)
                                                                                                                cTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTtg                                                                   <  2:1940541/118‑1 (MQ=255)
                                                                                                                cTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTtg                                                                   >  1:1940541/1‑118 (MQ=255)
                                                                                                                 tGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCtgtg                                                           <  2:369276/125‑1 (MQ=255)
                                                                                                                 tGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCtgtg                                                           >  1:369276/1‑125 (MQ=255)
                                                                                                                  ggATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCCGCTACAGGGTa                                 >  1:1451614/1‑150 (MQ=255)
                                                                                                                   gATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAGGGTAc                                <  1:518294/150‑1 (MQ=255)
                                                                                                                          ttGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAGGGTACGCCGTcg                         >  2:1342847/1‑150 (MQ=255)
                                                                                                                               gTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCAt                                              <  2:905728/124‑1 (MQ=255)
                                                                                                                               gTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCAt                                              >  1:905728/1‑124 (MQ=255)
                                                                                                                               gTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAGGGTACGCCGTCGCCGCt                    <  1:862820/150‑1 (MQ=255)
                                                                                                                                      agcagcCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGt                                                         <  2:279663/106‑1 (MQ=255)
                                                                                                                                      agcagcCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGt                                                         >  1:279663/1‑106 (MQ=255)
                                                                                                                                       gcagcCTGCCGGTGAGCCTGCCGTGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAGGGTACGCCGTCGCCGCTGCATGCgg            <  2:1340330/150‑1 (MQ=255)
                                                                                                                                           ccTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAGGGTACGCCGTCGCCGCTGCATGCGGCttt        <  1:1287994/150‑1 (MQ=255)
                                                                                                                                             tGCCGGTGAGCCTGCCGGGCATGccggcaggctgctgagaaactgcaaggcgatccagaacggcaacatctcagggcgcctcgaccagttgcccatcgccacccggcaccagccgcatcaacgcgccatggccaacctcgacctgcagca      <  1:420315‑M2/150‑128 (MQ=255)
                                                                                                                                             tGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAGGGTACGCCGTCGCCGCTGCATGc              >  1:1056431/1‑142 (MQ=255)
                                                                                                                                             tGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAGGGTACGCCGTCGCCGCTGCATGc              <  2:1056431/142‑1 (MQ=255)
                                                                                                                                              gCCGGTGAGCCTGCCGGGCATGccggcaggctgctgagaaactgcaaggcgatccagaacggcaacatctcagggcgcctcgaccagttgcccatcgccacccggcaccagccgcatcaacgcgccatggccaacctcgacctgcagcaa     >  2:420315‑M2/1‑22 (MQ=255)
                                                                                                                                                 ggTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTc                                                                                     <  1:1114654/67‑1 (MQ=255)
                                                                                                                                                 ggTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTc                                                                                     >  2:1114654/1‑67 (MQ=255)
                                                                                                                                                 ggTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAGGGTACGCCGTTGCCGCTGCATGCGGCTTTgctgct  <  2:74185/150‑1 (MQ=255)
                                                                                                                                                 ggTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAGGGTACGCCGTCGCCGCTGCATGCGGCTTTgctgct  >  2:1108982/1‑150 (MQ=255)
                                                                                                                                                  gTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAggg                                                                                                                >  2:204081/1‑39 (MQ=255)
                                                                                                                                                  gTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAggg                                                                                                                <  1:204081/39‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGGCTTGCCCAGGGAGCAGTTGCTGCAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAGGGTACGCCGTCGCCGCTGCATGCGGCTTTGCTGCT  >  NC_002947/1875251‑1875545

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: