Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,070,518 T→C 5.6% A123A (GCA→GCG PP_4466 ← TauD/TfdA family dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,070,5180TC 5.6% 264.4 / 3.0 88A123A (GCA→GCGPP_4466TauD/TfdA family dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (42/41);  new base C (2/3);  total (44/44)
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

CTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGC  >  NC_002947/5070370‑5070665
                                                                                                                                                    |                                                                                                                                                   
cTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGt                                                                                                                                                     >  2:307028/1‑149 (MQ=255)
cTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGt                                                                                                                                                     >  6:275827/1‑149 (MQ=255)
   ccgccgCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCt                                                                                                                                                  >  6:99728/1‑149 (MQ=255)
        gcgTCTTACACGCGCTGGAAGATGATCTGCCGGCCTTTGACTGCGTCCAGTTGTTAGGACGGCAGGGGCTCGTATGCCTGATACATCGATAGGCAATAGGTGTCACGACCAGTGGGCGGTAGTGGAATGGAGCACTGCTGCGCTCCatt                                                                                                                                             <  1:124610/149‑3 (MQ=255)
           tCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCtgg                                                                                                                                          >  5:411506/1‑149 (MQ=255)
           tCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCACTGCTGCGCTCCattgca                                                                                                                                          <  2:289697/149‑6 (MQ=255)
               acacGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGa                                                                                                                                      >  5:45398/1‑149 (MQ=255)
                cacGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCTTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGa                                                                                                                                      <  6:411506/148‑1 (MQ=255)
                     tctGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTg                                                                                                                                                    <  5:17339/129‑1 (MQ=255)
                     tctGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTg                                                                                                                                                    >  6:17339/1‑129 (MQ=255)
                       tGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCg                                                                                                                                >  3:149557/1‑147 (MQ=255)
                       tGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCg                                                                                                                                >  5:368632/1‑147 (MQ=255)
                         gaagaaGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAg                                                                                                                             >  6:309694/1‑148 (MQ=255)
                          aagaagATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCACTGCTGCGCTCCat                                                                                                                                              <  5:80116/130‑2 (MQ=255)
                          aagaagATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCACTGCTGCGCTCCat                                                                                                                                              >  6:80116/1‑129 (MQ=255)
                            gaagaTCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAAc                                                                                                                          >  2:105590/1‑148 (MQ=255)
                            gaagaTCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAAc                                                                                                                          >  2:122983/1‑148 (MQ=255)
                              agaTCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCa                                                                                                                        <  5:309694/148‑1 (MQ=255)
                              agaTCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAg                                                                                                                       <  4:149557/149‑1 (MQ=255)
                                  cTGCCGGCCTTTGACTGCGTCCAGGAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGGCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTAt                                                                                                                   <  8:6953/149‑1 (MQ=255)
                                    gCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCg                                                                                                                 >  2:401557/1‑149 (MQ=255)
                                     ccGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCg                                                                                                                 <  3:170193/148‑1 (MQ=255)
                                     ccGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCg                                                                                                                 <  1:122983/148‑1 (MQ=255)
                                       ggCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGTTACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAg                                                                                                                                              <  3:138999/117‑1 (MQ=255)
                                        gCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGgtgt                                                                                                             <  1:401557/149‑1 (MQ=255)
                                        gCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGgtgt                                                                                                             <  6:45398/149‑1 (MQ=255)
                                             tGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGcc                                                                                                          <  5:155662/147‑1 (MQ=255)
                                              gACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAgg                                                                                                       >  6:86793/1‑149 (MQ=255)
                                                   cGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGtcgc                                                                                                   <  5:86793/148‑1 (MQ=255)
                                                     tCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGtcgcttc                                                                                                >  2:37631/1‑149 (MQ=255)
                                                     tCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGtcgcttc                                                                                                <  1:105590/149‑1 (MQ=255)
                                                         gTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAg                                                                                                                             <  2:386047/116‑1 (MQ=255)
                                                         gTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAg                                                                                                                             >  1:386047/1‑116 (MQ=255)
                                                         gTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTg                                                                                            >  3:387470/1‑149 (MQ=255)
                                                            gATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTa                                                                                         <  3:297469/149‑1 (MQ=255)
                                                                  aCGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCaaa                                                                                   <  1:307028/149‑1 (MQ=255)
                                                                    ggCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAg                                                                                  >  5:274029/1‑148 (MQ=255)
                                                                      cAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCAc                                                                               <  1:29233/149‑1 (MQ=255)
                                                                          gTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGa                                                                           >  7:366524/1‑149 (MQ=255)
                                                                          gTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGa                                                                           >  1:30997/1‑149 (MQ=255)
                                                                          gTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGa                                                                           >  1:37126/1‑149 (MQ=255)
                                                                          gTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGa                                                                           <  1:378500/149‑1 (MQ=255)
                                                                           tctcGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCACGTCGCTTCGCTGTTACCCAAAGCACCCGAt                                                                          <  1:37631/149‑1 (MQ=255)
                                                                             tcGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATg                                                                         >  8:70573/1‑148 (MQ=255)
                                                                               gTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATg                                                                         <  7:70573/146‑1 (MQ=255)
                                                                               gTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCt                                                                      >  7:174800/1‑149 (MQ=255)
                                                                                tAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCtt                                                                     <  2:30997/149‑1 (MQ=255)
                                                                                tAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCtt                                                                     <  2:37126/149‑1 (MQ=255)
                                                                                tAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCAAGGTCGCTTCGCTGATACCCAAAGCACCCGATGGCt                                                                      >  8:268779/1‑148 (MQ=255)
                                                                                      tGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTc                                                               >  8:43084/1‑149 (MQ=255)
                                                                                         tACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCttt                                                            <  8:174800/149‑1 (MQ=255)
                                                                                         tACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCtt                                                             >  1:395915/1‑148 (MQ=255)
                                                                                          aCATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGTAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGgtgt                                                                                                             >  7:208498/1‑99 (MQ=255)
                                                                                          aCATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGTAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGgtgt                                                                                                             <  8:208498/99‑1 (MQ=255)
                                                                                           cATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGATACCCAAAGCACCCGATGGCTTACCGTCTTta                                                           <  7:268779/148‑1 (MQ=255)
                                                                                             tCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTtagta                                                        <  5:99728/149‑1 (MQ=255)
                                                                                                  aCCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAgtgt                                                    >  7:374743/1‑148 (MQ=255)
                                                                                                   cccAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGtt                                                   <  3:317177/148‑1 (MQ=255)
                                                                                                    ccAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAg                                                 <  2:395915/149‑1 (MQ=255)
                                                                                                    ccAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAg                                                 <  6:274029/149‑1 (MQ=255)
                                                                                                    ccAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAg                                                 >  4:148160/1‑149 (MQ=255)
                                                                                                     cAATACGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGc                                                <  3:6647/149‑1 (MQ=255)
                                                                                                          ggTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATaa                                           <  1:316681/149‑1 (MQ=255)
                                                                                                           gtgtCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAAc                                          >  1:315494/1‑149 (MQ=255)
                                                                                                                 ccaccaGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGc                                                <  7:45054/137‑1 (MQ=255)
                                                                                                                 ccaccaGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGc                                                >  8:45054/1‑137 (MQ=255)
                                                                                                                 ccaccaGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAAt                                    >  1:365416/1‑149 (MQ=255)
                                                                                                                 ccaccaGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAAt                                    >  7:307383/1‑149 (MQ=255)
                                                                                                                  caccaGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCa                                               >  2:352126/1‑137 (MQ=255)
                                                                                                                  caccaGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCa                                               <  1:352126/137‑1 (MQ=255)
                                                                                                                  caccaGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAAtt                                   >  3:8618/1‑149 (MQ=255)
                                                                                                                        tGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCtttt                             <  5:275827/149‑1 (MQ=255)
                                                                                                                        tGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCtttt                             >  3:352683/1‑149 (MQ=255)
                                                                                                                          ggCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTg                           <  8:366524/149‑1 (MQ=255)
                                                                                                                           gCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAGAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGc                          <  6:145366/149‑1 (MQ=255)
                                                                                                                           gCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAGAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGc                          >  5:145366/1‑149 (MQ=255)
                                                                                                                           gCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGc                          >  3:210465/1‑149 (MQ=255)
                                                                                                                               aaGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTg                      <  3:148160/149‑1 (MQ=255)
                                                                                                                                aGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTgg                     <  7:43084/149‑1 (MQ=255)
                                                                                                                                 gTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGa                    >  5:259661/1‑149 (MQ=255)
                                                                                                                                  tGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTta                                                           <  7:179161/109‑1 (MQ=255)
                                                                                                                                  tGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTta                                                           >  8:179161/1‑109 (MQ=255)
                                                                                                                                  tGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGAt                   >  2:146016/1‑149 (MQ=255)
                                                                                                                                    cAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTtagtag                                                       <  6:109409/111‑1 (MQ=255)
                                                                                                                                    cAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTtagtag                                                       >  5:109409/1‑111 (MQ=255)
                                                                                                                                     gctGGAGCACTGCTGCGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAg                >  2:136312/3‑149 (MQ=255)
                                                                                                                                     aaTGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAg                >  4:258262/1‑149 (MQ=255)
                                                                                                                                         gAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGc                                                >  4:41300/1‑113 (MQ=255)
                                                                                                                                         gAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGc                                                <  3:41300/113‑1 (MQ=255)
                                                                                                                                                   gTGCTCCAGCTGGTGGACGCTCGTAGAAGCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGc  >  1:233931/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGC  >  NC_002947/5070370‑5070665

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: