Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,012,598 C→G 6.7% G475A (GGC→GCC)  PP_1793 ← group 2 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,012,5980CG 6.7% 241.6 / 5.9 75G475A (GGC→GCC) PP_1793group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (32/37);  new base G (2/3);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TGTTTTTTCCGCCCCGAAAGAAATCCCGCCTTGATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGCTGCCTGCTCACGAGCATTGGGCATGGCATGGCTGACTTGGCTTT  >  NC_002947/2012450‑2012745
                                                                                                                                                    |                                                                                                                                                   
tgtTTTTTCCGCCCCGAAAGAAATCCCGCCTTGATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGc                                                                                                                                                     >  7:33233/1‑149 (MQ=255)
 gtTTTTTCCGCCCCGAAAGAAATCCCGCCTTGATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGcc                                                                                                                                                    >  6:174465/1‑149 (MQ=255)
       tCCGCCCCGAAAGAAATCCCGCCTTGATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTg                                                                                                                                              >  1:192648/1‑149 (MQ=255)
         cGCCCCGAAAGAAATCCCGCCTTGATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTgg                                                                                                                                             <  5:150865/148‑1 (MQ=255)
          gCCCCGAAAGAAATCCCGCCTTGATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTCGGCTACCTggat                                                                                                                                           <  7:225470/149‑3 (MQ=255)
          gCCCCGAAAGAAATCCCGCCTTGATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTCGGCTACCTggat                                                                                                                                           <  1:57493/149‑3 (MQ=255)
             ccGAATGAAATCCCGCCTTGATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTCGGCTACCTGGATGaa                                                                                                                                        <  6:55802/149‑2 (MQ=255)
              cGAAAGAAATCCCGCCTTGATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAgc                                                                                                                                       >  3:346939/1‑149 (MQ=255)
                 aagaaaTCCCGCCTTGATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCt                                                                                                                                    <  7:255724/149‑1 (MQ=255)
                        cccGCCTTGATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCtt                                                                                                                            <  1:92740/149‑1 (MQ=255)
                        cccGCCTTGATGCGACTCGAATACATTATGCGCGCCCTTGTTGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGGTTCTTCGCCATAACCTCTGGCCAGGTAAACATCATCCCGGTAGCCGACCTGGTTGAGCGCTGCACGCt                                                                                                                             <  7:176132/149‑1 (MQ=255)
                         ccGCCTTGATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCtt                                                                                                                            >  2:95753/1‑149 (MQ=255)
                         ccGCCTTGATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGGTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCtt                                                                                                                            >  8:162882/1‑149 (MQ=255)
                          cGCCTTGATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCttt                                                                                                                           <  2:326159/149‑1 (MQ=255)
                             cTTGATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTtata                                                                                                                        <  1:95753/149‑1 (MQ=255)
                                gATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTtatata                                                                                                                      >  5:353102/1‑148 (MQ=255)
                                gATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTtatata                                                                                                                      <  3:340926/148‑1 (MQ=255)
                                 aTGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATac                                                                                                                     <  6:353102/148‑1 (MQ=255)
                                 aTGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATac                                                                                                                     <  7:317508/148‑1 (MQ=255)
                                       cACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATacagacag                                                                                                               >  6:76555/1‑148 (MQ=255)
                                          gAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGc                                                                                                           <  4:346939/149‑1 (MQ=255)
                                           aaTACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGcc                                                                                                          >  4:150199/1‑149 (MQ=255)
                                            aTACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCg                                                                                                         >  8:22552/1‑149 (MQ=255)
                                                     gcgcgcCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGtc                                                                                                >  6:125692/1‑149 (MQ=255)
                                                       gcgcCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGACGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGACGCAGCCAGtctc                                                                                              >  8:272531/1‑149 (MQ=255)
                                                         gcCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAAAAGCCGCAGCCa                                                                                                   <  4:260927/142‑1 (MQ=255)
                                                         gcCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAAAAGCCGCAGCCa                                                                                                   >  3:260927/1‑142 (MQ=255)
                                                          cccATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGtctcgat                                                                                           >  3:261852/1‑149 (MQ=255)
                                                          cccATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGtctcgat                                                                                           <  1:120639/149‑1 (MQ=255)
                                                          cccATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGtctcgat                                                                                           <  8:244496/149‑1 (MQ=255)
                                                          cccATGGGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGGTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGtctcgat                                                                                           <  4:219943/149‑1 (MQ=255)
                                                           ccATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGtctcgatc                                                                                          >  2:386847/1‑149 (MQ=255)
                                                             aTGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGtctcgatctc                                                                                        >  2:29640/1‑149 (MQ=255)
                                                                    gccagccaAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACgg                                                                                   >  4:75350/1‑147 (MQ=255)
                                                                      cagccaAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGgaga                                                                               <  5:76555/149‑1 (MQ=255)
                                                                       agccaAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGgagag                                                                              <  1:386847/149‑1 (MQ=255)
                                                                        gccaAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGga                                                                                 >  3:93036/1‑145 (MQ=255)
                                                                        gccaAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGga                                                                                 >  5:77970/1‑145 (MQ=255)
                                                                           aaaGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCGTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGt                                                                          >  3:41662/1‑149 (MQ=255)
                                                                            aaGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGtt                                                                         <  6:77970/149‑1 (MQ=255)
                                                                             aGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGCGCGCTGCACGCTTTATATACAGAAAGAAGCCGCAGCAAGTCTCGATCTAGACGGGAGAGATGTTg                                                                        >  8:279914/1‑149 (MQ=255)
                                                                                   gcgcACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGtctcgatctcg                                                                                       <  1:207464/128‑1 (MQ=255)
                                                                                   gcgcACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGtctcgatctcg                                                                                       <  1:201135/128‑1 (MQ=255)
                                                                                   gcgcACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGtctcgatctcg                                                                                       >  2:207464/1‑128 (MQ=255)
                                                                                   gcgcACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGtctcgatctcg                                                                                       >  2:201135/1‑128 (MQ=255)
                                                                                        cgcaGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGc                                                             <  5:174465/149‑1 (MQ=255)
                                                                                         gcaGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGcc                                                            >  6:226113/1‑149 (MQ=255)
                                                                                           aGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAg                                                          <  6:52620/149‑1 (MQ=255)
                                                                                              cACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTc                                                       <  4:206936/149‑1 (MQ=255)
                                                                                                 cAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATc                                                    <  5:125692/149‑1 (MQ=255)
                                                                                                  aGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGt                                                         <  5:153854/143‑1 (MQ=255)
                                                                                                  aGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGt                                                         >  6:153854/1‑143 (MQ=255)
                                                                                                    tCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAg                                                 <  5:226113/149‑1 (MQ=255)
                                                                                                      ggTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAgcg                                               >  7:272631/1‑149 (MQ=255)
                                                                                                        tttCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGct                                             <  8:44539/149‑1 (MQ=255)
                                                                                                        tttCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGct                                             <  8:50649/149‑1 (MQ=255)
                                                                                                           cttcTCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGctgc                                           <  7:279914/148‑1 (MQ=255)
                                                                                                            ttcGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGctgcct                                         >  5:269255/1‑149 (MQ=255)
                                                                                                            ttcGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGctgcct                                         >  6:185716/1‑149 (MQ=255)
                                                                                                             tcGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGtc                                                                                                <  1:381381/93‑1 (MQ=255)
                                                                                                             tcGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGtc                                                                                                >  2:381381/1‑93 (MQ=255)
                                                                                                               gCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGctgcctgc                                       <  3:150199/148‑1 (MQ=255)
                                                                                                                    aaCCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGCTGCCTGCTc                                     <  4:305624/145‑1 (MQ=255)
                                                                                                                    aaCCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGCTGCCTGCTc                                     >  3:305624/1‑145 (MQ=255)
                                                                                                                     aCCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGCTGCCTGCTCACGa                                 <  2:192648/148‑1 (MQ=255)
                                                                                                                     aCCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGCTGCCTGCTCACGa                                 <  7:22552/148‑1 (MQ=255)
                                                                                                                             ccAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGCTGCCTGCTCACGAGCATTggg                         >  1:193445/1‑148 (MQ=255)
                                                                                                                                   aaaCTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGCTGCCTGCTCACGAGCATTGggcat                      <  2:210185/145‑1 (MQ=255)
                                                                                                                                   aaaCTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGCTGCCTGCTCACGAGCATTGggcat                      >  1:210185/1‑145 (MQ=255)
                                                                                                                                   aaaCTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGCTGCCTGCTAACGAGCATTGggcatggca                  >  1:269047/1‑149 (MQ=255)
                                                                                                                                    aaCTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGCTGCCTGCTCACGAGCATTGggcatggcat                 <  5:185716/149‑1 (MQ=255)
                                                                                                                                        tCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGCTGCCTGCTCACGAGCATTGGGCATGGCATGGCt             >  1:323077/1‑149 (MQ=255)
                                                                                                                                        tCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGCTGCCTGCTCACGAGCATTGGGCATGGCATGGCt             >  1:318528/1‑149 (MQ=255)
                                                                                                                                        tCAACCAGGTCGGCTACCTGGTTGAGAGCTGCCCGATTTATATACCGACAGCAGCCGCAGCAAAGCAAGATCTCGACGGGAGAGATGGTGGCATGTTTGTCCAGTTCATCCAGCGCTGCCTGCTCACGAGCATTGGGCATGGCATGGCt             >  1:311433/1‑149 (MQ=255)
                                                                                                                                        tCAACCAGGTCGGCTACCTGGATGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGCTGCCTGCTCACGAGCATTGGGCATGGCATGGCt             >  3:173628/1‑149 (MQ=255)
                                                                                                                                             cAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGCTGCCTGCTCACGAGCATTGGGCATGGCATGGCTGACtt        <  3:93165/149‑1 (MQ=255)
                                                                                                                                                   gCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGCTGCCTGCTCACGAGCATTGGGCATGGCATGGCTGACTTGGCttt  <  8:272631/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGTTTTTTCCGCCCCGAAAGAAATCCCGCCTTGATGCGCCACGAATACATTAGGCGCGCCCATGTGGCGCCAGCCAAAGCTGCGCGCACGCAGGCACCAGTCGGTTTCTTCGCCATAACCTCTGGCCAGGTAAACTTCATCCAGGTAGCCGACCTGGTTGAGCGCTGCACGCTTTATATACAGACAGAAGCCGCAGCCAGTCTCGATCTCGACGGGAGAGATGTTGGCATGTTTGGCCAGTTCATCCAGCGCTGCCTGCTCACGAGCATTGGGCATGGCATGGCTGACTTGGCTTT  >  NC_002947/2012450‑2012745

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: