Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,120,934 T→C 5.4% V125V (GTT→GTC)  PP_0982 → YjgP/YjgQ family permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,120,9340TC 5.4% 221.8 / 2.2 75V125V (GTT→GTC) PP_0982YjgP/YjgQ family permease
Reads supporting (aligned to +/- strand):  ref base T (31/39);  new base C (2/2);  total (34/41)
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.73e-01

CGGGCGCCTGTACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGACCGCATCAACCTGGGCGGGG  >  NC_002947/1120787‑1121079
                                                                                                                                                     |                                                                                                                                                 
cGGGCGCCTGTACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTg                                                                                                                                                  <  5:147385/149‑1 (MQ=255)
   gcgcCTGTACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGcc                                                                                                                                                >  4:77297/1‑148 (MQ=255)
    cgcCTGTACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGcc                                                                                                                                                >  5:125089/1‑147 (MQ=255)
     gcCTGTACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAg                                                                                                                                             <  3:77297/149‑1 (MQ=255)
      ccTGTACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCA‑ACGCCCAgg                                                                                                                                            <  2:144126/149‑1 (MQ=255)
        tGTACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTg                                                                                                                                          >  4:37288/1‑149 (MQ=255)
         gTACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTg                                                                                                                                          >  5:3088/1‑148 (MQ=255)
         gTACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTg                                                                                                                                          >  2:223001/1‑148 (MQ=255)
          tACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCA‑ACGCCCAggggg                                                                                                                                         <  6:272280/148‑4 (MQ=255)
            ccTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTgcag                                                                                                                                       >  8:162995/1‑148 (MQ=255)
                gAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTgcag                                                                                                                                       >  4:220967/1‑144 (MQ=255)
                        gaTGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATcagcc                                                                                                                          <  6:221459/149‑1 (MQ=255)
                        gaTGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATcagc                                                                                                                           >  1:282118/1‑148 (MQ=255)
                            aCCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCcagca                                                                                                                      <  1:172416/149‑1 (MQ=255)
                                     tCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGcc                                                                                                               <  2:282118/147‑1 (MQ=255)
                                     tCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGcc                                                                                                               <  6:125089/147‑1 (MQ=255)
                                     tCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGcc                                                                                                               <  3:220967/147‑1 (MQ=255)
                                     tCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCt                                                                                                             >  6:201138/1‑149 (MQ=255)
                                       ggCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGa                                                                                                           <  8:119528/149‑1 (MQ=255)
                                        gCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGAc                                                                                                          >  3:281302/1‑149 (MQ=255)
                                         ccACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGCTCCATGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGAc                                                                                                          <  1:102375/149‑1 (MQ=255)
                                          cACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTACTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGAcc                                                                                                         >  5:116531/1‑148 (MQ=255)
                                           aCTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGAcc                                                                                                         >  5:234363/1‑147 (MQ=255)
                                             tGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCCCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGATGAGCCAGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCAACAGCCAGAAGGATGCCCTGACCGAGt                                                                                                     >  3:296018/1‑149 (MQ=255)
                                                 aTGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGa                                                                                                 >  3:235927/1‑149 (MQ=255)
                                                  tGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGac                                                                                                >  3:149900/1‑149 (MQ=255)
                                                    ttccAGCATCGGCTTCTGGGGATGACCATG‑GCGGCGGCTGCCCTGGGTGCCCTGCTGGTGGCCTGGCTGTGCCTGAGCCTGGCCCCCCTTGGCGTTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGac                                                                                                <  1:76370/146‑1 (MQ=255)
                                                      ccAGCAGCGCCTGCTGTGCATGACCATG‑GCGCCGGCTGCCCTGGGTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCt                                                                                            <  5:97452/149‑1 (MQ=255)
                                                       ctgcagcGTCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGTCTGGCTGAGCCTGAGCCTTGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCtt                                                                                           <  7:255888/147‑1 (MQ=255)
                                                       cagcagCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCtt                                                                                           <  6:148027/149‑1 (MQ=255)
                                                        agcagcGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTg                                                                                          >  4:181447/1‑149 (MQ=255)
                                                         gcagcGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTACTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTgg                                                                                         <  6:116531/149‑1 (MQ=255)
                                                              gcCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCg                                                                                    <  1:291986/149‑1 (MQ=255)
                                                               ccTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACAcc                                                                                             <  2:181460/139‑1 (MQ=255)
                                                               ccTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACAcc                                                                                             >  1:181460/1‑139 (MQ=255)
                                                                 tgctgGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCATGGCCG‑TTGCCCAGGTGCAGCAGAAAAAAGGCCCCCAGGATGCCCTGACGGATTTCCACACCTTGGTGCCGGGc                                                                                 >  5:162101/1‑149 (MQ=255)
                                                                   ctgGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGc                                                                                 <  6:234363/147‑1 (MQ=255)
                                                                   ctgGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGc                                                                                 <  6:267830/147‑1 (MQ=255)
                                                                      ggCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGc                                                                              <  3:181447/147‑1 (MQ=255)
                                                                       gCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGATTTCGACACCTTGGTGCCGGGCCGCTTc                                                                           <  4:248433/149‑1 (MQ=255)
                                                                           gACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCtt                                                                                           <  7:260265/129‑1 (MQ=255)
                                                                           gACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCtt                                                                                           >  8:260265/1‑129 (MQ=255)
                                                                            aCCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCa                                                                         >  6:300328/1‑146 (MQ=255)
                                                                            aCCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCa                                                                         <  5:300328/146‑1 (MQ=255)
                                                                            aCCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGAc                                                                      >  3:158497/1‑149 (MQ=255)
                                                                                   gcgcCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGc                                                                 <  7:236057/148‑1 (MQ=255)
                                                                                    cgcCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGTATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGAGCTTGCGt                                                               <  8:216135/149‑1 (MQ=255)
                                                                                    cgcCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGt                                                               >  5:278671/1‑149 (MQ=255)
                                                                                     gcCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGt                                                               >  6:66404/1‑148 (MQ=255)
                                                                                      ccGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTTACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGa                                                             <  1:223001/149‑1 (MQ=255)
                                                                                      ccGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGa                                                             <  4:235927/149‑1 (MQ=255)
                                                                                       cGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGAt                                                            <  3:37288/149‑1 (MQ=255)
                                                                                        ggCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATg                                                           <  4:209570/149‑1 (MQ=255)
                                                                                                tGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCgg                                                   >  1:89423/1‑149 (MQ=255)
                                                                                                  gTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGc                                                                                 <  4:42333/117‑1 (MQ=255)
                                                                                                      cccTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAAcc                                             >  1:178446/1‑149 (MQ=255)
                                                                                                      cccTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAAc                                              <  4:281302/148‑1 (MQ=255)
                                                                                                       ccTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGa                                                             <  3:37819/132‑1 (MQ=255)
                                                                                                       ccTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGa                                                             >  4:37819/1‑132 (MQ=255)
                                                                                                       ccTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAAcc                                             >  1:228890/1‑148 (MQ=255)
                                                                                                              gtggCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCg                                      <  6:3088/148‑1 (MQ=255)
                                                                                                               tggCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAg                                    <  5:66404/149‑1 (MQ=255)
                                                                                                                 gCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGc                                   <  2:228890/148‑1 (MQ=255)
                                                                                                                          tgggtcagccagGCCCTCCTGTTCG‑TTGCCCAGGTGCTGCAGATCATCAGCCAGCAGGATGCGCTGACCGAGTTCGAGACGTTGGTGCCGGGGCTGTTGCAGTCCTTGCGTGATGGATCGCGGGTAACCTATACGGAGGATCTGtg                            <  4:92548/140‑2 (MQ=255)
                                                                                                                           gcctgagcctgGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGAt                        <  7:162995/149‑1 (MQ=255)
                                                                                                                           gcctgagcctgGCCCCCCTGGGCG‑TAGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCAATACCGCGCAACTTTCGGAt                        >  6:12802/1‑149 (MQ=255)
                                                                                                                                agcctgGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGAcc                    <  2:89423/148‑1 (MQ=255)
                                                                                                                                      ctgcaccTGGGCA‑ACGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACATTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGACCGCATCaa             >  8:249239/6‑149 (MQ=255)
                                                                                                                                       ccccccTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGACCGCATCAAc            >  4:67932/1‑149 (MQ=255)
                                                                                                                                        cccccTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGACCGCATCAAcc           >  3:138282/1‑149 (MQ=255)
                                                                                                                                         ccccTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAACAACTGTCGGATGACCGCATCAACCt          <  6:278671/149‑1 (MQ=255)
                                                                                                                                          cccTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGa                                                                       <  7:249619/87‑1 (MQ=255)
                                                                                                                                          cccTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGa                                                                       >  8:249619/1‑87 (MQ=255)
                                                                                                                                           ccTGGGCA‑ACGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGACCGCATCAACCTgg        >  6:285823/1‑149 (MQ=255)
                                                                                                                                                 cG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGAAGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGACCGCATCAACCTGGGCgggg  >  1:205546/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CGGGCGCCTGTACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATG‑GCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCG‑TTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGACCGCATCAACCTGGGCGGGG  >  NC_002947/1120787‑1121079

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