Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,552,678 C→T 6.5% A112V (GCC→GTC)  PP_3137 → group 2 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,552,6780CT 6.5% 227.2 / 3.6 77A112V (GCC→GTC) PP_3137group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (37/35);  new base T (3/2);  total (40/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 = 1.10e-01

ACGCGCCCGCAATGAAGGTTTTGCCGCCTTGCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCTGTGGCGGCGATGCGCTGATGCGGGCAGACGCCTTTGCTGCC  >  NC_002947/3552541‑3552826
                                                                                                                                         |                                                                                                                                                    
aCGCGCCCGCAATGAAGGTTTTGCCGCCTTGCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAgtgg                                                                                                                                           >  8:89045/1‑149 (MQ=255)
      ccGCAATGAAGGTTTTGCCGCCTTGCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGtgtgtg                                                                                                                                     >  7:257927/1‑149 (MQ=255)
      ccGCAATGAAGGTTTTGCCGCCTTGCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGtgtgtg                                                                                                                                     >  4:275412/1‑149 (MQ=255)
              aaGGTTTTGCCGCCTTGCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGtgtgtg                                                                                                                                     >  8:62706/1‑141 (MQ=255)
              aaGGTTTTGCCGCCTTGCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGtgtgtg                                                                                                                                     <  7:62706/141‑1 (MQ=255)
                 gTTTTGCCGCCTTGCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGt                                                                                                                          <  3:184064/149‑1 (MQ=255)
                   tttGCCGCCTTGCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGt                                                                                                                          >  7:79435/1‑147 (MQ=255)
                      gccgccTTGCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAgc                                                                                                                      <  5:85926/148‑1 (MQ=255)
                       ccgccTTGCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAgcgc                                                                                                                    >  5:315673/1‑149 (MQ=255)
                       ccgccTTGCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAgcgc                                                                                                                    >  7:165177/1‑149 (MQ=255)
                        cgccTTGCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCt                                                                                                                   >  3:326720/1‑149 (MQ=255)
                           cTTGCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTcc                                                                                                                <  8:76145/149‑1 (MQ=255)
                            ttGCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGTCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTccc                                                                                                               >  2:214168/1‑149 (MQ=255)
                              gCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGATGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCg                                                                                                              <  8:165177/148‑1 (MQ=255)
                              gCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTccc                                                                                                               <  5:193971/147‑1 (MQ=255)
                               cAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAgtgg                                                                                                                                           >  7:271921/1‑118 (MQ=255)
                               cAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAgtgg                                                                                                                                           <  8:271921/118‑1 (MQ=255)
                               cAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCg                                                                                                              >  8:227192/1‑147 (MQ=255)
                               cAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCa                                                                                                            <  6:315673/149‑1 (MQ=255)
                                 gcgcCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGc                                                                                                          >  8:87464/1‑149 (MQ=255)
                                             aTCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTAc                                                                                               >  6:89365/1‑148 (MQ=255)
                                             aTCGATGCAGCTGGTACAGTTCGTAGATGGCGACTGAGAGGTGAATGACAGCTGGCTGGCCACTGCGCAGGCTTTTCTTGAGCAGCACACAGCAGTGGCAGTGGTGTGTGGTCGCCTGCGTGAGCGCTTACAGCAGCGTTAGGTTTAc                                                                                               >  7:261520/1‑148 (MQ=255)
                                                gATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAAc                                                                                            <  5:89365/148‑1 (MQ=255)
                                                gATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAAc                                                                                            <  4:271288/148‑1 (MQ=255)
                                                 aTGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACtt                                                                                          >  1:177073/1‑149 (MQ=255)
                                                  tGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACtt                                                                                          >  3:24481/1‑148 (MQ=255)
                                                   gCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGc                                                                                        >  4:275428/1‑149 (MQ=255)
                                                    cAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCt                                                                                       <  3:217517/149‑1 (MQ=255)
                                                      gCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCtgt                                                                                     <  4:326720/149‑1 (MQ=255)
                                                      gCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCtgt                                                                                     <  8:257927/149‑1 (MQ=255)
                                                       cTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCtgtg                                                                                    >  7:58442/1‑149 (MQ=255)
                                                          gTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTTAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGa                                                                                 <  8:244822/149‑1 (MQ=255)
                                                          gTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGc                                                                                   <  7:87464/147‑1 (MQ=255)
                                                          gTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGc                                                                                   <  2:200338/147‑1 (MQ=255)
                                                          gTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGa                                                                                 >  7:244822/1‑149 (MQ=255)
                                                          gTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGa                                                                                 >  1:316186/1‑149 (MQ=255)
                                                          gTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGa                                                                                 >  4:54810/1‑149 (MQ=255)
                                                           taaaGTTCCTCGATGGCGACTTCGAGGTGAATGCCAGCTGTCTTGCCACGTCTCAGGCTTTTCTTGAGCATCACCCCGCCGTGGCAGTGGTGTGTGGTCTCAGGAGTGAGATCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTTCGAc                                                                                <  8:261520/146‑1 (MQ=255)
                                                                ttCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTg                                                                             >  6:128112/1‑147 (MQ=255)
                                                                   gtcgATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTg                                                                        <  7:89045/149‑1 (MQ=255)
                                                                    tcgATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGGGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTgg                                                                       <  7:227192/149‑1 (MQ=255)
                                                                      gATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTggg                                                                      <  8:151037/148‑1 (MQ=255)
                                                                      gATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGa                                                                     <  3:275412/149‑1 (MQ=255)
                                                                      gATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCACCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGc                                                                                                             >  7:69523/1‑109 (MQ=255)
                                                                      gATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCACCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGc                                                                                                             <  8:69523/109‑1 (MQ=255)
                                                                      gATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAACACCCCGTCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGa                                                                     >  2:44516/1‑149 (MQ=255)
                                                                      gATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAACACCCCGTCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGa                                                                     >  2:39303/1‑149 (MQ=255)
                                                                         ggCGACTGCGAGGTGAATGGCAGCTGGCTGGCCACGGCGGAGGCTGGTCTTGAGCAGCGCACCGCCGGGGGAGGGGGGGGGGGGCGGGGGCGTGGGGGCGTCCGGCGGCGTGCGGTGGGCTAGTTGCTGGGCGTGCTGGCGTggggggg                                                                  >  6:350176/1‑145 (MQ=255)
                                                                                 cGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCcggt                                                            <  7:57018/147‑1 (MQ=255)
                                                                                    ggTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCcggtcggt                                                        >  2:75947/1‑148 (MQ=255)
                                                                                     gTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCt                                                                                       <  1:34361/116‑1 (MQ=255)
                                                                                     gTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCt                                                                                       >  2:34361/1‑116 (MQ=255)
                                                                                     gTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGa                                                      <  8:79435/149‑1 (MQ=255)
                                                                                     gTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGa                                                      <  3:275428/149‑1 (MQ=255)
                                                                                     gTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGCTGTACAACTTGCTGTGAGACCTGGAGTGGGATACCCcggtcggt                                                        >  5:63702/1‑147 (MQ=255)
                                                                                      tGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAg                                                     <  4:24481/149‑1 (MQ=255)
                                                                                      tGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAg                                                     <  3:54810/149‑1 (MQ=255)
                                                                                          tGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTg                                                                             <  5:181631/121‑1 (MQ=255)
                                                                                          tGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTg                                                                             >  6:181631/1‑121 (MQ=255)
                                                                                            ccAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAg                                               >  1:301416/1‑149 (MQ=255)
                                                                                                   gctggcCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGa                                                                     <  8:212067/120‑1 (MQ=255)
                                                                                                   gctggcCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGa                                                                     >  7:212067/1‑120 (MQ=255)
                                                                                                   gctggcCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCtgt                                         <  2:301416/148‑1 (MQ=255)
                                                                                                       gcCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCTGTggcgg                                    <  5:128112/149‑1 (MQ=255)
                                                                                                        ccACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCAGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACACCGGTCGGTGAGGCCAAGGCCTGTggcggc                                   >  1:243083/1‑149 (MQ=255)
                                                                                                             gcgcAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCTGTGGCGGCGATgc                              >  4:92120/1‑149 (MQ=255)
                                                                                                               gcAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCtgtg                                        <  8:32906/137‑1 (MQ=255)
                                                                                                               gcAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCtgtg                                        >  7:32906/1‑137 (MQ=255)
                                                                                                                 aGGCTTTTCTTGAGCAACACCCCGTCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCTGTGGCGGCGATGCGCTg                          <  1:39303/149‑1 (MQ=255)
                                                                                                                 aGGCTTTTCTTGAGCAACACCCCGTCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCTGTGGCGGCGATGCGCTg                          <  1:44516/149‑1 (MQ=255)
                                                                                                                     ttttCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCTGTGGCGGCGATGCGCTGATGc                      >  1:35528/1‑149 (MQ=255)
                                                                                                                      tttCTTGAGCAGCACCACGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGAAGCGTTCGGTGTACAACTAGATGTGCGCCCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCTGTGCCGGCGATGCGCTGATGCg                     >  5:186707/1‑149 (MQ=255)
                                                                                                                         cTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCTGTGGCGGCGATGCGCTGATGCggg                   >  1:139936/1‑148 (MQ=255)
                                                                                                                          gtGAGCAGCACGCCGCGGTGGCAGTGGGGTGTGGTCGCCGTCGTGAGCGCGTCCCGCAGCGTTCGGTGTATAACTTTCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCTGTGGCGGCGATGCGCTGATGCGGGc                  <  6:186707/147‑1 (MQ=255)
                                                                                                                              gcagcaCCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCTGTGGCGGCGATGCGCt                           <  8:40695/135‑1 (MQ=255)
                                                                                                                              gcagcaCCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCTGTGGCGGCGATGCGCt                           >  7:40695/1‑135 (MQ=255)
                                                                                                                                   aCCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCTGTGGCGGCGATGCGCTGATGCGGGCAGACGCCttt        >  4:248355/1‑149 (MQ=255)
                                                                                                                                      ccgccgTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGCCCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCTGTGGCGGCGATGCGCTGATGCGGGCAGACGCCTTtgct     >  4:27804/1‑149 (MQ=255)
                                                                                                                                       cgccgTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCTGTGGCGGCGATGCGCTGATGCGGGCAGACGCCTTtgctg    >  7:39445/1‑149 (MQ=255)
                                                                                                                                         ccgTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCTGTGGCGGCGATGCGCTGATGCGGGCAGACGCCTTTGCTGcc  >  4:55938/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
ACGCGCCCGCAATGAAGGTTTTGCCGCCTTGCAGCGCCTGTTGCCATCGATGCAGCTGGTACAGTTCGTCGATGGCGACTGCGAGGTGAATGCCAGCTGGCTGGCCACGGCGCAGGCTTTTCTTGAGCAGCACCCCGCCGTGGCAGTGGTGTGTGGTCGCCGGCGTGAGCGCTTCCCGCAGCGTTCGGTGTACAACTTGCTGTGCGACCTGGAGTGGGATACCCCGGTCGGTGAGGCCAAGGCCTGTGGCGGCGATGCGCTGATGCGGGCAGACGCCTTTGCTGCC  >  NC_002947/3552541‑3552826

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: