Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 3,021,860 C→A 8.2% L210M (CTG→ATG)  ppnKB → inorganic polyphosphate/ATP‑NAD kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome3,021,8600CA 8.2% 234.7 / 8.8 73L210M (CTG→ATG) ppnKBinorganic polyphosphate/ATP‑NAD kinase
Reads supporting (aligned to +/- strand):  ref base C (30/37);  new base A (3/3);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

CCTATAACAAGTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGATCGAAATCAAACTGTCTAATAAGAAAATCAAAACAGTCAAACTAAAG  >  BS168_reference_genome/3021726‑3021994
                                                                                                                                      |                                                                                                                                      
ccTATAACAAGTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGAc                                                                                                                                    <  8:127017/139‑1 (MQ=255)
   ataACAACTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACAGAAAGCTGACCCt                                                                                                                                 >  2:62026/1‑139 (MQ=255)
     aaCAAGTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTc                                                                                                                               >  4:235231/1‑139 (MQ=255)
        aaGTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCGATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTTCCGGACACTCGGCTCGGCATTTGTCCTCAGCTCTGACCGAAAGCTGACCATTCgtg                                                                                                                            <  1:276006/139‑1 (MQ=255)
        aaGTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCgtg                                                                                                                            <  5:39522/139‑1 (MQ=255)
        aaGTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCgtg                                                                                                                            <  3:235231/139‑1 (MQ=255)
        aaGTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCgtg                                                                                                                            <  3:31082/139‑1 (MQ=255)
        aaGTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCgtg                                                                                                                            <  6:257710/139‑1 (MQ=255)
        aaGTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCgtg                                                                                                                            <  7:240822/139‑1 (MQ=255)
        aaGTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCgtg                                                                                                                            <  7:32580/139‑1 (MQ=255)
        aaGTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCgtg                                                                                                                            <  7:85876/139‑1 (MQ=255)
        aaGTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCgtg                                                                                                                            <  8:232272/139‑1 (MQ=255)
        aaGTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGACCTCAGCTCTGACCGAAAGCTGACCCTTCgtg                                                                                                                            <  4:275404/139‑1 (MQ=255)
        aaGTCTGTTGCAGGAGCTGGTGCTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCgtg                                                                                                                            <  2:266776/139‑1 (MQ=255)
         aGTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCgtgt                                                                                                                           >  1:106506/1‑139 (MQ=255)
          gTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGCCCGAAAGCTGACCCTTCGtgtg                                                                                                                          >  3:31112/1‑139 (MQ=255)
          gTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGtgtg                                                                                                                          >  8:97013/1‑139 (MQ=255)
          gTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGtgtg                                                                                                                          >  4:261012/1‑139 (MQ=255)
          gTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGtgtg                                                                                                                          >  3:212772/1‑139 (MQ=255)
          gTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGtgtg                                                                                                                          >  3:131899/1‑139 (MQ=255)
                 gCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAg                                                                                                                    <  6:207736/138‑1 (MQ=255)
                       gCTGTTGTTGATCCATTGCTTCCATTCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTAAAGGACACTCGGCTAGCCATATTTCCTCAGCTCTGACCGAAAGCAGACCCTTCGTGTGGTGAAGGACGGc                                                                                                              >  8:51152/1‑138 (MQ=255)
                          gttgttGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGGAGGACGGCAACg                                                                                                          <  1:62026/139‑1 (MQ=255)
                           ttgttgATCCATTGCTTCCATGCATGCAGGTATCAGAGCTAGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACg                                                                                                          <  5:170666/138‑1 (MQ=255)
                            tgttgATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGaa                                                                                                        >  1:227062/1‑139 (MQ=255)
                                     aTTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAat                                                                                                >  5:73092/1‑138 (MQ=255)
                                     aTTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAat                                                                                                >  8:89525/1‑138 (MQ=255)
                                      ttGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAatca                                                                                              <  2:106506/139‑1 (MQ=255)
                                            ccATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAatca                                                                                              <  1:186982/133‑1 (MQ=255)
                                            ccATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAatca                                                                                              >  2:186982/1‑133 (MQ=255)
                                            ccATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGc                                                                                         <  3:190876/138‑1 (MQ=255)
                                            ccATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCt                                                                                        >  5:89917/1‑139 (MQ=255)
                                            ccATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGGGTGGTGCAGGACGGCAACGAACACCCAATCATCGGc                                                                                         <  4:241117/138‑1 (MQ=255)
                                                             agagCTTGCATCCCTTAACAATACTACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACt                                                                        <  6:240163/138‑1 (MQ=255)
                                                             agagCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTc                                                                       <  1:5342/139‑1 (MQ=255)
                                                                    gCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGCCCCTTCGTGTGGTGCAGGCCGGCAACGCACACCCAATCATCGGCTTGGACAATGACGCACTCAGCAc                                                                  >  4:89149/1‑137 (MQ=255)
                                                                    gCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGAc                                                                                   >  2:106583/1‑120 (MQ=255)
                                                                    gCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGAc                                                                                   <  1:106583/120‑1 (MQ=255)
                                                                    gCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTCAGCAc                                                                  >  6:136210/1‑137 (MQ=255)
                                                                              aCAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGc                                                                                         <  8:157817/104‑1 (MQ=255)
                                                                              aCAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGc                                                                                         >  7:157817/1‑104 (MQ=255)
                                                                                    ataCCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGAt                                                <  3:261012/139‑1 (MQ=255)
                                                                                    ataCCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGAt                                                <  7:97013/139‑1 (MQ=255)
                                                                                    ataCCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGAt                                                >  7:151099/1‑139 (MQ=255)
                                                                                    ataCCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGAt                                                <  2:227652/139‑1 (MQ=255)
                                                                                    ataCCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTGGGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGAt                                                >  4:81681/1‑139 (MQ=255)
                                                                                    ataCCTACAGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGAt                                                <  4:23251/139‑1 (MQ=255)
                                                                                          accGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGCCCCTTCGTGTGGTGCAGGACGGCACCGAACCCCCAAGCATCGGCTTGGACAATGAAGCACACAGCACGAGAAATGTGAAAACGATCGaaa                                           >  6:258360/1‑138 (MQ=255)
                                                                                          accGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCgt                                                                                                                             >  8:124701/1‑56 (MQ=255)
                                                                                          accGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCgt                                                                                                                             <  7:124701/56‑1 (MQ=255)
                                                                                          accGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCgt                                                                                                                             <  7:128600/56‑1 (MQ=255)
                                                                                          accGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCgt                                                                                                                             >  8:128600/1‑56 (MQ=255)
                                                                                          accGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGc                                                                                                              >  8:32749/1‑71 (MQ=255)
                                                                                          accGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGc                                                                                                              <  7:32749/71‑1 (MQ=255)
                                                                                          accGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGa                                                                                    >  7:294116/1‑97 (MQ=255)
                                                                                          accGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGa                                                                                    <  8:294116/97‑1 (MQ=255)
                                                                                          accGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGATCGAAAt                                          >  2:33909/1‑139 (MQ=255)
                                                                                          accGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTGGGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGATCGAAAt                                          >  2:83193/1‑139 (MQ=255)
                                                                                                cacTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGCCCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGATCGAAATCAAACt                                    >  4:219396/1‑139 (MQ=255)
                                                                                                   tCGGATCGCAATTTGTCCTAAGATCTGCCCGAAAGCTGACCCGTAGTTTGGTGCAGGACGGCAAAGAACACCCAATCATCGGCTTGGACAAAGAAGCACTCAGAAAGAGAAATGTGAAAACGATCGAAATCAAACttt                                  >  8:277690/1‑136 (MQ=255)
                                                                                                     ggCTCGCCATTTGTCCTCAGCTCTGACCGAAAGATGAc                                                                                                                                    <  7:116726/38‑1 (MQ=255)
                                                                                                     ggCTCGCCATTTGTCCTCAGCTCTGACCGAAAGATGAc                                                                                                                                    >  8:120611/1‑38 (MQ=255)
                                                                                                     ggCTCGCCATTTGTCCTCAGCTCTGACCGAAAGATGAc                                                                                                                                    >  8:116726/1‑38 (MQ=255)
                                                                                                     ggCTCGCCATTTGTCCTCAGCTCTGACCGAAAGATGAc                                                                                                                                    <  7:120611/38‑1 (MQ=255)
                                                                                                     ggCTCGCCATTTGTCCTCAGCTCTGACAGAAAGATGAc                                                                                                                                    >  6:33948/1‑38 (MQ=255)
                                                                                                     ggCTCGCCATTTGTCCTCAGCTCTGACAGAAAGATGAc                                                                                                                                    <  5:33948/38‑1 (MQ=255)
                                                                                                                    cTCAGCTCTGACCGTAAGCTTACCCTTAGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTTGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGATCGAAATCAAACTGTCTAATAAGAAAATCaaaa                <  7:51152/139‑1 (MQ=255)
                                                                                                                    cTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGATCGAAATCAAACTGTCTAATAAGAAAATCaaaa                <  8:151099/139‑1 (MQ=255)
                                                                                                                        gCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTCAGCACGAGAAACGTGAAAACGATCGAAATCAAACTGTCTAATAAGAAAATCAAAACAGt            <  4:131899/139‑1 (MQ=255)
                                                                                                                            tGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGATCGAAATCAAACTGTCTAATAAGAAAATCAAAACAGTCa          <  5:136210/137‑1 (MQ=255)
                                                                                                                            tGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGATCGAAATCAAACTGTCTAATAAGAAAATCAAAACAGTCa          <  7:170360/137‑1 (MQ=255)
                                                                                                                               ccGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGATCGAAATCAAACTGTCTAATAAGAAAATCAAAACAGTCAAACTa     <  4:212772/139‑1 (MQ=255)
                                                                                                                                   aaGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGATCGAAATCAAACTGTCTAATAAGAAAATCAAAACAGTCAAACTAAAg  >  3:127062/1‑138 (MQ=255)
                                                                                                                                      |                                                                                                                                      
CCTATAACAAGTCTGTTGCAGGAGCTGTTGTTGATCCATTGCTTCCATGCATGCAGGTATCAGAGCTTGCATCCCTTAACAATAATACCTACCGGACACTCGGCTCGCCATTTGTCCTCAGCTCTGACCGAAAGCTGACCCTTCGTGTGGTGCAGGACGGCAACGAACACCCAATCATCGGCTTGGACAATGAAGCACTCAGCACGAGAAATGTGAAAACGATCGAAATCAAACTGTCTAATAAGAAAATCAAAACAGTCAAACTAAAG  >  BS168_reference_genome/3021726‑3021994

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