Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 446,384 C→G 5.2% A357G (GCC→GGC)  PP_0368 → acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947446,3840CG 5.2% 248.6 / 3.8 77A357G (GCC→GGC) PP_0368acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (35/38);  new base G (2/2);  total (37/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.63e-01

GAGGTGGCGTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGACCGAAGTCGGTTTC  >  NC_002947/446236‑446532
                                                                                                                                                     |                                                                                                                                                    
gAGGTGGCGTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGc                                                                                                                                                      <  8:255047/149‑1 (MQ=255)
gAGGTGGCGTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGc                                                                                                                                                      <  1:233069/149‑1 (MQ=255)
gAGGTGGCGTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGc                                                                                                                                                      >  2:109474/1‑149 (MQ=255)
    tGGCGTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGaa                                                                                                                                                  <  2:109819/149‑1 (MQ=255)
     ggCGTTCCAGGCTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTCACCCGCTCATCGTCCACCCGGACGTCCGTCGCATGCTGCTGACCATGAAGGCCTTCGcagcag                                                                                                                                                 >  6:376992/1‑144 (MQ=37)
      gCGTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAgg                                                                                                                                                >  6:29317/1‑149 (MQ=255)
      gCGTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAg                                                                                                                                                 <  1:266962/148‑1 (MQ=255)
       cGTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGaa                                                                                                                                                  >  2:158519/1‑146 (MQ=255)
       cGTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGaa                                                                                                                                                  <  1:158519/146‑1 (MQ=255)
        gTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGc                                                                                                                                                      >  5:79663/1‑141 (MQ=255)
        gTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGc                                                                                                                                                      <  6:79663/141‑1 (MQ=255)
            cAGGTTGGTCTGAAGTACGCCAGTGAGCTCCTGCAGATTCTTTCCCTGACCGTCCCGAAAGCACCGGACAAGGTCGCTGACCCGATAATCGTCCACCCGGACGTGCTGCTCATGCTGCTGACCATGAAGTCATTCGGCGAAGGccttca                                                                                                                                          <  2:260032/149‑6 (MQ=25)
            cAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAAcc                                                                                                                                           <  4:92759/148‑1 (MQ=255)
             aGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACcgc                                                                                                                                         >  4:143183/1‑149 (MQ=255)
                    ccTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATgg                                                                                                                                  >  1:377820/1‑149 (MQ=255)
                    ccTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATgg                                                                                                                                  >  8:132104/1‑149 (MQ=255)
                     cTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGt                                                                                                                                 >  1:56538/1‑149 (MQ=255)
                      tGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCGTGACGGGCCCGAAAGCACCGGACATGGCCGCTGACCCGTTCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCTTCCGCGCGATGGtt                                                                                                                                <  5:376992/149‑1 (MQ=37)
                                  gTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACtt                                                                                                                           <  5:335259/142‑1 (MQ=255)
                                  gTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACtt                                                                                                                           >  6:335259/1‑142 (MQ=255)
                                   tGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCa                                                                                                                   <  7:132104/149‑1 (MQ=255)
                                     aGCGCCTGCAGATGCGTTCCCTGACAGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAg                                                                                                                 <  5:29317/149‑1 (MQ=255)
                                      gcgcCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGc                                                                                                                <  5:292682/149‑1 (MQ=255)
                                         ccTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAgg                                                                                                             >  3:371777/1‑149 (MQ=255)
                                          cTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGt                                                                                                            >  5:188969/1‑149 (MQ=255)
                                          cTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGt                                                                                                            <  6:88770/149‑1 (MQ=255)
                                          cTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGt                                                                                                            >  3:365199/1‑149 (MQ=255)
                                             cAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGa                                                                                                         >  3:267166/1‑149 (MQ=255)
                                              aGATTCGTTCCCTGACCGGCCCGAAAGCACCGGACATTGCCGCTGACCCGATCATCGTCCACACGGACGTGCGTCGCATGCTGATGACCATGAAGGCCTTCGCCGAATGCAAAAGCTCGATGTTTTACTTCACCTCAAAGCTGGTGGAc                                                                                                        <  8:24036/149‑1 (MQ=25)
                                               gATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCAGGACGTGCGACGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCAACGCCAAGCAGGTGGAc                                                                                                        >  7:143939/1‑148 (MQ=255)
                                                   cGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGt                                                                                                   <  7:211737/149‑1 (MQ=255)
                                                    cttCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGACCACCCTGACGTGCGTCGCATGCTGCTGACCATAAAGGCCTTCGCCGAAGCCAACCCCGCGATGCTTTACATCACACCCAAGCAGGTCCACATAGt                                                                                                   >  2:70797/2‑148 (MQ=255)
                                                     ttCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCa                                                                                                 <  6:372046/149‑1 (MQ=255)
                                                       cccTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCaa                                                                                                <  1:75152/148‑1 (MQ=255)
                                                       cccTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCaa                                                                                                <  3:143183/148‑1 (MQ=255)
                                                        ccTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAg                                                                                               <  4:267166/148‑1 (MQ=255)
                                                         cTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTa                                                                                             <  8:162547/149‑1 (MQ=255)
                                                            aCCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAg                                                                                          <  6:188969/149‑1 (MQ=255)
                                                                  ccGAAAGCACCGGACAAGGCCGATGACCCGATCATCGTCCACCCGGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTg                                                                                                           >  1:79416/1‑127 (MQ=255)
                                                                  ccGAAAGCACCGGACAAGGCCGATGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTg                                                                                                           <  2:79416/126‑1 (MQ=255)
                                                                     aaaGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACga                                                                                 <  1:109474/149‑1 (MQ=255)
                                                                      aaGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACgaa                                                                                <  4:365199/149‑1 (MQ=255)
                                                                      aaGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACg                                                                                  >  8:90100/1‑147 (MQ=255)
                                                                         cACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACgaagaa                                                                             <  4:371777/149‑1 (MQ=255)
                                                                           ccGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAagaaga                                                                           <  1:70797/149‑1 (MQ=255)
                                                                           ccGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAagaaga                                                                           <  2:377820/149‑1 (MQ=255)
                                                                              gACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAgcg                                                                        <  2:56538/149‑1 (MQ=255)
                                                                                 aaGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAgcgc                                                                       >  3:368510/1‑147 (MQ=255)
                                                                                   ggCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCaag                                                                    >  3:60486/1‑148 (MQ=255)
                                                                                    gCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACgaag                                                                               >  4:169535/1‑136 (MQ=255)
                                                                                    gCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACgaag                                                                               <  3:169535/136‑1 (MQ=255)
                                                                                       gCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGc                                                               >  6:122352/1‑149 (MQ=255)
                                                                                         tGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTg                                                             >  8:121396/1‑149 (MQ=255)
                                                                                           aCCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGAc                                                           >  8:128325/1‑149 (MQ=255)
                                                                                             ccGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCaa                                                                                                >  4:323428/1‑110 (MQ=255)
                                                                                             ccGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCaa                                                                                                <  3:323428/110‑1 (MQ=255)
                                                                                                atcatcGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAgg                                                                >  7:35161/1‑139 (MQ=255)
                                                                                                atcatcGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAgg                                                                <  8:35161/139‑1 (MQ=255)
                                                                                                  catcGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCctg                                                     <  5:234387/148‑1 (MQ=255)
                                                                                                      gTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGt                                                                                              <  8:94157/103‑1 (MQ=255)
                                                                                                      gTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGt                                                                                              >  7:94157/1‑103 (MQ=255)
                                                                                                           cccGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGTTGGACATCGTCAAGTACAGCCAGGACGATGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTcc                                           >  6:14510/1‑149 (MQ=255)
                                                                                                           cccGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTcc                                           >  5:89559/1‑149 (MQ=255)
                                                                                                           cccGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTcc                                           >  8:231591/1‑149 (MQ=255)
                                                                                                           cccGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTcc                                           <  7:128325/149‑1 (MQ=255)
                                                                                                            ccGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCt                                          <  1:138514/149‑1 (MQ=255)
                                                                                                                       gTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATc                                >  6:47841/1‑149 (MQ=255)
                                                                                                                           cATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCa                            <  4:368510/149‑1 (MQ=255)
                                                                                                                                tgctgACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCa                       >  1:231749/1‑149 (MQ=255)
                                                                                                                                  ctgACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCt                                          >  8:99987/1‑128 (MQ=255)
                                                                                                                                    gACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGcc                             <  2:262370/139‑1 (MQ=255)
                                                                                                                                    gACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGcc                             >  1:262370/1‑139 (MQ=255)
                                                                                                                                    gACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCa                   >  7:179384/1‑149 (MQ=255)
                                                                                                                                       cATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACgaag                                                                               <  4:134236/86‑1 (MQ=255)
                                                                                                                                       cATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACgaag                                                                               >  3:134236/1‑86 (MQ=255)
                                                                                                                                          ggttgCCTTCGGCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCAt                      >  4:160055/5‑140 (MQ=255)
                                                                                                                                          ggttgCCTTCGGCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCAt                      <  3:160055/136‑1 (MQ=255)
                                                                                                                                           gttgCCTTCGGCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGACCg             >  6:299939/4‑148 (MQ=255)
                                                                                                                                            aGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGACCGaa           <  7:231591/149‑1 (MQ=255)
                                                                                                                                                 ttCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGACCGAAGTc        <  4:366363/147‑1 (MQ=255)
                                                                                                                                                 ttCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGACCGAAGTc        >  3:366363/1‑147 (MQ=255)
                                                                                                                                                  tCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGACCGAAGTCGGt     >  7:93756/1‑149 (MQ=255)
                                                                                                                                                   cGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGACCGAAGTCGGt     <  5:47841/148‑1 (MQ=255)
                                                                                                                                                     ccGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGACCGAAGTCGGTTTc  <  4:60486/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAGGTGGCGTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGACAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGACCGAAGTCGGTTTC  >  NC_002947/446236‑446532

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: