Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,787,719 T→G 6.8% S306R (AGC→CGC) ‡ PP_3350 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,787,7190TG 6.8% 224.3 / 6.2 73S306R (AGC→CGC) ‡PP_3350hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/31);  new base G (3/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 = 9.82e-01

TGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTG  >  NC_002947/3787572‑3787867
                                                                                                                                                   |                                                                                                                                                    
tggctTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTc                                                                                                                                                     <  7:186869/149‑1 (MQ=255)
    ttCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGa                                                                                                                                                   >  6:43826/1‑147 (MQ=255)
       gCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGAGCGACAgc                                                                                                                                               <  4:216957/148‑1 (MQ=255)
          aTAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGcc                                                                                                                                            >  2:169830/1‑148 (MQ=255)
            aGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCGCTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAg                                                                                                                                         >  3:200567/1‑149 (MQ=255)
             gACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAg                                                                                                                                         <  6:90937/148‑1 (MQ=255)
             gACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAg                                                                                                                                         >  5:90937/1‑148 (MQ=255)
              acacCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGgc                                                                                                                                       <  5:43826/149‑1 (MQ=255)
                acCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACa                                                                                                                                                 <  7:255131/137‑1 (MQ=255)
                acCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACa                                                                                                                                                 >  8:255131/1‑137 (MQ=255)
                  cATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGcgcg                                                                                                                                    >  5:167723/1‑148 (MQ=255)
                    ttGCTTGGGTAAAGAGGGTCGAAGGTGCCTACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGa                                                                                                                                  >  7:174891/1‑148 (MQ=255)
                    ttGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGAt                                                                                                                                 <  7:258843/149‑1 (MQ=255)
                      gCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGAGCGACAgcg                                                                                                                                              >  5:142016/1‑134 (MQ=255)
                      gCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGAGCGACAgcg                                                                                                                                              <  6:142016/134‑1 (MQ=255)
                            gTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGt                                                                                                                         >  7:61915/1‑149 (MQ=255)
                            gTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGAGCCCAGGCGCGGATCGACGCGt                                                                                                                         >  4:44502/1‑149 (MQ=255)
                                 gCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGcc                                                                                                                     >  4:336247/1‑148 (MQ=255)
                                   gggTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAg                                                                                                                  <  2:44100/149‑1 (MQ=255)
                                     gTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGt                                                                                                                >  5:201652/1‑149 (MQ=255)
                                     gTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGt                                                                                                                <  2:200860/149‑1 (MQ=255)
                                        gAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGc                                                                                                             >  3:157425/1‑149 (MQ=255)
                                               ccGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCta                                                                                                      >  2:118688/1‑148 (MQ=255)
                                                 gACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTAcc                                                                                                    <  4:157425/149‑1 (MQ=255)
                                                 gACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTAcc                                                                                                    <  1:118688/149‑1 (MQ=255)
                                                  aCTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTAccc                                                                                                   >  1:249672/1‑149 (MQ=255)
                                                        ccGTGGCCCAGGTCACTGTAGCCGCTTGCGGCATGGATGCGCATGGACAGGCGCGGTTGGCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGAGTATCGACGCGTTGCCCAGGTGGCCGGTCTGCCcgacga                                                                                              <  8:61915/148‑1 (MQ=255)
                                                        ccGTCGCCCAGGTCGCTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCcgacga                                                                                              <  4:200567/148‑1 (MQ=255)
                                                            cGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTc                                                                                         <  8:174891/149‑1 (MQ=255)
                                                               ccAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAg                                                                                      >  6:219484/1‑149 (MQ=255)
                                                                cAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGc                                                                                     <  6:167723/149‑1 (MQ=255)
                                                                 aGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCt                                                                                    >  6:243010/1‑149 (MQ=255)
                                                                  ggTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAg                                                                                             >  3:302154/1‑139 (MQ=255)
                                                                  ggTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAg                                                                                             <  4:302154/139‑1 (MQ=255)
                                                                   gTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTcc                                                                                        <  7:87724/143‑1 (MQ=255)
                                                                   gTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTcc                                                                                        >  8:87724/1‑143 (MQ=255)
                                                                    tCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCTACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCtcca                                                                                 <  2:249672/149‑1 (MQ=255)
                                                                     cACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGt                                                                                                                <  5:324803/117‑1 (MQ=255)
                                                                     cACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGt                                                                                                                >  6:324803/1‑117 (MQ=255)
                                                                     cACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCtccag                                                                                >  7:342569/1‑149 (MQ=255)
                                                                         gTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGAcgc                                                                                                                           >  8:227889/1‑102 (MQ=255)
                                                                          tCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCtccagtccag                                                                           >  6:189094/1‑149 (MQ=255)
                                                                              cgcTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGcc                                                                       <  3:336247/149‑1 (MQ=255)
                                                                                   gCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGAt                                                                   <  2:215655/148‑1 (MQ=255)
                                                                                   gCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGa                                                                    <  1:169830/147‑1 (MQ=255)
                                                                                   gCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATg                                                                  >  8:120826/1‑149 (MQ=255)
                                                                                      gCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATggcg                                                               >  5:233756/1‑149 (MQ=255)
                                                                                            aTGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGc                                                                        <  1:136146/134‑1 (MQ=255)
                                                                                            aTGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGc                                                                        >  2:136146/1‑134 (MQ=255)
                                                                                            aTGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGa                                                          >  6:114449/1‑148 (MQ=255)
                                                                                                 cATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGcc                                                    <  5:219484/149‑1 (MQ=255)
                                                                                                  aTGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCcg                                                   >  2:259684/1‑149 (MQ=255)
                                                                                                    ggCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCcgc                                                  <  8:342569/148‑1 (MQ=255)
                                                                                                     gCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGAACGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGt                                                >  6:346102/1‑149 (MQ=255)
                                                                                                      ccAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCtc                                                                                   <  6:31620/113‑1 (MQ=255)
                                                                                                      ccAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCtc                                                                                   >  5:31620/1‑113 (MQ=255)
                                                                                                          gcgcGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCg                                            >  6:36033/1‑148 (MQ=255)
                                                                                                            gcgGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGaaa                                         <  1:259684/149‑1 (MQ=255)
                                                                                                                   ccAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAg                                   <  6:201652/148‑1 (MQ=255)
                                                                                                                   ccAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAg                                   <  8:311003/148‑1 (MQ=255)
                                                                                                                    cAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAgcg                                 >  6:101710/1‑149 (MQ=255)
                                                                                                                     aGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAgcgc                                >  6:307824/1‑149 (MQ=255)
                                                                                                                        ggtgggtgAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCgctgc                             <  6:233756/149‑1 (MQ=255)
                                                                                                                          tgggtgAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTc                           <  5:101710/149‑1 (MQ=255)
                                                                                                                               gaaggtgaagCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCAc                      >  7:177522/1‑149 (MQ=255)
                                                                                                                               gaaggtgaagCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCAc                      >  8:193513/1‑149 (MQ=255)
                                                                                                                                 aggtgaagCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCAccc                    <  3:21215/149‑1 (MQ=255)
                                                                                                                                 aggtgaagCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCAccc                    >  2:599/1‑149 (MQ=255)
                                                                                                                                  ggtgaagCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCg                   >  7:320092/1‑149 (MQ=255)
                                                                                                                                   gtgaagCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGa                  >  3:181369/1‑149 (MQ=255)
                                                                                                                                        ggcGCTGTCGAGCGACAGCGGCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGCCGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGc             >  3:133333/3‑149 (MQ=255)
                                                                                                                                         gcGCTGTCGAGCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGcc            >  3:159172/2‑149 (MQ=255)
                                                                                                                                                 gCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAgg     <  5:243010/148‑1 (MQ=255)
                                                                                                                                                   tCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTg  >  1:196862/1‑149 (MQ=255)
                                                                                                                                                   tCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTACCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTg  >  5:225908/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGGGTG  >  NC_002947/3787572‑3787867

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