Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,216,060 C→T 5.0% A139A (GCG→GCA yehW ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,216,0600CT 5.0% 259.5 / 2.6 79A139A (GCG→GCAyehWputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (41/34);  new base T (2/2);  total (43/36)
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.17e-01

ATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGG  >  NC_000913/2215914‑2216205
                                                                                                                                                  |                                                                                                                                                 
aTCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGc                                                                                                                                                 >  1:352831/1‑149 (MQ=255)
   gtcgCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTcc                                                                                                                                              <  3:145532/149‑1 (MQ=255)
   gtcgCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTcc                                                                                                                                              <  2:35065/149‑1 (MQ=255)
   gtcgCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTcc                                                                                                                                              >  1:35065/1‑149 (MQ=255)
              aaTGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGc                                                                                                                                     >  3:48485/1‑147 (MQ=255)
              aaTGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCaa                                                                                                                                   <  6:107480/149‑1 (MQ=255)
              aaTGTTGATAATCACCGAAGTTCGCACGCCCGACAGAATCACCGGAGCCGACAGCGGTAGCTCGACCTGACGCACTCGCTGCCCACGACTCATTCCCATAACTTTCGCAACTTCTGGACCGCTGGCAGCAATCGCTCCCAGCCCGGCaa                                                                                                                                   >  5:236037/1‑149 (MQ=255)
                tGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAgt                                                                                                                                 <  4:60708/149‑1 (MQ=255)
                 gTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGTTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCCGCCCGGCAAgtg                                                                                                                                >  3:12954/1‑149 (MQ=255)
                    gATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTc                                                                                                                              <  3:58862/148‑1 (MQ=255)
                         tCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTg                                                                                                                         >  3:18161/1‑148 (MQ=255)
                         tCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTg                                                                                                                         >  1:120872/1‑148 (MQ=255)
                           aCCGAAGTTCGCCCGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCCCGCCTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAg                                                                                                                      >  5:234133/1‑149 (MQ=255)
                           aCCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAg                                                                                                                      >  8:32698/1‑149 (MQ=255)
                           aCCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAg                                                                                                                      >  7:343082/1‑149 (MQ=255)
                           aCCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAg                                                                                                                      >  6:195479/1‑149 (MQ=255)
                             cGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGCCCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGa                                                                                                                    >  1:175582/1‑149 (MQ=255)
                             cGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGCCCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGa                                                                                                                    >  1:69930/1‑149 (MQ=255)
                                 gTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACggg                                                                                                                >  6:206158/1‑149 (MQ=255)
                                 gTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACggg                                                                                                                >  6:66912/1‑149 (MQ=255)
                                  ttCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGgc                                                                                                               >  3:295816/1‑149 (MQ=255)
                                       aCGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGcac                                                                                                          <  2:174363/149‑1 (MQ=255)
                                       aCGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGcac                                                                                                          >  4:187178/1‑149 (MQ=255)
                                           ccgccAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCg                                                                                                      >  4:178287/1‑149 (MQ=255)
                                            cgccAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCa                                                                                                                       <  4:326151/131‑1 (MQ=255)
                                            cgccAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCa                                                                                                                       >  3:326151/1‑131 (MQ=255)
                                             gccAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGt                                                                                                     >  3:3380/1‑148 (MQ=255)
                                               cAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTa                                                                                                    >  4:19003/1‑147 (MQ=255)
                                                      aCCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCaa                                                                                           <  2:352831/149‑1 (MQ=255)
                                                      aCCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCa                                                                                            <  8:142428/148‑1 (MQ=255)
                                                        cGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAgg                                                                                         <  1:313234/149‑1 (MQ=255)
                                                        cGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAgg                                                                                         <  6:234133/149‑1 (MQ=255)
                                                         ggAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGc                                                                                        <  2:96665/149‑1 (MQ=255)
                                                               gccAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAAt                                                                                  >  3:166633/1‑149 (MQ=255)
                                                                ccAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATTGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATc                                                                                 >  5:167798/1‑149 (MQ=255)
                                                                ccAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATTGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATc                                                                                 >  7:268479/1‑149 (MQ=255)
                                                                  aGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGc                                                                               <  2:120872/149‑1 (MQ=255)
                                                                  aGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGc                                                                               <  8:343082/149‑1 (MQ=255)
                                                                   gCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCt                                                                              <  5:195479/149‑1 (MQ=255)
                                                                    cGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTTGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTg                                                                             <  7:136888/149‑1 (MQ=255)
                                                                      gTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGt                                                                           <  2:14775/149‑1 (MQ=255)
                                                                          cTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCa                                                                       <  3:19003/149‑1 (MQ=255)
                                                                           tCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGc                                                                               >  6:31006/1‑140 (MQ=255)
                                                                           tCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGc                                                                               <  5:31006/140‑1 (MQ=255)
                                                                           tCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAg                                                                      >  3:16721/1‑149 (MQ=255)
                                                                            cGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGa                                                                     <  7:67380/149‑1 (MQ=255)
                                                                             gACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAg                                                                      <  5:206158/147‑1 (MQ=255)
                                                                             gACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAg                                                                      <  5:66912/147‑1 (MQ=255)
                                                                                cTTACGCACTCGCTGCCCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCTCTGGTTGCAGAACa                                                                  >  2:206224/1‑148 (MQ=255)
                                                                                       aCTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGAt                                                          >  4:279871/1‑149 (MQ=255)
                                                                                           gCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCAcc                                                      <  8:95922/149‑1 (MQ=255)
                                                                                           gCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCAcc                                                      <  4:295816/149‑1 (MQ=255)
                                                                                               aCCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCCAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGaa                                                  <  4:18161/149‑1 (MQ=255)
                                                                                                 cACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTg                                                                         <  1:154452/124‑1 (MQ=255)
                                                                                                 cACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTg                                                                         >  2:154452/1‑124 (MQ=255)
                                                                                                 cACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACg                                                >  4:227911/1‑149 (MQ=255)
                                                                                                      cTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGc                                                              >  6:173452/1‑130 (MQ=255)
                                                                                                         aTTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCGCCGTTCCGGTGATCg                                         <  4:3380/148‑1 (MQ=255)
                                                                                                         aTTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGc                                        >  3:228741/1‑149 (MQ=255)
                                                                                                             ccATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGc                                                                                                                            <  6:266633/61‑1 (MQ=255)
                                                                                                             ccATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAg                                     >  6:173140/1‑148 (MQ=255)
                                                                                                              cATACCTTTCGCAACTTCTGTCACGCTGGCATCAATTGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCa                                   <  6:167798/149‑1 (MQ=255)
                                                                                                              cATACCTTTCGCAACTTCTGTCACGCTGGCATCAATTGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCa                                   <  8:268479/149‑1 (MQ=255)
                                                                                                                  ccTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGc                               >  1:52486/1‑149 (MQ=255)
                                                                                                                    tttCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCATGTGTCGCCTGCAGGACGGGCAGTACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCg                                         <  8:297560/137‑1 (MQ=255)
                                                                                                                    tttCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCATGTGTCGCCTGCAGGACGGGCAGTACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCg                                         >  7:297560/1‑137 (MQ=255)
                                                                                                                      tCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGAc                            <  2:175582/148‑1 (MQ=255)
                                                                                                                         cAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGc                        >  6:239716/1‑149 (MQ=255)
                                                                                                                            cTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCgg                      >  5:315042/1‑148 (MQ=255)
                                                                                                                                tGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCa                                                                       >  6:97172/1‑95 (MQ=255)
                                                                                                                                  tCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGaaaa                  <  1:247839/146‑1 (MQ=255)
                                                                                                                                  tCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGaaaa                  <  4:48485/146‑1 (MQ=255)
                                                                                                                                  tCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGGCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGt                <  7:272208/148‑1 (MQ=255)
                                                                                                                                    aCGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTg             >  8:151635/1‑149 (MQ=255)
                                                                                                                                     cGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGt            <  5:239716/149‑1 (MQ=255)
                                                                                                                                          gCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCCAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGCCGGGCGGAAAAGTCTGTCCAAc       >  1:69859/1‑149 (MQ=255)
                                                                                                                                            aTCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAAcgg     >  5:26143/1‑149 (MQ=255)
                                                                                                                                              cAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAAcggcg   >  4:226107/1‑149 (MQ=255)
                                                                                                                                                aTCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAAcggcgg  >  8:163275/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTTACGCACTCGCTGACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGG  >  NC_000913/2215914‑2216205

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: