Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,070,513 A→T 5.6% L125Q (CTG→CAG)  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,5130AT 5.6% 283.3 / 4.8 89L125Q (CTG→CAG) PP_4466TauD/TfdA family dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (43/41);  new base T (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 = 7.40e-01

CTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAG  >  NC_002947/5070370‑5070651
                                                                                                                                               |                                                                                                                                          
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                                                                                                                                            >  3:310301/1‑109 (MQ=255)
                                   tGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCgca                                                                                                                                            <  4:310301/109‑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                                                                                               <  6:45398/149‑1 (MQ=255)
                                        gCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGgtgt                                                                                               <  1:401557/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                                                                                                               >  1:386047/1‑116 (MQ=255)
                                                         gTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAg                                                                                                               <  2:386047/116‑1 (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                                                             >  1:30997/1‑149 (MQ=255)
                                                                          gTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGa                                                             >  7:366524/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:37126/149‑1 (MQ=255)
                                                                                tAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCtt                                                       <  2:30997/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                                                                                               <  8:208498/99‑1 (MQ=255)
                                                                                          aCATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGTAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGgtgt                                                                                               >  7:208498/1‑99 (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               <  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)
                                                                                                                                               |                                                                                                                                          
CTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAG  >  NC_002947/5070370‑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: