Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,070,513 A→T 7.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 7.6% 195.1 / 5.4 66L125Q (CTG→CAG) PP_4466TauD/TfdA family dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (27/34);  new base T (2/3);  total (29/37)
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.14e-01

AGACCTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTG  >  NC_002947/5070366‑5070661
                                                                                                                                                   |                                                                                                                                                    
aGACCTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCgcag                                                                                                                                                     >  8:204932/1‑149 (MQ=255)
        cgccgcGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTc                                                                                                                                             <  6:180683/149‑1 (MQ=255)
            gcgTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGc                                                                                                                                         <  7:255707/149‑1 (MQ=255)
            gcgTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCACTGCTGCGCTCCatt                                                                                                                                         <  5:159362/149‑3 (MQ=255)
             cgTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAg                                                                                                                                          >  1:37103/1‑147 (MQ=255)
             cgTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAg                                                                                                                                          <  2:37103/147‑1 (MQ=255)
             cgTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCt                                                                                                                                        >  8:188038/1‑149 (MQ=255)
                 aTACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCtggtg                                                                                                                                    <  5:332690/149‑1 (MQ=255)
                 aTACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCtggt                                                                                                                                     <  1:145731/148‑1 (MQ=255)
                 aTACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCtggt                                                                                                                                     <  2:98268/148‑1 (MQ=255)
                  tACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCtggtg                                                                                                                                    >  5:175504/1‑148 (MQ=255)
                   acacGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCACTGCTGCGCTCCatt                                                                                                                                         >  5:49887/1‑140 (MQ=255)
                   acacGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCACTGCTGCGCTCCatt                                                                                                                                         <  6:49887/142‑3 (MQ=255)
                    cacGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATAGATACCCAATAGGTGTCACCACCACTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGAc                                                                                                                                 >  6:246470/1‑149 (MQ=255)
                     acGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACg                                                                                                                                >  2:365641/1‑149 (MQ=255)
                     acGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACg                                                                                                                                >  6:238033/1‑149 (MQ=255)
                        ctctGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTc                                                                                                                             <  7:204932/149‑1 (MQ=255)
                             gaagaaGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTa                                                                                                                          <  3:172827/147‑1 (MQ=255)
                             gaagaaGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCCCCAGCTGGTGGACGCTc                                                                                                                             >  7:121271/1‑144 (MQ=255)
                             gaagaaGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCCCCAGCTGGTGGACGCTc                                                                                                                             <  8:121271/144‑1 (MQ=255)
                               agaagaTCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAAc                                                                                                                      >  2:233983/1‑149 (MQ=255)
                                gaagaTCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAAcc                                                                                                                     <  3:67456/149‑1 (MQ=255)
                                  agaTCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAg                                                                                                                   <  5:183185/149‑1 (MQ=255)
                                  agaTCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAg                                                                                                                   <  1:365641/149‑1 (MQ=255)
                                          cGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGgt                                                                                                           >  1:184970/1‑149 (MQ=255)
                                          cGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGgt                                                                                                           <  7:188038/149‑1 (MQ=255)
                                               tttGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCg                                                                                                             <  4:165915/142‑1 (MQ=255)
                                               tttGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCg                                                                                                             >  3:165915/1‑142 (MQ=255)
                                                            aGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGtcgcttcgct                                                                                         <  5:380725/149‑1 (MQ=255)
                                                                        ggCAGGGTCTCGTAAGCCTGATACATCGATAGCATATAGGTGTGACCTCCTGTGTGCGGAAGTGCAATTGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTTGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGc                                                                             <  5:238033/149‑1 (MQ=255)
                                                                           aGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGTTGTTACCCAAAGCAcc                                                                          <  1:343286/149‑1 (MQ=255)
                                                                            gggTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCAccc                                                                         >  3:5400/1‑149 (MQ=255)
                                                                                   gTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCt                                                                  >  8:51631/1‑149 (MQ=255)
                                                                                      aGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTa                                                                <  1:146691/148‑1 (MQ=255)
                                                                                       gCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGt                                                                                                                           >  5:15737/1‑88 (MQ=255)
                                                                                       gCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGt                                                                                                                           <  6:15737/88‑1 (MQ=255)
                                                                                            aTACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCtt                                                         <  1:233983/149‑1 (MQ=255)
                                                                                            aTACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCt                                                          <  4:5400/148‑1 (MQ=255)
                                                                                              aCATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGAt                                                                      >  4:80821/1‑134 (MQ=255)
                                                                                              aCATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCtt                                                         >  8:186751/1‑147 (MQ=255)
                                                                                              aCATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTAGCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGAt                                                                      <  3:80821/134‑1 (MQ=255)
                                                                                               cATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTta                                                       <  5:355605/148‑1 (MQ=255)
                                                                                               cATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTta                                                       <  2:184970/148‑1 (MQ=255)
                                                                                                aTCGATACCCAATAGGTGTCACCATCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTtag                                                      >  7:234344/1‑148 (MQ=255)
                                                                                                     tACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAgtgt                                                <  5:246470/149‑1 (MQ=255)
                                                                                                     tACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAgtgt                                                <  6:175504/149‑1 (MQ=255)
                                                                                                     tACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAgtgt                                                <  1:74113/149‑1 (MQ=255)
                                                                                                        ccAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCtggtgg                                                                                                                                   <  8:158230/63‑1 (MQ=255)
                                                                                                        ccAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCtggtgg                                                                                                                                   >  7:158230/1‑63 (MQ=255)
                                                                                                        ccAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAg                                             <  4:68897/149‑1 (MQ=255)
                                                                                                         cAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTta                                                       <  3:159239/138‑1 (MQ=255)
                                                                                                         cAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTta                                                       >  4:159239/1‑138 (MQ=255)
                                                                                                         cAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAg                                             >  2:291359/1‑148 (MQ=255)
                                                                                                          aaTAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTtagta                                                    <  1:77559/140‑1 (MQ=255)
                                                                                                          aaTAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTtagta                                                    >  2:77559/1‑140 (MQ=255)
                                                                                                                    accaccAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCaaaa                                 >  8:268713/1‑149 (MQ=255)
                                                                                                                      caccaGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAAtt                               >  3:6394/1‑149 (MQ=255)
                                                                                                                              ggCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCttttt                        <  8:234344/148‑1 (MQ=255)
                                                                                                                                       gCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGAt               <  7:186751/148‑1 (MQ=255)
                                                                                                                                       gCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGAt               >  8:146132/1‑148 (MQ=255)
                                                                                                                                        cAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGAt               >  6:169135/1‑147 (MQ=255)
                                                                                                                                         aaTGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCTGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAg            >  3:289008/1‑149 (MQ=255)
                                                                                                                                          ctGGAGCACTGCTGCGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCttttt                        >  7:310864/2‑136 (MQ=255)
                                                                                                                                          ctGGAGCACTGCTGCGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCttttt                        <  8:310864/135‑1 (MQ=255)
                                                                                                                                            ggAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTg                  <  6:340806/140‑1 (MQ=255)
                                                                                                                                            ggAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTg                  >  5:340806/1‑140 (MQ=255)
                                                                                                                                              aGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGccc       >  1:109884/1‑149 (MQ=255)
                                                                                                                                                   agcagTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGAt               <  6:128092/136‑1 (MQ=255)
                                                                                                                                                   agcagTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGAt               >  5:128092/1‑136 (MQ=255)
                                                                                                                                                   agcagTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGcc        <  4:129254/143‑1 (MQ=255)
                                                                                                                                                   agcagTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGcc        >  3:129254/1‑143 (MQ=255)
                                                                                                                                                   agcagTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTg  >  8:14101/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGACCTCCCGCCGCGTCATACACGCTCTGGAAGAAGATCTGCCGGCCTTTGACTGCGTCCAGTAGATAGGACGGCAGGGTCTCGTAAGCCTGATACATCGATACCCAATAGGTGTCACCACCAGTGGGCGGAAGTGCAATGGAGCGCAGCAGTGCTCCAGCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTG  >  NC_002947/5070366‑5070661

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: