Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,070,511 G→A 5.0% 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.0% 243.6 / 2.4 79R126C (CGC→TGC) PP_4466TauD/TfdA family dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (44/31);  new base A (2/2);  total (46/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GCAGACCTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCT  >  NC_002947/5070364‑5070657
                                                                                                                                                   |                                                                                                                                                  
gCAGACCTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAgcg                                                                                                                                                    >  3:84578/1‑148 (MQ=255)
      cTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGt                                                                                                                                             >  1:79448/1‑149 (MQ=255)
      cTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGt                                                                                                                                             >  4:314564/1‑149 (MQ=255)
        cccGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGc                                                                                                                                           >  2:41374/1‑149 (MQ=255)
             cgcgTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCat                                                                                                                                      <  3:314564/149‑2 (MQ=255)
               cgTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCt                                                                                                                                    >  4:216775/1‑149 (MQ=255)
                     acacGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCtggt                                                                                                                                 <  6:170460/146‑1 (MQ=255)
                     acacGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCtggt                                                                                                                                 >  5:170460/1‑146 (MQ=255)
                      cacGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGa                                                                                                                              <  2:79448/148‑1 (MQ=255)
                       acGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACg                                                                                                                            >  8:230381/1‑149 (MQ=255)
                         gCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCt                                                                                                                          >  2:360362/1‑149 (MQ=255)
                         gCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCt                                                                                                                          >  2:345215/1‑149 (MQ=255)
                                  gaagaTCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGAAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAAc                                                                                                                  >  4:31421/1‑148 (MQ=255)
                                    agaTCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGTAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAg                                                                                                               <  3:31421/149‑1 (MQ=255)
                                           ccGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGTGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCg                                                                                                         <  5:15067/148‑1 (MQ=255)
                                           ccGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATc                                                                                                          <  1:360362/147‑1 (MQ=255)
                                                cTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGc                                                                                                   >  5:176864/1‑149 (MQ=255)
                                                    gACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAgg                                                                                               >  7:6621/1‑149 (MQ=255)
                                                      cTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGtc                                                                                             <  1:345215/149‑1 (MQ=255)
                                                       tGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGtcg                                                                                            >  3:258812/1‑149 (MQ=255)
                                                        gCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGtcgc                                                                                           >  7:65644/1‑149 (MQ=255)
                                                         cGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGtcgct                                                                                          <  1:41374/149‑1 (MQ=255)
                                                             cAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGtcgcttcgc                                                                                      >  4:200081/1‑149 (MQ=255)
                                                             cAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGtcgcttcgc                                                                                      <  4:258812/149‑1 (MQ=255)
                                                              aGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGtcgcttcgct                                                                                     <  4:84578/149‑1 (MQ=255)
                                                               gTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTg                                                                                    >  8:1681/1‑149 (MQ=255)
                                                               gTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTg                                                                                    >  4:125878/1‑149 (MQ=255)
                                                                 aGATAGGACGGCAGGTTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCTGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGt                                                                                   <  5:166578/148‑1 (MQ=255)
                                                                  gATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTa                                                                                 >  3:69046/1‑149 (MQ=255)
                                                                   aTAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGCAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTAc                                                                                <  7:63442/149‑1 (MQ=255)
                                                                     aGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCg                                                                                                         >  1:214474/1‑122 (MQ=255)
                                                                     aGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCg                                                                                                         <  2:214474/122‑1 (MQ=255)
                                                                           gCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCa                                                                        >  4:99922/1‑149 (MQ=255)
                                                                            cAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGtcgcttcgct                                                                                     <  7:333254/135‑1 (MQ=255)
                                                                            cAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGtcgcttcgct                                                                                     >  8:333254/1‑135 (MQ=255)
                                                                            cAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCAc                                                                       >  3:101991/1‑149 (MQ=255)
                                                                                gTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGa                                                                   >  3:10830/1‑149 (MQ=255)
                                                                                   tcGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATgg                                                                <  7:230381/149‑1 (MQ=255)
                                                                                     gTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGc                                                               <  5:378697/148‑1 (MQ=255)
                                                                                      tAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCaaa                                                                           >  2:99846/1‑135 (MQ=255)
                                                                                      tAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCaaa                                                                           <  1:99846/135‑1 (MQ=255)
                                                                                          ccTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCAAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCATGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCg                                                         >  1:301960/1‑149 (MQ=255)
                                                                                            tGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTc                                                       >  6:322245/1‑149 (MQ=255)
                                                                                             gATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCt                                                      >  2:246377/1‑149 (MQ=255)
                                                                                                   tCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTtagta                                                <  3:367677/149‑1 (MQ=255)
                                                                                                      aTACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAgtg                                             >  4:178978/1‑149 (MQ=255)
                                                                                                       tACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAgtgt                                            <  4:69046/149‑1 (MQ=255)
                                                                                                       tACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAgtgt                                            <  3:99922/149‑1 (MQ=255)
                                                                                                         cccAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGtt                                           <  3:216775/148‑1 (MQ=255)
                                                                                                          ccAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTta                                                   >  6:227753/1‑139 (MQ=255)
                                                                                                          ccAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTta                                                   <  5:227753/139‑1 (MQ=255)
                                                                                                           cAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAg                                         >  4:108179/1‑148 (MQ=255)
                                                                                                            aaTAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGc                                        >  2:115757/1‑148 (MQ=255)
                                                                                                            aaTAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGc                                        >  4:48047/1‑148 (MQ=255)
                                                                                                              tAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCaa                                      >  2:128214/1‑148 (MQ=255)
                                                                                                                ggTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATaa                                   <  3:178978/149‑1 (MQ=255)
                                                                                                                 gtgtCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAAc                                  >  6:213923/1‑149 (MQ=255)
                                                                                                                         accaGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAAttt                          >  8:20979/1‑149 (MQ=255)
                                                                                                                           caGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCt                        <  3:48047/149‑1 (MQ=255)
                                                                                                                             gTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATaa                                   <  4:350806/136‑1 (MQ=255)
                                                                                                                             gTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATaa                                   >  3:350806/1‑136 (MQ=255)
                                                                                                                                 gCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGc                  >  8:161796/1‑149 (MQ=255)
                                                                                                                                 gCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGc                  <  1:246377/149‑1 (MQ=255)
                                                                                                                                   ggAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCaaaa                             <  2:50630/136‑1 (MQ=255)
                                                                                                                                   ggAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCaaaa                             >  1:50630/1‑136 (MQ=255)
                                                                                                                                   ggAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGt                >  4:281358/1‑149 (MQ=255)
                                                                                                                                    gAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGtt               <  7:20979/149‑1 (MQ=255)
                                                                                                                                    gAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGtt               <  3:125878/149‑1 (MQ=255)
                                                                                                                                     aaGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTg              <  1:128214/149‑1 (MQ=255)
                                                                                                                                        tGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGAt           >  7:348658/1‑149 (MQ=255)
                                                                                                                                          cAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGAt           >  3:254269/1‑147 (MQ=255)
                                                                                                                                           aaTGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAg        >  3:164500/1‑149 (MQ=255)
                                                                                                                                           aaTGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAg        >  4:184865/1‑149 (MQ=255)
                                                                                                                                            ctGGAGCACTGCTGCGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCttttt                    <  1:362435/135‑1 (MQ=255)
                                                                                                                                            ctGGAGCACTGCTGCGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCttttt                    >  2:362435/2‑136 (MQ=255)
                                                                                                                                            ctGGAGCACTGCTGCGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGt                >  6:196890/2‑140 (MQ=255)
                                                                                                                                            ctGGAGCACTGCTGCGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGt                <  5:196890/139‑1 (MQ=255)
                                                                                                                                            aTGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAgc       <  6:176864/149‑1 (MQ=255)
                                                                                                                                               gAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGcc    >  7:67543/1‑149 (MQ=255)
                                                                                                                                                aGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTGTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGGGTTGGATGAGCGccc   >  4:60805/1‑149 (MQ=255)
                                                                                                                                                 gcgcAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCct  <  8:65644/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCAGACCTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCT  >  NC_002947/5070364‑5070657

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