Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,901,484 A→G 5.5% R54R (AGA→AGG)  PP_3444 → glyoxalase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,901,4840AG 5.5% 217.1 / 2.1 73R54R (AGA→AGG) PP_3444glyoxalase family protein
Reads supporting (aligned to +/- strand):  ref base A (38/31);  new base G (2/2);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

CATCCCCGAAGGCCAGCACAGCATCACGCCCTACCTGGCCATCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGCGACAAGGTCTACGCCCAGGCCCTGGCCGCCGGGG  >  NC_002947/3901337‑3901632
                                                                                                                                                    |                                                                                                                                                    
cATCCCCGAAGGCCAGCACAGCATCACGCCCTACCTGGCCATCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGaa                                                                                                                                                     <  7:162128/149‑1 (MQ=255)
 aTCCCCGAAGGCCAGCACAGCATCACGCCCTACCTGGCCATCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGGAt                                                                                                                                                    >  2:90146/1‑149 (MQ=255)
 aTCCCCGAAGGCCAGCACAGCATCACGCCCTACCTGGCCATCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGGAt                                                                                                                                                    >  2:84476/1‑149 (MQ=255)
 aTCCCCGAAGGCCAGCACAGCATCACGCCCTACCTGGCCATCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAAt                                                                                                                                                    >  7:115856/1‑149 (MQ=255)
    cccGAAGGCCAGCACAGCATGACGCCCTACCTGGCCATCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCATGCCTTTGGTGCGGTGGAAATGTTGCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCg                                                                                                                                                  <  5:113643/148‑1 (MQ=255)
     ccGAAGGCCAGCACAGCATCACGCCCTACCTGGCCATCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCTAACTGAGAATCGGc                                                                                                                                                >  8:289267/1‑149 (MQ=255)
     ccGAAGGCCAGCACAGCATCACGCCCTACCTGGCCATCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGc                                                                                                                                                <  2:105879/149‑1 (MQ=255)
             cagcacagcaTCACGCCCTACCTGGCCATCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATtcgtc                                                                                                                                        >  2:394898/1‑149 (MQ=255)
               gcacagcaTCACGCCCTACCTGGCCATCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGc                                                                                                                                      <  7:289267/149‑1 (MQ=255)
                cacagcaTCACGCCCTACCTGGCCATCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCt                                                                                                                                     >  2:391094/1‑149 (MQ=255)
                   agcaTCACGCCCTACCTGGCCATCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGAt                                                                                                                                  >  5:1066/1‑149 (MQ=255)
                                  cTGGCCATCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTg                                                                                                                   >  8:156628/1‑149 (MQ=255)
                                     gCCATCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGTACTGAGGATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGa                                                                                                                <  1:84476/149‑1 (MQ=255)
                                     gCCATCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGGATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGa                                                                                                                <  1:90146/149‑1 (MQ=255)
                                        aTCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACAt                                                                                                             >  5:66834/1‑149 (MQ=255)
                                            acgacgCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATgg                                                                                                           >  1:4448/1‑147 (MQ=255)
                                            acgacgCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGaa                                                                                                         <  8:115856/149‑1 (MQ=255)
                                                cgccgcCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAggcg                                                                                                     <  1:391094/149‑1 (MQ=255)
                                                cgccgcCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAggcg                                                                                                     <  1:394898/149‑1 (MQ=255)
                                                 gccgccAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAggcgg                                                                                                    >  7:265144/1‑149 (MQ=255)
                                                  ccgccAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGAc                                                                                                               >  4:376265/1‑137 (MQ=255)
                                                  ccgccAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGAc                                                                                                               <  3:376265/137‑1 (MQ=255)
                                                  ccgccAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAggcggc                                                                                                   <  2:4448/149‑1 (MQ=255)
                                                  ccgccAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAggcggc                                                                                                   <  6:66834/149‑1 (MQ=255)
                                                    gccAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCt                                                                                                 >  5:327522/1‑149 (MQ=255)
                                                       aaGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCAc                                                                                              >  1:331198/1‑149 (MQ=255)
                                                            cATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCgc                                                                                           >  7:82091/1‑147 (MQ=255)
                                                              tCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCg                                                                                                                 >  5:228128/1‑123 (MQ=255)
                                                              tCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCg                                                                                                                 <  6:228128/123‑1 (MQ=255)
                                                                 aGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCgccagcca                                                                                     <  3:166412/148‑1 (MQ=255)
                                                                 aGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCgccagc                                                                                       >  4:166412/1‑146 (MQ=255)
                                                                    tCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAg                                                                                 >  4:343471/1‑149 (MQ=255)
                                                                     cTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGa                                                                                                                                   >  2:249198/1‑98 (MQ=255)
                                                                     cTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGa                                                                                                                                   <  1:249198/98‑1 (MQ=255)
                                                                      tACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAggcg                                                                                                     >  8:168589/1‑127 (MQ=255)
                                                                      tACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAggcg                                                                                                     <  7:168589/127‑1 (MQ=255)
                                                                        caagcaagCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTc                                                                              >  4:159569/1‑148 (MQ=255)
                                                                        caagcaagCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCa                                                                             <  3:343471/149‑1 (MQ=255)
                                                                           gcaagcCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACg                                                                          >  7:322611/1‑149 (MQ=255)
                                                                           gcaagcCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACg                                                                          >  5:225665/1‑149 (MQ=255)
                                                                             aagcCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGt                                                                        >  2:26481/1‑149 (MQ=255)
                                                                                    tGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGgt                                                                  <  6:1066/148‑1 (MQ=255)
                                                                                      gTGCGGTGGAAATGTTCCGCCTCGATGGCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGgtgt                                                                >  1:212828/1‑149 (MQ=255)
                                                                                       tGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTg                                                              <  2:331198/149‑1 (MQ=255)
                                                                                           gTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCt                                                          <  3:171909/149‑1 (MQ=255)
                                                                                                 aTGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACtt                                                    >  4:388606/1‑149 (MQ=255)
                                                                                                   gTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGccagccag                                                                                    >  2:412923/1‑115 (MQ=255)
                                                                                                   gTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGccagccag                                                                                    <  1:412923/115‑1 (MQ=255)
                                                                                                    ttCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGc                                                                      >  7:143159/1‑128 (MQ=255)
                                                                                                    ttCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGc                                                                      <  8:143159/128‑1 (MQ=255)
                                                                                                    tgcctcgtCGCT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGt                                                  <  1:141229/141‑1 (MQ=255)
                                                                                                           tCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAg                                          >  6:351217/1‑149 (MQ=255)
                                                                                                             gAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGa                                        <  8:265144/149‑1 (MQ=255)
                                                                                                                   cccAGGCGGTCGCGTTGGCCATGTCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGc                                    >  1:105896/1‑148 (MQ=255)
                                                                                                                    ccAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGa                                                                                                                                   <  8:312473/52‑1 (MQ=255)
                                                                                                                    ccAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGa                                                                                                                                   >  7:312473/1‑52 (MQ=255)
                                                                                                                    ccAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGa                                           <  1:242067/140‑1 (MQ=255)
                                                                                                                    ccAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGa                                           >  2:242067/1‑140 (MQ=255)
                                                                                                                      aGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGCGAc                                 >  8:64512/1‑148 (MQ=255)
                                                                                                                        gCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGCGACAAg                              >  1:27600/1‑149 (MQ=255)
                                                                                                                        gCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGCGACAAg                              <  7:64512/149‑1 (MQ=255)
                                                                                                                            tCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGCGACAAGGTc                           >  7:160707/1‑148 (MQ=255)
                                                                                                                              gcgTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGCGACAAGGTCTAc                        <  8:82091/149‑1 (MQ=255)
                                                                                                                                   ggCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGCGACAAGGTCTACGcc                     >  1:176177/1‑147 (MQ=255)
                                                                                                                                     ccATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGt                                                  >  2:392656/1‑116 (MQ=255)
                                                                                                                                     ccATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGt                                                  <  1:392656/116‑1 (MQ=255)
                                                                                                                                         gCCGAACAGAGAATCGGAGATTCGTCCATGATGCACGCCCACCCTTGCGACATCGAAGGCGGCCTCACCCCAAGCCAGAAGCTCAACGGTGCGGGTGTTTGCGTGCACTTGTATATTGAGGACTGCGACAAGGTCTACGCCCAGGCCCt             >  4:61212/1‑149 (MQ=255)
                                                                                                                                           cGAACTGAGAATCGGCGATTCGTCGCTGATGCTgg                                                                                                                             <  3:202787/35‑1 (MQ=255)
                                                                                                                                           cGAACTGAGAATCGGCGATTCGTCGCTGATGCTgg                                                                                                                             >  4:202787/1‑35 (MQ=255)
                                                                                                                                           cGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGCGACAAGGTCTACGCCCAGGCCCTgg           <  6:225665/149‑1 (MQ=255)
                                                                                                                                            gAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGCGACAAGGTCTACGCCCAGGCCCTGgc          <  6:327522/149‑1 (MQ=255)
                                                                                                                                             aaCTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGCGACAAGGTCTACGCCCAGGCCCTGgcc         <  3:388606/149‑1 (MQ=255)
                                                                                                                                                tGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGc                                    >  2:19692/1‑119 (MQ=255)
                                                                                                                                                tGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGc                                    <  1:19692/119‑1 (MQ=255)
                                                                                                                                                    aaTCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGCGACAAGGTCTACGCCCAGGCCCTGGCCGCCgggg  <  2:27600/149‑1 (MQ=255)
                                                                                                                                                    aaTCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTACGTTGAGGACTGCGACAAGGTCTACGCCCAGGCCCTGGCCGCCgggg  <  3:159569/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATCCCCGAAGGCCAGCACAGCATCACGCCCTACCTGGCCATCAACGACGCCGCCAAGGCCATCGAGTTCTACAAGCAAGCCTTTGGTGCGGTGGAAATGTTCCGCCTCGAT‑GCCCCAGGCGGTCGCGTTGGCCATGCCGAACTGAGAATCGGCGATTCGTCGCTGATGCTGGGCGACCCTTGCGACATGGAAGGCGGCCTCACCGCCAGCCAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGCGACAAGGTCTACGCCCAGGCCCTGGCCGCCGGGG  >  NC_002947/3901337‑3901632

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: ‑