Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,907,550 G→A 7.2% G77G (GGC→GGT PP_1707 ← D‑isomer specific 2‑hydroxyacid dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,907,5500GA 7.2% 191.6 / 5.0 69G77G (GGC→GGTPP_1707D‑isomer specific 2‑hydroxyacid dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base G (30/34);  new base A (2/3);  total (32/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 = 7.97e-01

AACGAAGCGGTCACCGGTGCAGATCCCCCGGCACAAGCCGATCAGCAGGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATGCTGGCTGACGTAGGCCTCCTTGTCGGCCTGTTCGTAG  >  NC_002947/1907403‑1907694
                                                                                                                                                   |                                                                                                                                                
aacgaaGCGGTCACCGGTGCAGATCCCCCGGCACAAGCCGATCAGCAGGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGcc                                                                                                                                                >  1:822578/1‑150 (MQ=255)
      gCGGTCACCGGTGCAGATCCCCCGGCACAAGCCGATCAGCAGGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGCACCGACGcg                                                                                                                                          <  2:1274825/150‑2 (MQ=255)
       cGGTCACCGGTGCAGATCCCCCGGCACAAGCCGATCAGCAGGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCCCCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCa                                                                                                                                         <  1:959907/150‑1 (MQ=255)
        ggTCACCGGTGCAGATCCCCCGGCACAAGCCGATCAGCAGGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCAt                                                                                                                                        >  2:1168526/1‑150 (MQ=255)
               ggTGCAGATCCCCCGGCACAAGCCGATCAGCAGGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGc                                                                                                                                 <  1:169614/150‑1 (MQ=255)
                gTGCAGATCCCACGGCACAAGCCGATGAGCAGGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATGACGGGCATAGGCCAGGTATACCGCGTCGGTGCCGACGCCATTGACggag                                                                                                                                >  2:460446/1‑148 (MQ=255)
                    aGATCCCCCGGCACAAGCCGATCAGCAGGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCa                                                                                                                            >  2:459223/1‑150 (MQ=255)
                     gATCGCCCGGCACAAGCCGACCAGCAGGCCGATGCCACGATCGCCTGCAGCCCCGAAGAGCGCAGCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCAGAGGGCCGGATCACCGCGTCGAAGCCGACGCCATTGACGGCGATCAc                                                                                                                           >  2:145556/1‑150 (MQ=255)
                         ccccGGCACAAGCCGATCAGCAGGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCg                                                                                                                       >  1:658136/1‑150 (MQ=255)
                          cccGGCACAAGCCGATCAGCAGGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCCCGATGTTGCCGTCCGTCACCCGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGa                                                                                                                      <  1:460446/150‑1 (MQ=255)
                           ccGGCACAAGCCGATCAGCAGGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAg                                                                                                                     <  2:229938/150‑1 (MQ=255)
                            cGGCACAAGCCGATCAGCAGTCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAgg                                                                                                                    <  1:225112/150‑1 (MQ=255)
                               cacaAGCCGATCAGCAGGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTg                                                                                                                 <  1:660020/150‑1 (MQ=255)
                                  aaGCCGATCAGCAGGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGc                                                                                                              >  2:1024025/1‑150 (MQ=255)
                                     ccGATCAGCAGGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAgg                                                                                                           >  1:928571/1‑150 (MQ=255)
                                            gcagGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCa                                                                                                    >  1:731144/1‑150 (MQ=255)
                                            gcagGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCa                                                                                                    >  2:654724/1‑150 (MQ=255)
                                                gCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTGAGGTTGGGCAGGCGCGCCATCAgc                                                                                               >  1:686210/1‑150 (MQ=255)
                                                       gCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTgg                                                                                         <  2:1412301/150‑1 (MQ=255)
                                                           gATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCAAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGa                                                                                     >  1:1194497/1‑150 (MQ=255)
                                                            aTCGGCCACGTCGTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATACGCTGTGTCTAATGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGAt                                                                                    <  1:100828/150‑1 (MQ=255)
                                                                gCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCg                                                                                >  2:635526/1‑150 (MQ=255)
                                                                 cGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCTGGCGCGCCATCAGCGCCTGGCTGATGCCgg                                                                               <  2:595255/150‑1 (MQ=255)
                                                                    ctgccgcGGTGAGCGCACCGATGGTCGCCGTCCCCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCAATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGgtgt                                                                            <  1:453227/144‑1 (MQ=255)
                                                                     gTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGACGAGATCCCGGGCATATGCCAGGTCTACCGCGTGGGTGCCGACGCCAGTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGtgtg                                                                           <  1:674375/150‑1 (MQ=255)
                                                                        cTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGCACCGACGcgg                                                                                                                                         >  2:624187/1‑83 (MQ=255)
                                                                        cTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGCACCGACGcgg                                                                                                                                         <  1:624187/85‑3 (MQ=255)
                                                                              gAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGt                                                                  <  1:1024025/150‑1 (MQ=255)
                                                                               aGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTg                                                                 <  1:459223/150‑1 (MQ=255)
                                                                                gcgcACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCgg                                                                               <  2:1193700/135‑1 (MQ=255)
                                                                                gcgcACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCgg                                                                               >  1:1193700/1‑135 (MQ=255)
                                                                                    aCCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCAc                                                            <  1:941931/150‑1 (MQ=255)
                                                                                     ccGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATg                                                                                   >  2:1165955/1‑126 (MQ=255)
                                                                                     ccGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATg                                                                                   <  1:1165955/126‑1 (MQ=255)
                                                                                       gATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGc                                                                                                                                                 <  2:1382170/62‑1 (MQ=255)
                                                                                       gATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGc                                                                                                                                                 >  1:1382170/1‑62 (MQ=255)
                                                                                            tCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGaa                                                    >  1:1025032/1‑150 (MQ=255)
                                                                                                    aCCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCAc                                                            >  1:836482/1‑134 (MQ=255)
                                                                                                    aCCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCCAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCAc                                                            <  2:836482/134‑1 (MQ=255)
                                                                                                     ccTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAg                                           <  1:654724/150‑1 (MQ=255)
                                                                                                      cTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGAt                                                               <  1:1310052/129‑1 (MQ=255)
                                                                                                      cTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGAt                                                               >  2:1310052/1‑129 (MQ=255)
                                                                                                      cTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGc                                          >  1:1321307/1‑150 (MQ=255)
                                                                                                            gCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATgctg                                    <  1:635526/150‑1 (MQ=255)
                                                                                                             ccGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAACATTGCCAGCATgctgg                                   <  2:1321307/150‑1 (MQ=255)
                                                                                                               gcgATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATgctggct                                 >  1:874181/1‑150 (MQ=255)
                                                                                                               ccgATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATgctggct                                 <  1:606304/149‑1 (MQ=255)
                                                                                                                     aCGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATGCTGGCTGACGTa                           <  1:1193194/150‑1 (MQ=255)
                                                                                                                       gggCATAGGCCAGGTCTACAGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATGCTGGCTGACGTAgg                         >  1:45013/1‑150 (MQ=255)
                                                                                                                        ggCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATGCTGGCTGACGTactg                        <  1:1342895/150‑4 (MQ=255)
                                                                                                                        ggCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATGCTGGCTGACGTAgg                         <  2:1254724/149‑1 (MQ=255)
                                                                                                                        ggCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATGCTGGCTGACGTAgg                         >  1:1254724/1‑149 (MQ=255)
                                                                                                                        ggCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATGCTGGCTGACGTAGGc                        <  2:822578/150‑1 (MQ=255)
                                                                                                                          cATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGGTCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATGCTGGCTGACGTAGGcct                      <  2:1025032/150‑1 (MQ=255)
                                                                                                                                   ggTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGgtgt                                                                            >  1:535662/1‑87 (MQ=255)
                                                                                                                                   ggTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGgtgt                                                                            <  2:535662/87‑1 (MQ=255)
                                                                                                                                    gTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATGCTGGCTGACGTAGGCCTCCTTGTCGGc            >  2:800255/1‑150 (MQ=255)
                                                                                                                                     tCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCa                                                                                                 >  2:1281907/1‑64 (MQ=255)
                                                                                                                                     tCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCa                                                                                                 <  1:1281907/64‑1 (MQ=255)
                                                                                                                                      cTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATGCTGGCTGACGTAGGCCTCCTTGTCGGCCt          <  2:928571/150‑1 (MQ=255)
                                                                                                                                       tACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATGCTGGCTGACGTa                           >  2:21215/1‑132 (MQ=255)
                                                                                                                                       tACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATGCTGGCTGACGTa                           <  1:21215/132‑1 (MQ=255)
                                                                                                                                         tggcgTCGGCACCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCa                                         >  2:404415/3‑116 (MQ=255)
                                                                                                                                         tggcgTCGGCACCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCa                                         <  1:404415/114‑1 (MQ=255)
                                                                                                                                         ccGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATgctggcaata                              <  1:202388/127‑5 (MQ=255)
                                                                                                                                         ccGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATgctggcaata                              >  2:202388/1‑123 (MQ=255)
                                                                                                                                         ccGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATGCTGGCTGACGTAGGCCTCCTTGTGGGCATGtt       >  1:40133/1‑150 (MQ=255)
                                                                                                                                           gcgTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATGCTGGCTGACGTAGGCCTCCTTGTCGGCCTGTTCg     >  2:170268/1‑150 (MQ=255)
                                                                                                                                              tCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATGCTGGCTGACGTAGGCCTCCTTGTCGGCCTGTTCgtag  <  2:584242/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
AACGAAGCGGTCACCGGTGCAGATCCCCCGGCACAAGCCGATCAGCAGGCCGATGGCAAGATCGGCGACGTCCTCGGTGAGCGCACCGATGGTCGCCGTCACCTGGATGCCGCGATCACGGGCATAGGCCAGGTCTACCGCGTCGGTGCCGACGCCATTGACGGCGATCACTTCGAGGTTGGGCAGGCGCGCCATCAGCGCCTGGCTGATGCCGGTGTGGCCGCCGGTGATCACGCCGCGAATATTGCCAGCATGCTGGCTGACGTAGGCCTCCTTGTCGGCCTGTTCGTAG  >  NC_002947/1907403‑1907694

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