Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 854,663 G→A 6.9% A88T (GCC→ACC)  moeB → molybdopterin synthase adenylyltransferase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947854,6630GA 6.9% 212.0 / 4.3 72A88T (GCC→ACC) moeBmolybdopterin synthase adenylyltransferase subunit
Reads supporting (aligned to +/- strand):  ref base G (37/30);  new base A (2/3);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

TGGCCTGGGCTCGCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTACCCGCGAGGCGGTC  >  NC_002947/854515‑854809
                                                                                                                                                    |                                                                                                                                                  
tGGCCTGGGCTCGCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCgg                                                                                                                                                    <  7:70025/149‑1 (MQ=255)
 ggCCTGGGCTCGCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGc                                                                                                                                                   >  3:15367/1‑149 (MQ=255)
  gCCTGGGCTCGCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGcc                                                                                                                                                  <  8:315655/149‑1 (MQ=255)
            gCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGCCGAGTCGACCTTg                                                                                                                                               <  1:291804/142‑1 (MQ=255)
            gCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGCCGAGTCGACCTTg                                                                                                                                               >  2:291804/1‑142 (MQ=255)
            gCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGCCGAGTCGACCTTGctcat                                                                                                                                          <  2:238417/147‑5 (MQ=255)
             ccAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCtt                                                                                                                                       >  2:430203/1‑149 (MQ=255)
                gTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCa                                                                                                                                    <  6:186795/149‑1 (MQ=255)
                   gCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCCGCGCTTGCAgg                                                                                                                                  >  4:266362/1‑148 (MQ=255)
                   gCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAgg                                                                                                                                  <  7:21512/148‑1 (MQ=255)
                          aCCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAAc                                                                                                                          >  8:172124/1‑149 (MQ=255)
                            cTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAAccc                                                                                                                        >  2:6963/1‑149 (MQ=255)
                               gccgccGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCgg                                                                                                                      <  4:15367/148‑1 (MQ=255)
                               gccgccGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCgg                                                                                                                      <  6:208316/148‑1 (MQ=255)
                               gccgccGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGa                                                                                                                     <  8:404599/149‑1 (MQ=255)
                                ccgccgCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGaa                                                                                                                    <  8:327684/149‑1 (MQ=255)
                                ccgccgCTGGCGTGGGTCAGCTGCACCTGGCCGACTTCGACCCCGTCGACCTGACCAACCTGCATCGCCCGGTCATCCACGACAGGGCCAGCGGTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGTaa                                                                                                                    <  4:312477/149‑1 (MQ=255)
                                    cgcTTGGGTGGTTGAGCTGCACCTGGGGGACTTCGTCACAGTCGACCTGACCAAGCTGCAACGCCAGGTCATCCACGACAGTGCAAGCGTTGGAATTAGCTAGGTCGACTAGGCCTTGCAGCGCTTGCAGGCAATCAACCCTGAAGTCa                                                                                                                <  8:195145/149‑1 (MQ=255)
                                         gCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATc                                                                                                                             <  6:192825/131‑1 (MQ=255)
                                         gCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATc                                                                                                                             >  5:192825/1‑131 (MQ=255)
                                           gTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGt                                                                                                         >  7:52193/1‑149 (MQ=255)
                                             gggTGAGCTGCACCTGGCCGACTTCGACACCGGCGACCTGACCACCCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTg                                                                                                       <  8:165099/149‑1 (MQ=255)
                                               gTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCt                                                                                                     >  7:62336/1‑149 (MQ=255)
                                                       cACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCa                                                                                             <  2:248803/149‑1 (MQ=255)
                                                       cACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCa                                                                                             <  5:440868/149‑1 (MQ=255)
                                                        aCCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAg                                                                                            >  2:234325/1‑149 (MQ=255)
                                                          cTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAgg                                                                                           >  7:326104/1‑148 (MQ=255)
                                                             gCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACt                                                                                       >  4:223255/1‑149 (MQ=255)
                                                               cGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTg                                                                                      >  1:228074/1‑148 (MQ=255)
                                                                gACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGa                                                                                    <  1:234325/149‑1 (MQ=255)
                                                                   ttCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGa                                                                                 <  7:59403/149‑1 (MQ=255)
                                                                   ttCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTCCCAGCGTTGGCAAGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAACTCAGCCTGGTTGCTCACCGCCAGGCACTGGAc                                                                                   >  8:135437/1‑147 (MQ=255)
                                                                     cGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAg                                                                                <  8:326104/148‑1 (MQ=255)
                                                                      gACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGa                                                                              >  7:80348/1‑149 (MQ=255)
                                                                        cacCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAgg                                                                                           <  5:211198/134‑1 (MQ=255)
                                                                        cacCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAgg                                                                                           >  6:211198/1‑134 (MQ=255)
                                                                                 ccTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCtggtgg                                                                   >  5:177690/1‑149 (MQ=255)
                                                                                       cAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGAc                                                                             <  1:173712/133‑1 (MQ=255)
                                                                                       cAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCtggtg                                                                    >  4:99255/1‑142 (MQ=255)
                                                                                       cAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCtggtg                                                                    <  3:99255/142‑1 (MQ=255)
                                                                                       cAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCgg                                                             >  5:225555/1‑149 (MQ=255)
                                                                                       cAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGAATGAGCAAGGTCGACTCGGCCTTGAAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGCTGGCTGCgg                                                             >  7:13438/1‑149 (MQ=255)
                                                                                        aaCCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGt                                                            >  8:154641/1‑149 (MQ=255)
                                                                                           cTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTggc                                                         >  5:4944/1‑149 (MQ=255)
                                                                                           cTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTggc                                                         >  4:180991/1‑149 (MQ=255)
                                                                                                 cGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGc                                                                                               >  2:401543/1‑105 (MQ=255)
                                                                                                 cGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGc                                                                                               <  1:401543/105‑1 (MQ=255)
                                                                                                  tccAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGTCAATCAACCCGGAAGTCAGCATGGTTGCTCACCGCAAGGCACTGTACGAGGACTCGATGGTGTCTGCGTTGGCGGCGGTc                                                  >  2:127557/2‑149 (MQ=255)
                                                                                                  gCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGACCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTc                                                  >  6:233317/1‑149 (MQ=255)
                                                                                                    cAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGa                                                >  5:273975/1‑149 (MQ=255)
                                                                                                     aGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGAc                                               >  2:406573/1‑149 (MQ=255)
                                                                                                         cATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTgg                                           >  7:62458/1‑149 (MQ=255)
                                                                                                          aTCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTgg                                           <  6:273975/148‑1 (MQ=255)
                                                                                                          aTCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGt                                          <  1:430203/149‑1 (MQ=255)
                                                                                                           tCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTg                                         <  2:361598/149‑1 (MQ=255)
                                                                                                           tCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTg                                         >  1:361598/1‑149 (MQ=255)
                                                                                                            ccACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTTGCGGCGGTCGACCTGGTGc                                        >  8:78573/1‑149 (MQ=255)
                                                                                                            ccACGACAGTGCCAGCGTTGGCATGAGCAAGCTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGAGc                                        >  4:96420/1‑149 (MQ=255)
                                                                                                              aCGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAg                                                                                <  6:243491/107‑1 (MQ=255)
                                                                                                              aCGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAg                                                                                >  5:243491/1‑107 (MQ=255)
                                                                                                              aCGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTg                                      <  8:80348/149‑1 (MQ=255)
                                                                                                                gACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGa                                    <  7:154641/149‑1 (MQ=255)
                                                                                                                      gCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGCCCTGGTGCTGGACTGCTc                              >  2:67933/1‑149 (MQ=255)
                                                                                                                         aGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGTGGTCGACCTGGTGCTGGACTGCTCCGa                           >  4:253397/1‑149 (MQ=255)
                                                                                                                          gCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGATTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGa                           >  5:79511/1‑148 (MQ=255)
                                                                                                                             ttGGCATGAGCAAGGTCGACTCGACCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAAt                       <  5:233317/149‑1 (MQ=255)
                                                                                                                                 cATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTtct                   <  3:223255/149‑1 (MQ=255)
                                                                                                                                 cATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTtct                   <  1:406573/149‑1 (MQ=255)
                                                                                                                                       cAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTAccc             >  4:117808/1‑149 (MQ=255)
                                                                                                                                       cAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCAGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCAGGAGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTAccc             >  6:73522/1‑149 (MQ=255)
                                                                                                                                            tCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTACCCGCGAg        <  8:62458/149‑1 (MQ=255)
                                                                                                                                               aCTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTACCCGCGAGGCg     <  8:227719/149‑1 (MQ=255)
                                                                                                                                               aCTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTACCCGCGAGGCg     >  6:428982/1‑149 (MQ=255)
                                                                                                                                                   ggCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTACCCGCGAGGCGGTc  <  5:428982/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGGCCTGGGCTCGCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTACCCGCGAGGCGGTC  >  NC_002947/854515‑854809

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