Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,271,825 T→A 8.1% intergenic (‑55/+42) garK ← / ← garR glycerate kinase I/tartronate semialdehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,271,8250TA 8.1% 241.9 / 11.2 74intergenic (‑55/+42)garK/garRglycerate kinase I/tartronate semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base A (4/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-01

AATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCATCATCTCCATAA  >  NC_000913/3271681‑3271973
                                                                                                                                                |                                                                                                                                                    
aaTTTCCCGAAATCCTTTTTCTATCGCCTGCGTAACCTCGCTGCCCGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGTTTTTCATACCTAGGCCTGTTACCACATAACGCCGGAGGGCGTGTCTCTTATTCGGCCTGGATTCCa                                                                                                                                                  <  2:751931/149‑1 (MQ=255)
aaTTTCCCGAAATCCTGTTTCTATCGGCTGCGCAACCTCGCTGGCAGATAAACTTTCTCTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATTACTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCa                                                                                                                                                  <  2:613304/149‑1 (MQ=255)
 aTTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAg                                                                                                                                                 <  3:562533/149‑1 (MQ=255)
       cGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCgg                                                                                                                                           >  1:669485/1‑149 (MQ=255)
        gAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGa                                                                                                                                          >  1:99940/1‑149 (MQ=255)
            tCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGCGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTGGGCCTGGATTCCAGGCCCGGATTGc                                                                                                                                      >  1:596068/1‑149 (MQ=255)
              cTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCaa                                                                                                                                    <  2:601118/149‑1 (MQ=255)
               tttttCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAAt                                                                                                                                   >  4:191397/1‑149 (MQ=255)
                         gCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCg                                                                                                                         >  1:458694/1‑149 (MQ=255)
                          ccTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCgg                                                                                                                        >  3:148192/1‑149 (MQ=255)
                             gcgcAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCa                                                                                                                     >  1:405159/1‑149 (MQ=255)
                               gcAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCaa                                                                                                                                    >  4:514765/1‑132 (MQ=255)
                               gcAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCaa                                                                                                                                    <  3:514765/132‑1 (MQ=255)
                                 aaCCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCacgac                                                                                                                 >  3:433763/1‑149 (MQ=255)
                                   ccTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGt                                                                                                               >  1:623263/1‑149 (MQ=255)
                                    cTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTc                                                                                                              <  2:173279/149‑1 (MQ=255)
                                          ggCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAAc                                                                                                        <  3:191397/149‑1 (MQ=255)
                                            cAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGa                                                                                                      <  3:129836/149‑1 (MQ=255)
                                            cAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGa                                                                                                      <  4:190299/149‑1 (MQ=255)
                                            cAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGa                                                                                                      <  2:638077/149‑1 (MQ=255)
                                            cAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGa                                                                                                      <  3:726512/149‑1 (MQ=255)
                                             aGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCAGCCTGTATTCCAGGCCCGGATTGCAATACGCCATCCGGGCacga                                                                                                                  <  4:625750/136‑1 (MQ=255)
                                             aGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCAGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCacga                                                                                                                  >  3:625750/1‑136 (MQ=255)
                                              gATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGt                                                                                                    <  3:97242/149‑1 (MQ=255)
                                              gATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGt                                                                                                    >  3:17642/1‑149 (MQ=255)
                                                tAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGccc                                                                                                                                             <  1:321674/106‑1 (MQ=255)
                                                tAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGccc                                                                                                                                             >  2:321674/1‑106 (MQ=255)
                                                   aCTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGAATCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACtt                                                                                               >  3:707641/1‑149 (MQ=255)
                                                    ctttctttATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTc                                                                                              >  1:20441/1‑149 (MQ=255)
                                                       tctttATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGAc                                                                                           <  1:568157/149‑1 (MQ=255)
                                                         tttATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGACCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACtt                                                                                         >  1:572476/1‑149 (MQ=255)
                                                           tataAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTc                                                                                       <  1:4213/149‑1 (MQ=255)
                                                             taAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGc                                                                                     >  4:94159/1‑149 (MQ=255)
                                                              aaGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGcc                                                                                    >  1:384467/1‑149 (MQ=255)
                                                               agagTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCa                                                                                   >  2:287743/1‑149 (MQ=255)
                                                                 agTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGt                                                                                 >  1:88009/1‑149 (MQ=255)
                                                                        ggCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCgt                                                                          <  2:20441/149‑1 (MQ=255)
                                                                              tACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCTGCCTGGATTCCAGGCCCGGATTGGAATGCGCCATCCGGGCCCGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGc                                                                    >  4:519960/1‑149 (MQ=255)
                                                                                  aTTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATa                                                                                                                                  <  1:140511/83‑1 (MQ=255)
                                                                                  aTTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATa                                                                                                                                  >  2:140511/1‑83 (MQ=255)
                                                                                   ttttCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCa                                                                                   <  3:115727/129‑1 (MQ=255)
                                                                                   ttttCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCa                                                                                   >  4:115727/1‑129 (MQ=255)
                                                                                            cTATGCCTGTTACCACGTGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTa                                                      >  4:110093/1‑149 (MQ=255)
                                                                                                 ccTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATc                                                 <  4:148192/149‑1 (MQ=255)
                                                                                                  cTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATcc                                                >  1:660874/1‑149 (MQ=255)
                                                                                                   tGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATccg                                               >  2:741854/1‑149 (MQ=255)
                                                                                                    gTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATccgc                                              <  3:441863/149‑1 (MQ=255)
                                                                                                       aCCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAg                                                             >  3:333577/1‑131 (MQ=255)
                                                                                                       aCCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAg                                                             <  4:333577/131‑1 (MQ=255)
                                                                                                         cacaTGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTc                                         >  3:584460/1‑149 (MQ=255)
                                                                                                              gACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTaaa                                    <  2:95256/149‑1 (MQ=255)
                                                                                                              gACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTaaa                                    <  1:655352/149‑1 (MQ=255)
                                                                                                                   cGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATc                               <  4:633454/149‑1 (MQ=255)
                                                                                                                   cGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATc                               <  1:607757/149‑1 (MQ=255)
                                                                                                                   cGGAGGGCGTTTCTCTTATTCGGCCTAGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATc                               <  4:96348/149‑1 (MQ=255)
                                                                                                                    ggAGGGCGTTTCTCTTATGCGGCCTGGATTCCACGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCCCGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCt                              <  3:150914/149‑1 (MQ=255)
                                                                                                                      aGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGc                                                                >  3:284127/1‑113 (MQ=255)
                                                                                                                      aGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGc                                                                <  4:284127/113‑1 (MQ=255)
                                                                                                                      aGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGc                            >  2:683588/1‑149 (MQ=255)
                                                                                                                        ggCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTc                          <  1:683588/149‑1 (MQ=255)
                                                                                                                        ggCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTc                          >  1:703353/1‑149 (MQ=255)
                                                                                                                                tctTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGcc                  >  2:467717/1‑149 (MQ=255)
                                                                                                                                 ctTATTCGGCCTGGAATCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCt                 <  4:707641/149‑1 (MQ=255)
                                                                                                                                  ttATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTg                >  4:318151/1‑149 (MQ=255)
                                                                                                                                    tccgggCCTGGAATCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTg                >  4:36070/5‑147 (MQ=255)
                                                                                                                                    tccgggCCTGGAATCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTg                <  3:36070/143‑1 (MQ=255)
                                                                                                                                    tccgggCCTGGAATCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGca              >  3:701755/5‑149 (MQ=255)
                                                                                                                                       gggCCTGGAATCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGcatca           >  1:27258/2‑149 (MQ=255)
                                                                                                                                         gCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCatcatc         <  1:734610/149‑1 (MQ=255)
                                                                                                                                         gCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCatcatc         >  3:394849/1‑149 (MQ=255)
                                                                                                                                         gCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCatcatc         <  2:623263/149‑1 (MQ=255)
                                                                                                                                         gCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCatcatc         <  2:669485/149‑1 (MQ=255)
                                                                                                                                           cTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCATCAtctc       <  1:287743/149‑1 (MQ=255)
                                                                                                                                           cTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCATCAcctc       <  4:205625/149‑1 (MQ=255)
                                                                                                                                                ttCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCATCATCTCCATaa  >  2:618553/1‑149 (MQ=255)
                                                                                                                                                ttCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCATCATCTCCATaa  >  1:515824/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCATCATCTCCATAA  >  NC_000913/3271681‑3271973

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: