Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,006,845 G→C 5.5% A344G (GCC→GGC)  PP_3535 ← gamma‑glutamyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,006,8450GC 5.5% 240.6 / 2.9 73A344G (GCC→GGC) PP_3535gamma‑glutamyltransferase
Reads supporting (aligned to +/- strand):  ref base G (34/35);  new base C (2/2);  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 = 9.94e-01

GGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGCCCAAGGTCTGTAGCACGGTCACCCCACCCGA  >  NC_002947/4006700‑4006993
                                                                                                                                                  |                                                                                                                                                     
ggCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGa                                                                                                                                                    <  1:137024/149‑1 (MQ=255)
 gcgcgcTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGAt                                                                                                                                                   >  3:146521/1‑149 (MQ=255)
  cgcgcTTCATGCTGTGCTCGCCGATCAGGGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCTATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATc                                                                                                                                                  <  3:338358/149‑1 (MQ=255)
  cgcgcTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCAGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGGCAGGTACAGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGAt                                                                                                                                                   >  4:338358/1‑148 (MQ=255)
      cTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGTCCGGGACGGGCTCGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGTGATCAGGt                                                                                                                                              <  3:101036/149‑1 (MQ=255)
          aTGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCAc                                                                                                                                          >  8:360499/1‑149 (MQ=255)
               gtgCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGcc                                                                                                                                      >  5:368996/1‑148 (MQ=255)
                tgCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCaa                                                                                                                                    <  2:317630/149‑1 (MQ=255)
                 gcgcGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAAc                                                                                                                                   <  8:216982/146‑1 (MQ=255)
                   tCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAAcgg                                                                                                                                 >  6:343958/1‑149 (MQ=255)
                     gCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAAcggcg                                                                                                                              >  3:318137/1‑149 (MQ=255)
                      ccGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAAcggcgg                                                                                                                             <  6:352175/149‑1 (MQ=255)
                             gcgcCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGCCGAGGCCGAGCGCCTGGCCTACGCCGATcgc                                                                                                                                                <  2:5970/124‑3 (MQ=25)
                             gcgcCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGCCGAGGCCGAGCGCCTGGCCTACGCCGATcgc                                                                                                                                                >  1:5970/1‑122 (MQ=25)
                                  gcgcgGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGCCGAGGCCGAGCGCCTGGCCTACGCCGATc                                                                                                                                                  <  4:280782/117‑1 (MQ=37)
                                  gcgcgGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGCCGAGGCCGAGCGCCTGGCCTACGCCGATc                                                                                                                                                  >  3:280782/1‑117 (MQ=37)
                                   cgcgGCTGGCCAGGTAAGCAGAATCGGTCCGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCGCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCgg                                                                                                                <  6:257854/149‑1 (MQ=255)
                                   cgcgGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCgg                                                                                                                <  2:50453/149‑1 (MQ=255)
                                      ggCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCa                                                                                                                                                 <  8:184951/114‑1 (MQ=255)
                                      ggCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCa                                                                                                                                                 >  7:184951/1‑114 (MQ=255)
                                       gCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCg                                                                                                            <  8:136935/149‑1 (MQ=255)
                                           gCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCggg                                                                                                          <  5:343958/147‑1 (MQ=255)
                                           gCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCggg                                                                                                          <  1:214038/147‑1 (MQ=255)
                                           gCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCggg                                                                                                          <  7:360499/147‑1 (MQ=255)
                                           gCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGa                                                                                                        <  6:368996/149‑1 (MQ=255)
                                            ccAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGAc                                                                                                       >  3:257796/1‑149 (MQ=255)
                                                 tAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGt                                                                                                  >  3:270900/1‑149 (MQ=255)
                                                 tAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGt                                                                                                  >  3:265922/1‑149 (MQ=255)
                                                          tCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATAGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCgg                                                                                         >  8:236114/1‑149 (MQ=255)
                                                             gTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCg                                                                                      >  1:86614/1‑149 (MQ=255)
                                                              tCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCg                                                                                      >  5:145967/1‑148 (MQ=255)
                                                                  ggCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAgcg                                                                                 <  4:146521/149‑1 (MQ=255)
                                                                    ccTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCt                                                                               >  2:78390/1‑149 (MQ=255)
                                                                        gACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGc                                                                             >  3:313126/1‑147 (MQ=255)
                                                                        gACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGc                                                                             >  7:117290/1‑147 (MQ=255)
                                                                           cGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCgctgc                                                                                      >  6:307810/1‑131 (MQ=255)
                                                                                 gggCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGtt                                                                  <  4:313126/149‑1 (MQ=255)
                                                                                  ggCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTg                                                                                                 <  6:42609/117‑1 (MQ=255)
                                                                                  ggCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTg                                                                                                 >  5:42609/1‑117 (MQ=255)
                                                                                    cACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCAGTCGCAGCGCTGCCAGGTCGCGTTgcg                                                               <  4:257796/149‑1 (MQ=255)
                                                                                           tCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCa                                                                           <  2:227541/130‑1 (MQ=255)
                                                                                           tCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTcgc                                                                     >  1:122931/1‑136 (MQ=255)
                                                                                           tCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTcgc                                                                     <  2:122931/136‑1 (MQ=255)
                                                                                           tCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGt                                                        >  5:256366/1‑149 (MQ=255)
                                                                                           tCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGt                                                        <  4:270900/149‑1 (MQ=255)
                                                                                           tCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGt                                                        <  4:265922/149‑1 (MQ=255)
                                                                                           tCGTTGTCGGCCAGGTACTGAGCGCGATCGACGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCa                                                                           >  1:227541/1‑130 (MQ=255)
                                                                                             gTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTgc                                                      >  6:235219/1‑149 (MQ=255)
                                                                                                   ggCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCa                                                <  6:196732/149‑1 (MQ=255)
                                                                                                        ggTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAgg                                                                                                                   <  6:36851/77‑1 (MQ=255)
                                                                                                        ggTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAgg                                                                                                                   >  5:36851/1‑77 (MQ=255)
                                                                                                         gTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGt                                          >  3:163959/1‑149 (MQ=255)
                                                                                                          tACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTc                                         <  1:78390/149‑1 (MQ=255)
                                                                                                               aGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCa                                    <  7:236114/149‑1 (MQ=255)
                                                                                                                gcgcGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAg                                   <  2:86614/149‑1 (MQ=255)
                                                                                                                 cgcgATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAgg                                  <  4:318137/149‑1 (MQ=255)
                                                                                                                   cgATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGcc                                >  1:244344/1‑149 (MQ=255)
                                                                                                                    gATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGcc                                >  6:235103/1‑148 (MQ=255)
                                                                                                                     aTCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCg                                        <  4:377487/139‑1 (MQ=255)
                                                                                                                     aTCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCg                                        >  3:377487/1‑139 (MQ=255)
                                                                                                                     aTCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCg                                        <  3:97006/139‑1 (MQ=255)
                                                                                                                     aTCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCg                                        >  4:97006/1‑139 (MQ=255)
                                                                                                                     aTCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGccc                               >  5:38558/1‑148 (MQ=255)
                                                                                                                       cGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGCCCAAg                            >  6:214890/1‑149 (MQ=255)
                                                                                                                         gCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGCCCAAGGt                          >  4:188921/1‑149 (MQ=255)
                                                                                                                            tAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGCCCAAGGTCTg                       <  5:214890/149‑1 (MQ=255)
                                                                                                                             aGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGCCCAAGGTCt                        >  6:385726/1‑147 (MQ=255)
                                                                                                                             aGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGCCCAAGGTCt                        <  5:385726/147‑1 (MQ=255)
                                                                                                                                 cAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGCCCAAGGTCt                        >  4:148269/1‑143 (MQ=255)
                                                                                                                                 cAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGCCCAAGGTCt                        <  3:148269/143‑1 (MQ=255)
                                                                                                                                   ggCGCGCGGCCGCCGCGAGCAGGTGCACGGCCAACGGCGGGGCTTCCACGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGCCCAAGGTCTGTAGCACg                <  6:256366/149‑1 (MQ=255)
                                                                                                                                      gcTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGAc                                                                                                       >  8:345025/1‑59 (MQ=255)
                                                                                                                                      gcTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGAc                                                                                                       <  7:345025/59‑1 (MQ=255)
                                                                                                                                        tcggcctcggcGATCAGGTGCACGGCCAACGGCGGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGCCCAAGGTCTGTAGCACg                >  1:253663/1‑145 (MQ=255)
                                                                                                                                        tcggcctcggcGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGCCCAAGGTCTGTAGCACg                <  2:253663/144‑1 (MQ=255)
                                                                                                                                        tcggcctcggcGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGCCCAAGGTCTGTAGCACGGTcac           >  1:217260/1‑149 (MQ=255)
                                                                                                                                                 ggcGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGCCCAAGGTCTGTAGCACGGTcaccccacc     <  6:38558/147‑1 (MQ=255)
                                                                                                                                                  gcGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGCCCAATGTCTGTAGCACGGTCACCCCACCCGa  <  6:100257/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGGCCAGGTAAGCAGAATCGGTCAGGGCCTTGACCGGGACGGGCACGTAGTCGTTGTCGGCCAGGTACTGAGCGCGATCGGCGTAGGCCAGGCGCTCGGCCTCGGCGATCAGGTGCACGGCCAACGGCGGGGCTTCCAGGCGGGCCGGGGACGAGGTGGCCAGCGGTCGCAGCGCTGCCAGGTCGCGTTGCGGCGCGGTGCGCTGCACGGCGTCGAGCAGGCCCAAGGTCTGTAGCACGGTCACCCCACCCGA  >  NC_002947/4006700‑4006993

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