Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,875,404 T→C 5.4% L29P (CTG→CCG)  cobS → adenosylcobinamide‑GDP ribazoletransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,875,4040TC 5.4% 279.0 / 3.9 93L29P (CTG→CCG) cobSadenosylcobinamide‑GDP ribazoletransferase
Reads supporting (aligned to +/- strand):  ref base T (40/48);  new base C (3/2);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

CCAGGGAGCAGTTGCTGCAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAGGGTACGCCGTCGCCGCTGCATGCGGCTTTGCTGCTGACGTTGT  >  NC_002947/1875259‑1875553
                                                                                                                                                 |                                                                                                                                                     
ccAGGGAGCAGTTGCTGCAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCg                                                                                                                                                   >  2:1026419/1‑150 (MQ=255)
  aGGTAGCAGTTGCTGCAGGGCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGACCGAATCGCACTGCGAGGTCGCCTTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGGCggg                                                                                                                                                 >  1:718666/1‑150 (MQ=255)
       gCAGTTGCTGCAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGc                                                                                                                                            <  2:1799329‑M2/150‑2 (MQ=255)
          gTTGCTGCAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGccgg                                                                                                                                         >  1:942025‑M2/1‑146 (MQ=255)
          gTTGCTGCAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGccgg                                                                                                                                         <  1:1291472‑M2/150‑5 (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:292423/150‑1 (MQ=255)
                    gTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGa                                                                                                                               <  1:508529/150‑1 (MQ=255)
                     tCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAg                                                                                                                              >  1:520130/1‑150 (MQ=255)
                        aGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTa                                                                                                                           <  2:1394211/150‑1 (MQ=255)
                        aGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTa                                                                                                                           <  1:1202961/150‑1 (MQ=255)
                        aGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTa                                                                                                                           >  1:1394211/1‑150 (MQ=255)
                          gTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAgg                                                                                                                         <  1:1728227/150‑1 (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:1103916‑M2/144‑38 (MQ=255)
                                                 gCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGCATGCCCGGCAGGCTCAccggcaggctgctgagaaactgcaaggcgatccagaa                                                                                                        >  2:134779‑M2/1‑107 (MQ=255)
                                                 gCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGCATGCCCGGCAGGCTCAccggcaggctgctgagaaactgcaaggcgatccagaa                                                                                                        <  1:134779‑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                                                                                      >  1:1356928/1‑150 (MQ=255)
                                                             tGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGgctgc                                                                                      <  2:442283/150‑1 (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:74185/1‑150 (MQ=255)
                                                                   tGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCg                                                                                <  1:1130987/150‑1 (MQ=255)
                                                                    gggCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCgg                                                                               <  1:1078859/150‑1 (MQ=255)
                                                                     ggCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGc                                                                              <  2:1356928/150‑1 (MQ=255)
                                                                     ggCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGc                                                                              <  1:1820989/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                                                                    <  2:1924013/138‑1 (MQ=255)
                                                                                           agaTGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCtgt                                                                    >  1:1924013/1‑138 (MQ=255)
                                                                                            gaTGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTtgc                                                                          <  1:813857/131‑1 (MQ=255)
                                                                                            gaTGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTtgc                                                                          >  2:813857/1‑131 (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                                                                           >  1:1940541/1‑118 (MQ=255)
                                                                                                        cTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTtg                                                                           <  2:1940541/118‑1 (MQ=255)
                                                                                                         tGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCtgtg                                                                   >  1:369276/1‑125 (MQ=255)
                                                                                                         tGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCtgtg                                                                   <  2:369276/125‑1 (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                                                                 >  1:279663/1‑106 (MQ=255)
                                                                                                                              agcagcCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGt                                                                 <  2:279663/106‑1 (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                      <  2:1056431/142‑1 (MQ=255)
                                                                                                                                     tGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAGGGTACGCCGTCGCCGCTGCATGc                      >  1:1056431/1‑142 (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)
                                                                                                                                                 tGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAGGGTACGCCGTCGCCGCTGCATGCGTCTTTGCTGCTGACGTtgt  <  1:1799329/150‑1 (MQ=255)
                                                                                                                                                 tGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAGGGTACGCCGTCGCCGCTGCATGCGGCTTTGCTGCTGACGTtgt  <  1:1342847/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCAGGGAGCAGTTGCTGCAGGTCGAGGTTGGCCATGGCGCGTTGATGCGGCTGGTGCCGGGTGGCGATGGGCAACTGGTCGAGGCGCCCTGAGATGTTGCCGTTCTGGATCGCCTTGCAGTTTCTCAGCAGCCTGCCGGTGAGCCTGCCGGGCATGCCGGCACCGCGTGAGGTAGGGCGTTCGCTGCTGTACTACCCGCTGGTCGGGCTGCTGTTCGGCTTGCTGCTGTGGTTGGCCAGTCATCTGCTACAGGGTACGCCGTCGCCGCTGCATGCGGCTTTGCTGCTGACGTTGT  >  NC_002947/1875259‑1875553

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: