Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,862,899 G→T 6.2% L6L (CTG→CTT cobD → Cobalamin biosynthesis protein CobD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,862,8990GT 6.2% 204.9 / 3.6 65L6L (CTG→CTTcobDCobalamin biosynthesis protein CobD
Reads supporting (aligned to +/- strand):  ref base G (37/24);  new base T (2/2);  total (39/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

GCCTTGGGCCGGTCGCCGAGTTCTATCCGGTGCCGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGCGGGCGGGCGCGGCTGGCGCAGCCACGG  >  ACB122/1862770‑1863033
                                                                                                                                 |                                                                                                                                      
gCCTTGGGCCGGTCGCCGAGTTCTATCCGGTGCCGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCg                                                                                                                                    <  1:284638/134‑1 (MQ=255)
gCCTAGGGCCGGTCGCCGAGGCCCATCCGGGGCCGATGTAGGTGCAGGAAAGCGGGGCCGAAGAGCGGCCCCTGAGTGAAATGCTTAACGAAAACCACTGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGt                                                                                                                                   <  1:344648/135‑1 (MQ=255)
 ccTTGGGCCGGTCGCCGAGTTCTATCCGGTGCCGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGt                                                                                                                                   >  1:144379/1‑134 (MQ=255)
 ccTTGGGCCGGTCGCCGAGTTCTATCCGGTGCCGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGt                                                                                                                                   <  2:144379/134‑1 (MQ=255)
    tGGGCCGGTCGCCGAGTTCTATCCGGTGCCGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCg                                                                                                                              >  7:112950/1‑136 (MQ=255)
       gCCGGTCGCCGAGTTCTATCCGGTGCCGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCGTGCTGCCCGGGGCCggg                                                                                                                            >  4:243553/1‑135 (MQ=255)
            tCGCCGAGTTCTATCCGGTGCCGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCggg                                                                                                                            <  7:256520/130‑1 (MQ=255)
            tCGCCGAGTTCTATCCGGTGCCGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCggg                                                                                                                            <  7:250169/130‑1 (MQ=255)
            tCGCCGAGTTCTATCCGGTGCCGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCggg                                                                                                                            >  8:256520/1‑130 (MQ=255)
            tCGCCGAGTTCTATCCGGTGCCGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCggg                                                                                                                            >  8:250169/1‑130 (MQ=255)
              gcggAGTTCTTTCCTGTGCGGATGTTTGTGCAGGAAAGCTGGGCGGAAGAGCGTCCCCTGTGTTAAATGCTTTACGAAAACCAATGGGTTGAGCGTCAATGAGCTTGGCCTTGCTGAGGGTGGCGGGGTTGGccc                                                                                                                     <  6:321814/132‑1 (MQ=255)
              gCCGAGTTCTATCCGGTGCCGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGccc                                                                                                                     >  6:214749/1‑135 (MQ=255)
                           cGGTGCCGATGTTGGTGCAGGAAAGCTGGGGCGAAGAGCGGCACCTGAGTGAAATGCTTTAGGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCt                                                                                                        <  3:74779/135‑1 (MQ=255)
                            ggTGCCGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTgg                                                                                                      <  7:330008/136‑1 (MQ=255)
                            ggTGCCGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTgg                                                                                                      <  7:397465/136‑1 (MQ=255)
                             gggCCGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTgg                                                                                                      <  5:55057/133‑1 (MQ=255)
                                 cGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGaa                                                                                                 >  7:289627/1‑136 (MQ=255)
                                 cGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGaa                                                                                                 >  7:113441/1‑136 (MQ=255)
                                 cGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGaa                                                                                                 >  7:343066/1‑136 (MQ=255)
                                 cGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGa                                                                                                  >  3:314711/1‑135 (MQ=255)
                                 cGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGa                                                                                                  >  5:12783/1‑135 (MQ=255)
                                           gCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTTACCGTGGCCGGGGTGGCCCTGGCCGCTTTGCTGGGCGAACCACAgcgc                                                                                        >  3:37819/1‑135 (MQ=255)
                                           gCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGCACCACAGCGcc                                                                                       >  7:315442/1‑136 (MQ=255)
                                           gCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAgcgc                                                                                        <  6:61315/135‑1 (MQ=255)
                                           gCAGGAAAGCTGGCCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTTACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAgcgc                                                                                        >  3:37815/1‑135 (MQ=255)
                                            cAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCg                                                                                      <  7:382736/136‑1 (MQ=255)
                                               gAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTc                                                                                    <  1:25143/135‑1 (MQ=255)
                                                 aaGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCAc                                                                                  >  4:85623/1‑135 (MQ=255)
                                                          cGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTtggtgg                                                                        >  7:377635/1‑136 (MQ=255)
                                                          cGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTtggtg                                                                         >  3:262748/1‑135 (MQ=255)
                                                                            tGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGcc                                                       >  1:52412/1‑135 (MQ=255)
                                                                               aaTGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCgg                                                     >  4:295939/1‑134 (MQ=255)
                                                                                 tGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCaa                                                  >  6:99045/1‑135 (MQ=255)
                                                                                       aCGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTg                                              <  5:161431/133‑1 (MQ=255)
                                                                                       aCGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGa                                            <  8:113441/135‑1 (MQ=255)
                                                                                       aCGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGa                                            >  5:107607/1‑135 (MQ=255)
                                                                                        cGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAg                                           <  8:341050/135‑1 (MQ=255)
                                                                                         gAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAgcg                                         >  8:244023/1‑136 (MQ=255)
                                                                                         gAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAgcg                                         >  7:80473/1‑136 (MQ=255)
                                                                                         gAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAgc                                          >  5:202439/1‑135 (MQ=255)
                                                                                          aaaaCCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAgcgc                                        <  7:292026/136‑1 (MQ=255)
                                                                                             aCCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGcgccg                                      <  6:202439/135‑1 (MQ=255)
                                                                                                aaTGGGGAGAGCGTCAATGAGCGAGGCCTTGCTGACCGTGGCCGCGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGcgccg                                      >  5:194651/1‑132 (MQ=255)
                                                                                                            gTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGcgggcggg                      >  7:196614/1‑136 (MQ=255)
                                                                                                              cAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGCgggcgggc                     >  6:206071/1‑135 (MQ=255)
                                                                                                               actGAGCGTGGCCTGGCTGCCCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGCGggcgggcg                    <  4:278806/133‑1 (MQ=255)
                                                                                                               aaTGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGCGGGCGGgcgc                   >  8:276418/1‑136 (MQ=255)
                                                                                                                    gCGTGGCGTTGCTGACGGTGGCGGGGGGGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGAAATATGGCCGGCAACCTGGCGCGCCGCCTGAAGGCGGGCGGGCGCGGCTg              <  8:377635/136‑1 (MQ=255)
                                                                                                                      gTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCCCCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGCGGGCGGGCGCGGCTGgc            >  8:245969/1‑136 (MQ=255)
                                                                                                                      gTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGCGGGCGGGCGCGGCTgg             >  3:261218/1‑135 (MQ=255)
                                                                                                                        ggCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAgcgc                                                                                        <  7:395069/58‑1 (MQ=255)
                                                                                                                        ggCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAgcgc                                                                                        >  8:395069/1‑58 (MQ=255)
                                                                                                                        ggCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGa                                            <  6:358642/102‑1 (MQ=255)
                                                                                                                        ggCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGa                                            >  5:358642/1‑102 (MQ=255)
                                                                                                                        ggCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGCGGGCGGGCGCGGCTGgcg           >  1:349654/1‑135 (MQ=255)
                                                                                                                        ggCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGCGGGCGGGCGCGGCTGgcg           >  1:240971/1‑135 (MQ=255)
                                                                                                                        ggCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGCGGGCGGGCGCGGCTGgc            >  4:18697/1‑134 (MQ=255)
                                                                                                                        ggCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGCGGGCGGGCGCGGCTGgc            >  1:223700/1‑134 (MQ=255)
                                                                                                                         gCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGCGGGCGGGCGCGGCTGgc            >  1:376402/1‑133 (MQ=255)
                                                                                                                          ccTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGCGGGCGGGCGCGGCTGGCGCa         >  6:52395/1‑135 (MQ=255)
                                                                                                                           cTTGCTTCCCGGGGCCGGGGGGGCCCTGGACGCTGTGCGGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATAGGGCCGGCAACCTGGAGCGCCGCCTGAATGCGGGCGGGCGCGGCTGCCGCAg        <  4:37819/135‑1 (MQ=255)
                                                                                                                           cTTGCTTACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGGGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGCGGGCGGGCGCGGCTGGCGCAg        <  4:37815/135‑1 (MQ=255)
                                                                                                                           cTTGCTGCCCGTGGCCGGGGGGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGGCACCCTTTGGTGGCGTGCGGCAATAGGGCCGGCAACCTGGAGCGCCGCCTGAATGCGGGCGGGCGCGGCTGGCGCAg        <  2:286532/135‑1 (MQ=255)
                                                                                                                           cTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGCGGGCGGGCGCGGCTGGCGCAg        <  3:295939/135‑1 (MQ=255)
                                                                                                                            atGCTGACCGTGGCCGGGGTGGCCCTGGACGCTATGCTGGGCGACCCACAGCGCCGACCCCCTTAGGTTGCGTTCGGCAATATGGCCGGGAACCTGGAGCGCCCCCGGAATGCGGGCGGGCCCGGCGGGCGCAGc       >  6:8420/2‑135 (MQ=255)
                                                                                                                                 gACCGTGGCCGGGGTGGCCCTGGACGCTTTTCTGGGCGAACCACAGCGCCCTCACCCTTTTGTGGCCTTCGGCAATATTGCCGGCAACCCGGCGCCCCCCCAGGAAGCGGGCGGGCCCCGCTGGCGCAggccccg  >  1:159052/1‑129 (MQ=255)
                                                                                                                                 gACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGCGGGCGGGCGCGGCTGGCGCAGCCACgg  >  4:240190/1‑135 (MQ=255)
                                                                                                                                 |                                                                                                                                      
GCCTTGGGCCGGTCGCCGAGTTCTATCCGGTGCCGATGTTGGTGCAGGAAAGCTGGGCCGAAGAGCGTCCCCTGAGTGAAATGCTTTACGAAAACCAATGGGGAGAGCGTCAATGAGCGTGGCCTTGCTGACCGTGGCCGGGGTGGCCCTGGACGCTTTGCTGGGCGAACCACAGCGCCGTCACCCTTTGGTGGCGTTCGGCAATATGGCCGGCAACCTGGAGCGCCGCCTGAATGCGGGCGGGCGCGGCTGGCGCAGCCACGG  >  ACB122/1862770‑1863033

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: