Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,120,934 T→C 7.2% 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 7.2% 197.1 / 4.0 69V125V (GTT→GTC) ‡PP_0982YjgP/YjgQ family permease
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base C (2/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGCGCCTGTACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGACCGCATCAACCTGGGCGG  >  NC_002947/1120787‑1121077
                                                                                                                                                    |                                                                                                                                               
ctggCTCCTGTACCTGGAAAGCGAGATGACCGTGCTGGCGGCCACTGTCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGtt                                                                                                                                                 <  6:127452/146‑1 (MQ=255)
 gggCGCCTGTACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGc                                                                                                                                               >  2:273113/1‑149 (MQ=255)
      ccTGTACCTGGAAAGCGCGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCAACGCCCAgg                                                                                                                                          <  2:295263/149‑1 (MQ=255)
      ccTGTACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAgg                                                                                                                                          >  4:45625/1‑149 (MQ=255)
       cTGTACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCAACGCCCAggg                                                                                                                                         <  4:279972/149‑2 (MQ=255)
          tACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTgc                                                                                                                                       <  5:143188/148‑1 (MQ=255)
            ccTGGATAGAGAGATGACCGTGCGGTGGGCCACTGGCATGAGCCATCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCGCTGGATGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGGCCAGGTgcag                                                                                                                                     <  6:270382/148‑1 (MQ=255)
            ccTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTgcag                                                                                                                                     >  6:198114/1‑148 (MQ=255)
             cTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTgcag                                                                                                                                     >  3:364988/1‑147 (MQ=255)
                 aaaGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGatc                                                                                                                               <  7:342512/149‑1 (MQ=255)
                 aaaGCGAGATGACCGTGCTGGCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGatca                                                                                                                              <  1:315664/149‑1 (MQ=255)
                  aaGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGatca                                                                                                                              >  3:13950/1‑149 (MQ=255)
                              cGTGCTGTCGGCCTCTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCTTGGCGCCGGCTGCCCTGGTCGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCGTGGCCCCCCTGGGCAACGCCCAgggg                                                                                                                                        <  5:271860/127‑3 (MQ=255)
                              cGTGCTGTCGGCCACTGGCATGAGACAGCAGCGCCTCCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCGTGCTGGTGGGCTGGCTGAGCCTGAGCCTGGCCCCCCTGCGCAACGCCCAgggg                                                                                                                                        >  6:271860/1‑125 (MQ=255)
                               gTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGa                                                                                                                 >  6:181884/1‑149 (MQ=255)
                                     tCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGcc                                                                                                             <  4:253807/147‑1 (MQ=255)
                                     tCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGGTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGccc                                                                                                            <  6:225381/148‑1 (MQ=255)
                                       ggCCACTGGCATGAGCCCGCAGCGCCAGCTGGGCATGACCATGGCGCCGGCGGCCCGGGTTTCCCAGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCAGGGGGTTGCCCAGGTGCAGCAGGTAAAAAGCCAGCAGGATTCCCTGa                                                                                                         >  8:71187/1‑149 (MQ=255)
                                        gCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGcccggcc                                                                                                        >  2:346157/1‑145 (MQ=255)
                                        gCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGAc                                                                                                        <  5:178024/149‑1 (MQ=255)
                                         ccACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGAcc                                                                                                       <  7:126325/149‑1 (MQ=255)
                                         ccACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGAcc                                                                                                       <  4:364988/149‑1 (MQ=255)
                                           aCTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGAcc                                                                                                       >  7:60443/1‑147 (MQ=255)
                                               gCATGAGCCAGCAGCGCCTGCTGGTCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGt                                                                                                   <  3:146209/147‑1 (MQ=255)
                                               gCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTc                                                                                                 <  5:99882/149‑1 (MQ=255)
                                                 aTGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCGTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGAGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGACTTCGa                                                                                               >  1:121608/1‑149 (MQ=255)
                                                      ccAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCt                                                                                          <  3:232607/149‑1 (MQ=255)
                                                      ccAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCt                                                                                          <  3:68709/149‑1 (MQ=255)
                                                      ccAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCt                                                                                          >  4:68709/1‑149 (MQ=255)
                                                      ccAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGCGCGTTGCCCAGGTGCAGCACATCATCAGACAGCAGGATGCCCTGCCCGAGTTCGACACCt                                                                                          >  8:179185/1‑149 (MQ=255)
                                                       cagcagCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCtt                                                                                         <  5:198114/149‑1 (MQ=255)
                                                        agcagcGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTg                                                                                        <  8:92482/149‑1 (MQ=255)
                                                         gcagcGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTg                                                                                        <  2:279674/149‑1 (MQ=255)
                                                                ctgctgGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTc                                                                                                 >  6:230772/1‑132 (MQ=255)
                                                                ctgctgGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTc                                                                                                 <  5:230772/132‑1 (MQ=255)
                                                                ctgctgGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCggg                                                                                >  3:30703/1‑149 (MQ=255)
                                                                        cATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTc                                                                         >  3:98979/1‑148 (MQ=255)
                                                                           gACCTTGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCt                                                                                          <  2:151831/128‑1 (MQ=255)
                                                                           gACCTTGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCt                                                                                          >  1:151831/1‑128 (MQ=255)
                                                                           gACCTTGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCt                                                                                          >  1:147400/1‑128 (MQ=255)
                                                                           gACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCTCTTCCAGa                                                                     >  8:327398/1‑149 (MQ=255)
                                                                           gACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGa                                                                     >  7:27137/1‑149 (MQ=255)
                                                                              cATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCt                                                                  >  2:34428/1‑149 (MQ=255)
                                                                               aTGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCtt                                                                 <  4:98979/149‑1 (MQ=255)
                                                                                      cGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGAt                                                          <  8:60443/149‑1 (MQ=255)
                                                                                       ggCTGCCCTGGGTGCCCTGCGGCGGGCCGGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGATCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATg                                                         <  2:121608/149‑1 (MQ=255)
                                                                                          tGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCt                                                      >  3:286517/1‑149 (MQ=255)
                                                                                              cTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTcgcg                                                  >  7:156733/1‑149 (MQ=255)
                                                                                               tGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCgg                                                 >  7:39066/1‑149 (MQ=255)
                                                                                                 gTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCgg                                                 <  4:286517/147‑1 (MQ=255)
                                                                                                     cccTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAAc                                            <  2:298109/148‑1 (MQ=255)
                                                                                                     cccTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAAc                                            <  7:71187/148‑1 (MQ=255)
                                                                                                     cccTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAAc                                            <  1:35211/148‑1 (MQ=255)
                                                                                                          ctgGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATAc                                      <  1:273113/149‑1 (MQ=255)
                                                                                                            ggtggCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATAcc                                     >  1:56792/1‑148 (MQ=255)
                                                                                                             gtggCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATAcc                                     >  8:17594/1‑147 (MQ=255)
                                                                                                             gtggCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGa                                   <  8:27137/149‑1 (MQ=255)
                                                                                                                         agcctgagcctgGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGCCCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGa                       >  4:176852/1‑149 (MQ=255)
                                                                                                                         agcctgagcctgGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGa                       <  4:108818/149‑1 (MQ=255)
                                                                                                                         agcctgagcctgGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCGTCCCGACCTTGCGGGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGa                       >  4:182658/1‑149 (MQ=255)
                                                                                                                          gcctgagcctgGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGAt                      <  5:181884/149‑1 (MQ=255)
                                                                                                                            ctgagcctgGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGa                    >  2:78252/1‑149 (MQ=255)
                                                                                                                            ctgagcctgGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGa                    >  2:191891/1‑149 (MQ=255)
                                                                                                                             tgagcctgGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGAc                   >  7:281987/1‑149 (MQ=255)
                                                                                                                               agcctgGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCTGGTAACCTATACCGAGCAACTGTCGGATGAc                   <  3:45625/147‑1 (MQ=255)
                                                                                                                                       ccccccTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGACCGCATCAAc          <  8:156733/149‑1 (MQ=255)
                                                                                                                                       ccccccTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGACCGCATCAAc          <  4:13950/149‑1 (MQ=255)
                                                                                                                                         ccccTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGCCCGAGTTCGACACCGTGGTGCCGGGCCGCGTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATAACGAGCAACTGTCGGATGCCCGCAGCATccg        >  6:229721/1‑148 (MQ=255)
                                                                                                                                         caccTGGGCAACGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGACCGCATCAACCt        >  8:345331/3‑149 (MQ=255)
                                                                                                                                               ggCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGACCGCATCAACCTGGGCgg  <  2:56792/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CGGGCGCCTGTACCTGGAAAGCGAGATGACCGTGCTGTCGGCCACTGGCATGAGCCAGCAGCGCCTGCTGGGCATGACCATGGCGCCGGCTGCCCTGGTTGCCCTGCTGGTGGCCTGGCTGAGCCTGAGCCTGGCCCCCCTGGGCGTTGCCCAGGTGCAGCAGATCATCAGCCAGCAGGATGCCCTGACCGAGTTCGACACCTTGGTGCCGGGCCGCTTCCAGACCTTGCGTGATGGCTCGCGGGTAACCTATACCGAGCAACTGTCGGATGACCGCATCAACCTGGGCGG  >  NC_002947/1120787‑1121077

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: