Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,619,471 A→G 7.0% C344R (TGT→CGT)  PP_4938 ← glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,619,4710AG 7.0% 251.4 / 6.5 87C344R (TGT→CGT) PP_4938glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (40/40);  new base G (2/4);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCCTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAGCCCCAG  >  NC_002947/5619325‑5619615
                                                                                                                                                  |                                                                                                                                                 
gccgccTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAgc                                                                                                                                                 >  7:218269/1‑149 (MQ=255)
gccgccTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAg                                                                                                                                                  <  5:60160/148‑1 (MQ=255)
  cgccTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGgcgc                                                                                                                                               >  6:257272/1‑149 (MQ=255)
     cTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGGCGCTca                                                                                                                                            <  4:266467/149‑2 (MQ=255)
        cGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGGCGCTcagt                                                                                                                                          <  5:257272/148‑4 (MQ=255)
        cGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGGCGCTcag                                                                                                                                           <  1:210452/147‑3 (MQ=255)
        cGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGa                                                                                                                                           <  5:311759/147‑1 (MQ=255)
            gCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGcc                                                                                                                                      >  6:128710/1‑148 (MQ=255)
            gCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCt                                                                                                                                     >  8:53250/1‑149 (MQ=255)
                      tGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGGCGCTcagt                                                                                                                                          >  5:29217/1‑131 (MQ=255)
                      tGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGGCGCTcagt                                                                                                                                          <  6:29217/134‑4 (MQ=255)
                       gACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCTCCGTACAGCGCTCGATCGCCTCAGCCCAGtcg                                                                                                                         >  8:308952/1‑149 (MQ=255)
                       gACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATc                                                                                                                                         >  4:121201/1‑134 (MQ=255)
                       gACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATc                                                                                                                                         <  3:121201/134‑1 (MQ=255)
                       gACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGgcg                                                                                                                         >  5:78462/1‑149 (MQ=255)
                           cGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCttt                                                                                                                     <  4:124196/149‑1 (MQ=255)
                             cTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTtcg                                                                                                                   >  6:122486/1‑149 (MQ=255)
                                aaGGTCAGCTGCGCAACAGGCGCTGCGAGCCGCTGCAGAACGCCGGCCAGGCCTGCTCCGGTTCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATTCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGAGCTttcttc                                                                                                                 <  2:101025/148‑1 (MQ=255)
                                 aGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCgg                                                                                                               <  5:48532/149‑1 (MQ=255)
                                  tgTCAGCTGCGCATCAGGCTCTGCGAGGCGCTGCTGAACGGCTGCCAGGCCTGCTCCGGTGCCAGTGCCTTGCGCTGCTTAGCGGCGATTGTCGTTGGACTGAGCGCCGTAGAGCTGTCGATCGCCTCAGGCCAGGCGATTTCGTCgg                                                                                                               <  4:270599/147‑1 (MQ=255)
                                       gCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCg                                                                                                         >  4:143438/1‑149 (MQ=255)
                                        cTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCgg                                                                                                        <  5:122486/149‑1 (MQ=255)
                                           cgcAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGc                                                                                                       <  5:110387/147‑1 (MQ=255)
                                           cgcAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAg                                                                                                     <  4:265378/149‑1 (MQ=255)
                                           cgcAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAg                                                                                                     <  4:260398/149‑1 (MQ=255)
                                           cgcAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAg                                                                                                     <  7:53250/149‑1 (MQ=255)
                                             cAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAgg                                                                                                    >  8:157247/1‑148 (MQ=255)
                                              aaCAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTa                                                                                                  <  2:79512/149‑1 (MQ=255)
                                                   gCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCAc                                                                                             <  8:218269/149‑1 (MQ=255)
                                                   gCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCAc                                                                                             <  5:128710/149‑1 (MQ=255)
                                                   gCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCAc                                                                                             <  4:69976/149‑1 (MQ=255)
                                                    ctctGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCAcc                                                                                            >  2:52049/1‑149 (MQ=255)
                                                     tctGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCCGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCAccc                                                                                           >  2:171709/1‑149 (MQ=255)
                                                     tctGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCAccc                                                                                           >  6:364601/1‑149 (MQ=255)
                                                       tGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCgt                                                                                         >  6:13517/1‑149 (MQ=255)
                                                         cGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCTACCTGCCGGTTGccgacgaaa                                                                                        <  6:78462/148‑8 (MQ=255)
                                                         cGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCgt                                                                                         <  1:38049/147‑1 (MQ=255)
                                                          gAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTc                                                                                      <  1:52049/149‑1 (MQ=255)
                                                            gCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTcgc                                                                                    <  8:187785/149‑1 (MQ=255)
                                                            gCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGAGAcgc                                                                                    >  4:61406/1‑149 (MQ=255)
                                                                    aGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGc                                                                            >  6:1373/1‑149 (MQ=255)
                                                                        cGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTcgc                                                                        <  6:169339/149‑1 (MQ=255)
                                                                          ccTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTcgcg                                                                       >  2:215100/1‑148 (MQ=255)
                                                                              ccAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCAGCCGTGTCGCGCAACTGCTCGCGaaa                                                                    <  3:61406/147‑1 (MQ=255)
                                                                              ccAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAaca                                                                  >  1:352235/1‑149 (MQ=255)
                                                                                  gCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTc                                                              >  6:206058/1‑149 (MQ=255)
                                                                                    cTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTc                                                              >  2:147826/1‑147 (MQ=255)
                                                                                    cTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCa                                                            <  5:364601/149‑1 (MQ=255)
                                                                                     tGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAg                                                           >  2:79607/1‑149 (MQ=255)
                                                                                        tCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTc                                                        >  7:2077/1‑149 (MQ=255)
                                                                                          cGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTc                                                        <  1:147826/147‑1 (MQ=255)
                                                                                          cGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGc                                                      >  2:61077/1‑149 (MQ=255)
                                                                                          cGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACGGCTCGCGAAACACTTCCAGGTCGc                                                      >  4:279445/1‑149 (MQ=255)
                                                                                           ggTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCg                                                       >  7:223238/1‑147 (MQ=255)
                                                                                           ggTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCg                                                       <  8:223238/147‑1 (MQ=255)
                                                                                            gTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTg                                                    <  2:352235/149‑1 (MQ=255)
                                                                                                 aGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCATGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACtt                                                               <  7:286929/133‑1 (MQ=255)
                                                                                                 aGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCacac                                               >  6:223300/1‑149 (MQ=255)
                                                                                                 aGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTAGCGCAACTGCTCGCGAAACACtt                                                               >  8:286929/1‑133 (MQ=255)
                                                                                                     ccTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACaccccc                                           >  5:363743/1‑149 (MQ=255)
                                                                                                         gcgcTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGt                                        <  5:206058/148‑1 (MQ=255)
                                                                                                         gcgcTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGt                                        <  1:61077/148‑1 (MQ=255)
                                                                                                                ttAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCCAGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATAGGc                                >  5:89645/1‑149 (MQ=255)
                                                                                                                    cggcgATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATc                            <  7:157247/149‑1 (MQ=255)
                                                                                                                      gcgATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTg                                                    >  2:192324/1‑123 (MQ=255)
                                                                                                                      gcgATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCAc                          >  5:324777/1‑149 (MQ=255)
                                                                                                                      gcgATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCAc                          >  1:237577/1‑149 (MQ=255)
                                                                                                                        gATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCg                        >  6:380191/1‑149 (MQ=255)
                                                                                                                          ttGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGcc                      <  8:2077/149‑1 (MQ=255)
                                                                                                                             tCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCctct                   <  5:223300/149‑1 (MQ=255)
                                                                                                                             tCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCctct                   <  1:215100/149‑1 (MQ=255)
                                                                                                                               gATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCtt                  <  1:171709/148‑1 (MQ=255)
                                                                                                                               gATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCCGTACATCGGCCATCACCGCCTCTTg                 >  1:168867/1‑149 (MQ=255)
                                                                                                                                aTGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTg                 >  2:185393/1‑148 (MQ=255)
                                                                                                                                    aCTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCtt                  >  6:21873/1‑143 (MQ=255)
                                                                                                                                    aCTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCtt                  <  5:21873/143‑1 (MQ=255)
                                                                                                                                     cTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGTCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGcaccac           <  5:380191/149‑1 (MQ=255)
                                                                                                                                      tGAGCGCCGTACAGCGCTCGATCGCCTCAGCGCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCCGCAcaccac                                            >  1:40626/1‑116 (MQ=255)
                                                                                                                                      tGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGcaccac           >  7:97965/1‑148 (MQ=255)
                                                                                                                                      tGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCCGCAcaccac                                            <  2:40626/115‑1 (MQ=255)
                                                                                                                                       gAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGcaccac           >  2:64598/1‑147 (MQ=255)
                                                                                                                                          cgcCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCCTCGGCAACCGGCAGGTAGCCACCCGTGTCTCGCAACTGCTCGCCAAACACTACCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAGcc      >  2:283322/1‑149 (MQ=255)
                                                                                                                                             cGTACCGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCGCAGCACCGGTACATCGGCCAGCCCCGCCTCTTGGACCACCCGccc     >  4:130825/1‑147 (MQ=255)
                                                                                                                                             cGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAGccc     >  3:65631/1‑147 (MQ=255)
                                                                                                                                             cGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAGCCCCa   <  7:308952/149‑1 (MQ=255)
                                                                                                                                              gTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAGCCCCAg  <  7:27495/149‑1 (MQ=255)
                                                                                                                                              gTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAGCCCCAg  <  3:143438/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCCGCCTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAGCCCCAG  >  NC_002947/5619325‑5619615

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: