Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,703,593 A→G 5.8% I239T (ATC→ACC)  ydgJ ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,703,5930AG 5.8% 516.4 / 14.8 154I239T (ATC→ACC) ydgJputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (75/70);  new base G (5/4);  total (80/74)
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.34e-01

Only 100 of 165 total aligned reads displayed.

GACACGTTCCTCACCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACC  >  NC_000913/1703445‑1703741
                                                                                                                                                    |                                                                                                                                                    
gACACGTTCCTCACCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGa                                                                                                                                                      >  3:478955/1‑149 (MQ=255)
 acacGTTCCTCACCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGat                                                                                                                                                     >  2:21505/1‑149 (MQ=255)
 acacGTTCCTCACCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGat                                                                                                                                                     >  4:1169179/1‑149 (MQ=255)
     gTTCCTCACCTTCCACGCGGGTAAGTACGCCATCACGCATCTCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATACACGATATAc                                                                                                                                                 <  4:1230110/149‑1 (MQ=255)
        ccTCACCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGGTATATCGt                                                                                                                                              >  3:567596/1‑149 (MQ=255)
            aCCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGGTATATCGTGcat                                                                                                                                          <  4:567596/149‑3 (MQ=255)
            aCCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTg                                                                                                                                          <  2:506593/149‑1 (MQ=255)
              cTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGATGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGa                                                                                                                                                      >  4:538114/1‑135 (MQ=255)
              cTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGa                                                                                                                                                      >  2:1315843/1‑135 (MQ=255)
              cTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGa                                                                                                                                                      <  1:1315843/135‑1 (MQ=255)
               ttCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGGTATATCGTGcatgga                                                                                                                                       <  1:319394/149‑6 (MQ=255)
                   aCGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGc                                                                                                                                   >  2:306745/1‑149 (MQ=255)
                         gTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGc                                                                                                                             >  2:1879951/1‑149 (MQ=255)
                         gTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGc                                                                                                                             <  3:1660970/149‑1 (MQ=255)
                           aaGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCa                                                                                                                           <  1:34310/149‑1 (MQ=255)
                                   cATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGa                                                                                                                                         <  4:1075595/127‑1 (MQ=255)
                                   cATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGa                                                                                                                                         >  3:1075595/1‑127 (MQ=255)
                                   cATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTa                                                                                                                   >  1:1540157/1‑149 (MQ=255)
                                       aCGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGata                                                                                                                                                    >  3:56566/1‑112 (MQ=255)
                                        cGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTa                                                                                                                   >  2:777245/1‑144 (MQ=255)
                                         gCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACTTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGa                                                                                                                                         <  2:1402120/121‑1 (MQ=255)
                                         gCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCTGGATCCATGCACGATATATCGTGCTGa                                                                                                                                         >  1:1402120/1‑121 (MQ=255)
                                         gCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCACGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTa                                                                                                            <  2:1365344/149‑1 (MQ=255)
                                               cGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTATTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATg                                                                                                       >  1:287588/1‑149 (MQ=255)
                                                     cccAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTcgc                                                                                                 <  1:21505/149‑1 (MQ=255)
                                                      ccAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTcgcc                                                                                                >  1:590867/1‑149 (MQ=255)
                                                      ccAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTcgcc                                                                                                >  2:1437787/1‑149 (MQ=255)
                                                        aGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTcgccgc                                                                                              >  1:1207960/1‑149 (MQ=255)
                                                          tCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCtg                                                                                            >  1:1308159/1‑149 (MQ=255)
                                                            ttCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATAcc                                                                                                                                                >  1:763555/1‑95 (MQ=255)
                                                                   ggTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGa                                                                                   <  4:1844319/149‑1 (MQ=255)
                                                                    gTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTc                                                                                                   >  3:1690543/1‑132 (MQ=255)
                                                                     tAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGa                                                                                                      >  1:1245625/1‑128 (MQ=255)
                                                                      aGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACaa                                                                                <  1:1437787/149‑1 (MQ=255)
                                                                          gCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGa                                                                                                      >  1:1036482/1‑123 (MQ=255)
                                                                            tCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGa                                                                              >  2:1731200/1‑145 (MQ=255)
                                                                            tCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGa                                                                              <  1:1731200/145‑1 (MQ=255)
                                                                                  tttttCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATg                                                                            >  4:1466397/1‑141 (MQ=255)
                                                                                  tttttCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATg                                                                            <  3:1466397/141‑1 (MQ=255)
                                                                                         gACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCg                                                             >  4:416414/1‑149 (MQ=255)
                                                                                         gACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCg                                                             >  2:492531/1‑149 (MQ=255)
                                                                                                ttcCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACt                                                      >  2:381184/1‑149 (MQ=255)
                                                                                                 tcCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTg                                                     >  1:211613/1‑149 (MQ=255)
                                                                                                   cTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATCCCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGCCAAGCTGGCGTGGCAATAATCGGTCGACTgcg                                                   >  2:1205067/1‑149 (MQ=255)
                                                                                                    tgtgGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTgcgc                                                  >  2:737058/1‑149 (MQ=255)
                                                                                                      tgGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTc                                                >  3:306437/1‑149 (MQ=255)
                                                                                                       ggATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTcc                                               <  3:416414/149‑1 (MQ=255)
                                                                                                       ggATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCACGGTACCGTGTAAAATGACTCGCCTCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTcc                                               <  4:1376945/149‑1 (MQ=255)
                                                                                                        gATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCg                                              >  4:718241/1‑149 (MQ=255)
                                                                                                         aTCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCgg                                             >  1:54253/1‑149 (MQ=255)
                                                                                                         aTCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCgg                                             >  3:291004/1‑149 (MQ=255)
                                                                                                         aTCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCgg                                             >  3:1186574/1‑149 (MQ=255)
                                                                                                              ggCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGt                                        >  3:94598/1‑149 (MQ=255)
                                                                                                              ggCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGt                                        <  1:737058/149‑1 (MQ=255)
                                                                                                              ggCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGt                                        >  2:1334190/1‑149 (MQ=255)
                                                                                                                   tatTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGc                                                                          >  1:463999/1‑110 (MQ=255)
                                                                                                                   tatTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGc                                                                          <  4:493747/110‑1 (MQ=255)
                                                                                                                   tatTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGc                                                                          >  3:493747/1‑110 (MQ=255)
                                                                                                                    atTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtgt                                  <  4:291004/149‑1 (MQ=255)
                                                                                                                     tttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTc                                                          <  2:1031626/124‑1 (MQ=255)
                                                                                                                     tttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTc                                                          >  1:1031626/1‑124 (MQ=255)
                                                                                                                     tttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtgtg                                 >  3:1208643/1‑149 (MQ=255)
                                                                                                                     tttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtgtg                                 >  4:1377517/1‑149 (MQ=255)
                                                                                                                     tttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtgtg                                 >  3:1564005/1‑149 (MQ=255)
                                                                                                                     tttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtgtg                                 >  3:1087580/1‑149 (MQ=255)
                                                                                                                      ttCATATACCTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGGGc                                <  1:1205067/149‑1 (MQ=255)
                                                                                                                         acaTGACTGCCGCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCaa                             >  2:222345/1‑149 (MQ=255)
                                                                                                                         acaTAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATa                                                                                      >  4:726162/1‑92 (MQ=255)
                                                                                                                          caTAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCaaa                            <  2:318145/149‑1 (MQ=255)
                                                                                                                           aTAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAAc                                     <  2:166729/139‑1 (MQ=255)
                                                                                                                           aTAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAAc                                     >  1:166729/1‑139 (MQ=255)
                                                                                                                           aTAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAAt                           >  3:1605886/1‑149 (MQ=255)
                                                                                                                           aTAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAAt                           <  4:1605886/149‑1 (MQ=255)
                                                                                                                             aaCTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGt                                                           >  3:1680893/1‑115 (MQ=255)
                                                                                                                             aaCTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGt                                                           <  4:1680893/115‑1 (MQ=255)
                                                                                                                              aCTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGcc                               >  4:358024/1‑142 (MQ=255)
                                                                                                                               cTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTAc                       <  3:874440/149‑1 (MQ=255)
                                                                                                                                tGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTAcc                      <  2:1039999/149‑1 (MQ=255)
                                                                                                                                   cTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGa                                                                              <  2:1238933/90‑1 (MQ=255)
                                                                                                                                   cTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTc                   >  3:1162599/1‑149 (MQ=255)
                                                                                                                                     cGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCAt                 >  2:1416513/1‑149 (MQ=255)
                                                                                                                                      gggATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATg                >  4:1852041/1‑149 (MQ=255)
                                                                                                                                       ggATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATg                                                                            <  3:475067/88‑1 (MQ=255)
                                                                                                                                       ggATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATg                                                                            >  4:475067/1‑88 (MQ=255)
                                                                                                                                         gagtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGTATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTg             >  4:1100649/7‑149 (MQ=255)
                                                                                                                                           gtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGAc           >  2:211395/5‑149 (MQ=255)
                                                                                                                                            cATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGAcc          <  4:1037989/149‑1 (MQ=255)
                                                                                                                                             aTGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtgt                                  >  1:1578508/1‑124 (MQ=255)
                                                                                                                                             aTGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtgt                                  <  2:1578508/124‑1 (MQ=255)
                                                                                                                                             aTGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCg         >  3:828860/1‑149 (MQ=255)
                                                                                                                                              tGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtg                                   >  4:1577046/1‑122 (MQ=255)
                                                                                                                                              agCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCgg        >  2:1854349/2‑149 (MQ=255)
                                                                                                                                               gCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGt       <  1:1879951/149‑1 (MQ=255)
                                                                                                                                               gCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGt       >  2:995498/1‑149 (MQ=255)
                                                                                                                                                cACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAAGAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTa      <  3:1100649/149‑1 (MQ=255)
                                                                                                                                                cACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTa      >  3:1392651/1‑149 (MQ=255)
                                                                                                                                                   gATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAAc   <  4:932/149‑1 (MQ=255)
                                                                                                                                                   gATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAAc   <  3:1852041/149‑1 (MQ=255)
                                                                                                                                                   gATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAAc   <  3:1521713/149‑1 (MQ=255)
                                                                                                                                                    atatACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAAcc  <  3:718241/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GACACGTTCCTCACCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACC  >  NC_000913/1703445‑1703741

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: