Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,212,523 G→A 8.2% intergenic (+39/+55) PP_4593 → / ← PP_4594 hypothetical protein/cystathionine gamma‑synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,212,5230GA 8.2% 206.8 / 7.4 73intergenic (+39/+55)PP_4593/PP_4594hypothetical protein/cystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base G (33/34);  new base A (3/3);  total (36/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 = 8.92e-01

GTAGCCGCCGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGC  >  NC_002947/5212374‑5212672
                                                                                                                                                     |                                                                                                                                                      
gtAGCCGCCGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGACATCTATGTAATCCACCGGCCTCTTAGCGGGAAAg                                                                                                                                                        >  2:1089210/1‑150 (MQ=255)
       ccgAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTc                                                                                                                                                 <  2:642461/150‑1 (MQ=255)
           aGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTac                                                                                                                                             <  1:1556017/150‑1 (MQ=255)
              cAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑Acc                                                                                                                                         >  1:225767/1‑150 (MQ=255)
                       cAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTc                                                                                                                                                 <  2:71313/134‑1 (MQ=255)
                       cAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTc                                                                                                                                                 >  1:71313/1‑134 (MQ=255)
                           gacgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTACTGAATCCACCGGCCCTTTCGCGGGTAAACCCGCTCCTACGACGGT‑GTAGGAGCgg                                                                                                                            <  1:1109753/150‑1 (MQ=37)
                            acgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCggg                                                                                                                           <  1:319165/150‑1 (MQ=255)
                             cgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGt                                                                                                                          <  1:553803/150‑1 (MQ=255)
                              gaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGtt                                                                                                                         <  2:422140/150‑1 (MQ=255)
                              gaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGtt                                                                                                                         >  1:756445/1‑150 (MQ=255)
                                       aaCCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTc                                                                                                                                                 >  1:1009099/1‑118 (MQ=255)
                                       aaCCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTc                                                                                                                                                 <  2:1009099/118‑1 (MQ=255)
                                          cGCCTGAGCAATGAAAGGTGGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑CCCGTCGTAGGAGCGGGTTTACCCGCGAAAg                                                                                                             <  2:756445/150‑1 (MQ=255)
                                                   aaGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCgg                                                                                                    >  2:1564329/1‑150 (MQ=255)
                                                   aaGCAAAGGTTGCTGATCAAGGACAGGATCGAACGGATCATTCATCCGCCCAAGGCCTGACAGGCATCTATTGAATCCCCCTGCCTCTTCGCTGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCgg                                                                                                    <  2:819825/150‑1 (MQ=38)
                                                     ggAAAGGTTGCTGATCAAGTACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCTACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGcc                                                                                                  <  1:714155/150‑1 (MQ=255)
                                                     ggAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAAGCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGcc                                                                                                  <  1:1089210/150‑1 (MQ=255)
                                                      gAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGccc                                                                                                 >  1:59194/1‑150 (MQ=255)
                                                       aaaGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCt                                                                                                >  1:301702/1‑150 (MQ=255)
                                                        aaGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTc                                                                                               <  2:1478191/150‑1 (MQ=255)
                                                            ttGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcacac                                                                                           >  2:237522/1‑150 (MQ=255)
                                                                tGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGTCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCg                                                                                       >  2:1622835/1‑150 (MQ=255)
                                                                tGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGcc                                                                                       >  1:1478477/1‑149 (MQ=255)
                                                                   tCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGT‑GTAGGAGCGGGCTTGCCCGCGaa                                                                                                               <  1:494871/123‑1 (MQ=21)
                                                                   tCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑GACGGTGTAGGAGCGGGCTTGCCCGCGaa                                                                                                               >  2:494871/1‑123 (MQ=21)
                                                                    cAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTc                                                                                               >  1:333968/1‑138 (MQ=255)
                                                                    cAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTc                                                                                               <  2:333968/138‑1 (MQ=255)
                                                                      aGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACTCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgc                                                                                 >  1:1320839/1‑150 (MQ=255)
                                                                       ggACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGc                                                                                                                                                   >  2:1213754/1‑84 (MQ=255)
                                                                       ggACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGc                                                                                                                                                   <  1:1213754/84‑1 (MQ=255)
                                                                       ggACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCt                                                                                                <  2:702054/134‑1 (MQ=255)
                                                                       ggACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCt                                                                                                >  1:702054/1‑134 (MQ=255)
                                                                        gACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑ACCGTCGTAGGAGc                                                                                                                              >  2:1550927/1‑103 (MQ=255)
                                                                        gACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑ACCGTCGTAGGAGc                                                                                                                              <  1:1550927/103‑1 (MQ=255)
                                                                          cAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCtt                                                                             >  1:966403/1‑150 (MQ=255)
                                                                           aaGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCg                                                                                                     >  2:564346/1‑125 (MQ=255)
                                                                           aaGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCg                                                                                                     <  1:564346/125‑1 (MQ=255)
                                                                             gATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTCCa                                                                          >  1:1147622/1‑150 (MQ=255)
                                                                                 gAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTAccc                                                                                                                    >  1:959991/1‑104 (MQ=255)
                                                                                 gAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTAccc                                                                                                                    <  2:959991/104‑1 (MQ=255)
                                                                                  aaCGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCttt                                                                            <  2:53234/143‑1 (MQ=255)
                                                                                  aaCGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCttt                                                                            >  1:53234/1‑143 (MQ=255)
                                                                                      gATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCAt                                                                 >  1:1282189/1‑150 (MQ=255)
                                                                                      gATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCAt                                                                 >  1:525819/1‑150 (MQ=255)
                                                                                       aTCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATc                                                                >  2:562142/1‑150 (MQ=255)
                                                                                          attcatCCGCCCAAATCCTTACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAg                                                             <  2:673537/150‑1 (MQ=255)
                                                                                          attcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcac                                                                                             >  1:1295429/1‑118 (MQ=255)
                                                                                          attcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcac                                                                                             <  2:1295429/118‑1 (MQ=255)
                                                                                          attcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAg                                                             <  1:562142/150‑1 (MQ=255)
                                                                                          attcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAg                                                             <  2:525819/150‑1 (MQ=255)
                                                                                              atCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGaaa                                                                                                              >  2:307132/1‑97 (MQ=255)
                                                                                              atCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGaaa                                                                                                              <  1:307132/97‑1 (MQ=255)
                                                                                              atCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTcc                                                         >  2:1322522/1‑150 (MQ=255)
                                                                                               tCCGCCCAAATCCTGACAGGCGTCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAgg                                                                                                            <  2:650629/98‑1 (MQ=255)
                                                                                               tCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAgg                                                                                                            >  1:650629/1‑98 (MQ=255)
                                                                                                ccGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTc                                                       <  2:1147622/150‑1 (MQ=255)
                                                                                                    ccAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGAcc                                                   >  1:215655/1‑150 (MQ=255)
                                                                                                     cAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGcc                                                                     <  2:572757/131‑1 (MQ=255)
                                                                                                     cAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGcc                                                                     >  1:572757/1‑131 (MQ=255)
                                                                                                      aaaCCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGAAAGCCCCCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAg                                                 <  1:509715/150‑1 (MQ=255)
                                                                                                       aaTCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAgg                                                <  2:225767/150‑1 (MQ=255)
                                                                                                                 ggCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAAt                                      <  2:1320839/150‑1 (MQ=255)
                                                                                                                     tCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTc                                  >  1:322302/1‑150 (MQ=255)
                                                                                                                        aTTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCt                                                                                                        <  1:178650/77‑1 (MQ=255)
                                                                                                                        aTTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCt                                                                                                        >  2:178650/1‑77 (MQ=255)
                                                                                                                              tCCACCGGCTCCGTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTAc                         <  1:887831/150‑1 (MQ=255)
                                                                                                                                 aCCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGa                      >  2:967085/1‑150 (MQ=255)
                                                                                                                                     gCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAAc                  <  2:322302/150‑1 (MQ=255)
                                                                                                                                     gCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTAGATCCCTACCGACAAc                  >  1:928606/1‑150 (MQ=255)
                                                                                                                                             gcgGGTAAACCCGCTCCTACGACGGT‑GTAGGAGCGGGCTTGCCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGcagcag          >  2:262677/1‑150 (MQ=37)
                                                                                                                                                   aaGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTtctc    <  2:59194/150‑1 (MQ=255)
                                                                                                                                                   aaGCCCGCTCCTAC‑ACCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGc                                                                                        <  2:1187035/66‑1 (MQ=255)
                                                                                                                                                   aaGCCCGCTCCTAC‑ACCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGc                                                                                        >  1:1187035/1‑66 (MQ=255)
                                                                                                                                                     gCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGcccc                                                                   <  2:1071656/85‑1 (MQ=255)
                                                                                                                                                     gCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGcccc                                                                   >  1:1071656/1‑85 (MQ=255)
                                                                                                                                                     gCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGc  <  1:1322522/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GTAGCCGCCGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGC  >  NC_002947/5212374‑5212672

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: ‑