Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,530,631 T→C 5.0% T341A (ACT→GCT)  sgcX ← putative endoglucanase with Zn‑dependent exopeptidase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,530,6310TC 5.0% 245.8 / 2.6 79T341A (ACT→GCT) sgcXputative endoglucanase with Zn‑dependent exopeptidase domain
Reads supporting (aligned to +/- strand):  ref base T (42/33);  new base C (2/2);  total (44/35)
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

GTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACC  >  NC_000913/4530487‑4530773
                                                                                                                                                |                                                                                                                                              
gTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCaaa                                                                                                                                            >  5:106833/1‑149 (MQ=255)
gTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCaaa                                                                                                                                            <  7:97359/149‑1 (MQ=255)
      gTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTc                                                                                                                                               <  8:367471/140‑1 (MQ=255)
      gTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTc                                                                                                                                               >  7:367471/1‑140 (MQ=255)
      gTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGCCAAATCACg                                                                                                                                       >  7:103221/1‑148 (MQ=255)
      gTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGCCAAATCACGc                                                                                                                                      >  5:31896/1‑149 (MQ=255)
       tCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCa                                                                                                                                     <  3:343527/149‑1 (MQ=255)
              gCCGGTACCGCATGCCACAAGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTgg                                                                                                                               >  6:144756/1‑148 (MQ=255)
                cGGTACCGCATGCCACAAGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCaa                                                                                                                            >  7:10826/1‑149 (MQ=255)
                  gTACCGCATGCCACAAGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAAcc                                                                                                                          >  2:22942/1‑149 (MQ=255)
                            ccACAAGGATCTTTTTCATAATGGATGTGACTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCTGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGga                                                                                                                 <  2:287461/148‑1 (MQ=255)
                             cacaAGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGc                                                                                                                                  <  4:68335/130‑1 (MQ=255)
                             cacaAGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGc                                                                                                                                  >  3:68335/1‑130 (MQ=255)
                              acaAGGATCTTTTTCAAAATGGATGAGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATTGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGaa                                                                                                              >  5:142106/1‑149 (MQ=255)
                               caAGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAAt                                                                                                             <  5:144756/149‑1 (MQ=255)
                                 aGGCTCTTTATCATAATGGATGTGCCTCTTGGGTAGGGCCTGACTCAGGCTCAACGGGGAAATGGGCGGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATgg                                                                                                           <  6:106833/149‑1 (MQ=255)
                                 aGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATgg                                                                                                           <  4:318727/149‑1 (MQ=255)
                                 aGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATgg                                                                                                           <  5:94866/149‑1 (MQ=255)
                                  ggATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATggg                                                                                                          <  6:250954/149‑1 (MQ=255)
                                  ggATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATggg                                                                                                          <  1:22942/149‑1 (MQ=255)
                                     tCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGta                                                                                                        >  7:347302/1‑148 (MQ=255)
                                      cTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGta                                                                                                        >  5:318504/1‑147 (MQ=255)
                                       tttttCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGta                                                                                                        >  5:221322/1‑146 (MQ=255)
                                         tttCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGtata                                                                                                      <  6:221322/146‑1 (MQ=255)
                                             aTAATGGATGTGCCGCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATTGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACaa                                                                                               <  6:142106/149‑1 (MQ=255)
                                              tAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAg                                                                                              >  8:323358/1‑149 (MQ=255)
                                              tAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAg                                                                                              >  1:91125/1‑149 (MQ=255)
                                               aaTGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATAAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAg                                                                                              >  5:176287/1‑148 (MQ=255)
                                                 tGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATTGCGACAAGGaa                                                                                           <  1:199832/149‑1 (MQ=255)
                                                   gATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAAt                                                                                          <  8:347302/148‑1 (MQ=255)
                                                     tgtgCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATaga                                                                                       <  8:52913/149‑1 (MQ=255)
                                                     tgtgCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATaga                                                                                       <  1:183372/149‑1 (MQ=255)
                                                       tgCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATagaga                                                                                     <  2:91125/149‑1 (MQ=255)
                                                       tgCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATagaga                                                                                     <  5:52452/149‑1 (MQ=255)
                                                       tgCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATagag                                                                                      <  5:102823/148‑1 (MQ=255)
                                                        gCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATaga                                                                                       >  7:265220/1‑146 (MQ=255)
                                                         ccTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAgagaga                                                                                   <  6:318504/149‑1 (MQ=255)
                                                                   aGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAgg                                                                         >  5:210000/1‑149 (MQ=255)
                                                                    gTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAgg                                                                         >  6:353682/1‑148 (MQ=255)
                                                                        cTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATa                                                                    >  2:10366/1‑149 (MQ=255)
                                                                          gAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATAc                                                                   >  7:136655/1‑148 (MQ=255)
                                                                           aaTCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCa                                                                 <  6:176287/149‑1 (MQ=255)
                                                                             tCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATc                                                               <  5:353682/149‑1 (MQ=255)
                                                                             tCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATc                                                               >  4:338438/1‑149 (MQ=255)
                                                                                     aaCGGGGAAATGTGCTGCTGCCATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATATCATCCTGCTCaa                                                       >  4:156903/1‑149 (MQ=255)
                                                                                     aaCGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCaa                                                       >  2:130683/1‑149 (MQ=255)
                                                                                     aaCGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCaa                                                       >  8:364489/1‑149 (MQ=255)
                                                                                         gggAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATa                                                                                         >  4:104885/1‑111 (MQ=255)
                                                                                          ggAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGt                                                                                                         >  8:1224/1‑94 (MQ=255)
                                                                                          ggAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGt                                                                                                         <  7:1224/94‑1 (MQ=255)
                                                                                           gAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGa                                                 >  8:83562/1‑149 (MQ=255)
                                                                                           gAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGCCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGa                                                 <  8:103221/149‑1 (MQ=255)
                                                                                           gAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGCCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGa                                                 <  6:31896/149‑1 (MQ=255)
                                                                                            aaaTGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGat                                                <  3:338438/149‑1 (MQ=255)
                                                                                              aTGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGata                                               >  8:340300/1‑148 (MQ=255)
                                                                                                 tgctgcTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAAc                                           >  7:24499/1‑149 (MQ=255)
                                                                                                    tgctgACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATaga                                                                                       <  2:91020/102‑1 (MQ=255)
                                                                                                    tgctgACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATaga                                                                                       >  1:91020/1‑102 (MQ=255)
                                                                                                         aCATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCt                                   <  8:136655/149‑1 (MQ=255)
                                                                                                           aTACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGt                                 >  7:154587/1‑149 (MQ=255)
                                                                                                           aTACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGt                                 <  8:10826/149‑1 (MQ=255)
                                                                                                                 gCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCAc                           >  7:126785/1‑149 (MQ=255)
                                                                                                                    agagCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGcc                        >  6:94442/1‑149 (MQ=255)
                                                                                                                          gTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGc                  >  4:329934/1‑149 (MQ=255)
                                                                                                                           tCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGcc                 >  7:144162/1‑149 (MQ=255)
                                                                                                                           tCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGcc                 >  2:26616/1‑149 (MQ=255)
                                                                                                                             aaTAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCAc               <  6:210000/149‑1 (MQ=255)
                                                                                                                             aaTAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCAc               >  3:74172/1‑149 (MQ=255)
                                                                                                                                  gCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGt                                 >  4:45128/1‑126 (MQ=255)
                                                                                                                                  gCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGt                                 <  3:45128/126‑1 (MQ=255)
                                                                                                                                  gCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTtctc          <  7:83562/149‑1 (MQ=255)
                                                                                                                                     gATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGGTATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCt       <  1:10366/149‑1 (MQ=255)
                                                                                                                                          aaTCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATc                                                               >  5:355313/1‑88 (MQ=255)
                                                                                                                                          aaTCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATc                                                               <  6:355313/88‑1 (MQ=255)
                                                                                                                                          aaTCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGat                                                >  4:101912/1‑103 (MQ=255)
                                                                                                                                          aaTCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCAcc  >  8:118535/1‑149 (MQ=255)
                                                                                                                                          aaTCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCAcc  >  6:107726/1‑149 (MQ=255)
                                                                                                                                          aaTCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCCGCAcc  >  8:215531/1‑149 (MQ=255)
                                                                                                                                          aaTCACTCAAATCACGCAGGCTGGCAAACTCCGCCGGAGAATGGGTATAGCGACAATGAATAGAGAGACTGGAGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCAcc  >  1:62968/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                              
GTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATGGATGTGCCTCTTGTGTAGTGCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACC  >  NC_000913/4530487‑4530773

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: