Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,070,511 G→A 5.6% R126C (CGC→TGC)  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,5110GA 5.6% 275.4 / 3.8 89R126C (CGC→TGC) PP_4466TauD/TfdA family dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (43/41);  new base A (2/3);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GCAGACCTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAG  >  NC_002947/5070364‑5070651
                                                                                                                                                   |                                                                                                                                            
gCAGACCTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGTGCGGAAGTGCAATGGAgcgc                                                                                                                                             <  5:209169/149‑1 (MQ=255)
gCAGACCTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAgcg                                                                                                                                              <  4:27295/148‑1 (MQ=255)
gCAGACCTCCCGACGCGTAATACACGCGCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAgcgc                                                                                                                                             <  4:281754/149‑1 (MQ=255)
      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                                                                                                                                      >  6:17339/1‑129 (MQ=255)
                           tctGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTg                                                                                                                                      <  5:17339/129‑1 (MQ=255)
                             tGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCg                                                                                                                  >  5:368632/1‑147 (MQ=255)
                             tGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCg                                                                                                                  >  3:149557/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)
                                         tGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCgca                                                                                                                                            <  4:310301/109‑1 (MQ=255)
                                         tGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCgca                                                                                                                                            >  3:310301/1‑109 (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:378500/149‑1 (MQ=255)
                                                                                gTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGa                                                             >  1:37126/1‑149 (MQ=255)
                                                                                gTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGa                                                             >  1:30997/1‑149 (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                                 <  1:352126/137‑1 (MQ=255)
                                                                                                                        caccaGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCa                                 >  2:352126/1‑137 (MQ=255)
                                                                                                                        caccaGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAAtt                     >  3:8618/1‑149 (MQ=255)
                                                                                                                              tGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCtttt               >  3:352683/1‑149 (MQ=255)
                                                                                                                              tGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCtttt               <  5:275827/149‑1 (MQ=255)
                                                                                                                                ggCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTg             <  8:366524/149‑1 (MQ=255)
                                                                                                                                 gCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAGAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGc            >  5:145366/1‑149 (MQ=255)
                                                                                                                                 gCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAGAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGc            <  6:145366/149‑1 (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                                  <  3:41300/113‑1 (MQ=255)
                                                                                                                                               gAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGc                                  >  4:41300/1‑113 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
GCAGACCTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAG  >  NC_002947/5070364‑5070651

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: