Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,787,718 C→A 6.8% S306I (AGC→ATC) ‡ PP_3350 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,787,7180CA 6.8% 232.1 / 5.9 73S306I (AGC→ATC) ‡PP_3350hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (37/31);  new base A (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 = 5.06e-01

GCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGG  >  NC_002947/3787570‑3787864
                                                                                                                                                    |                                                                                                                                                  
gctggctTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGc                                                                                                                                                    >  4:21215/1‑149 (MQ=255)
  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                                                                                                                                      >  5:90937/1‑148 (MQ=255)
               gACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAg                                                                                                                                      <  6:90937/148‑1 (MQ=255)
                acacCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGgc                                                                                                                                    <  5:43826/149‑1 (MQ=255)
                  acCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACa                                                                                                                                              >  8:255131/1‑137 (MQ=255)
                  acCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACa                                                                                                                                              <  7:255131/137‑1 (MQ=255)
                    cATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGcgcg                                                                                                                                 >  5:167723/1‑148 (MQ=255)
                      ttGCTTGGGTAAAGAGGGTCGAAGGTGCCTACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGa                                                                                                                               >  7:174891/1‑148 (MQ=255)
                      ttGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGAt                                                                                                                              <  7:258843/149‑1 (MQ=255)
                        gCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGAGCGACAgcg                                                                                                                                           <  6:142016/134‑1 (MQ=255)
                        gCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGAGCGACAgcg                                                                                                                                           >  5:142016/1‑134 (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                                                                                                 <  1:118688/149‑1 (MQ=255)
                                                   gACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTAcc                                                                                                 <  4:157425/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                                                                                                             >  6:324803/1‑117 (MQ=255)
                                                                       cACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGt                                                                                                             <  5:324803/117‑1 (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)
                                                                                                                                                    |                                                                                                                                                  
GCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGCCAGGCGCGGTTGCCAGGGGTGGGTGAAGGTGAAGCCGCTGTCGCTCGACAGCGCCCAGGCGCGGATCGACGCGTTGCCCAGGTGGCCGGTCTGCCCGACGAGGTTCCAGCTCCAGTCCAGTGCCTGATGGCGGCCGAACAGCCGCGTGCCGAAAGTGTAGCGCTGCTCGTCACCCGACTCGCCGGCCAGG  >  NC_002947/3787570‑3787864

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