Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,274,920 T→C 5.3% R379R (CGT→CGC ychO → putative invasin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,274,9200TC 5.3% 226.8 / 3.2 75R379R (CGT→CGCychOputative invasin
Reads supporting (aligned to +/- strand):  ref base T (34/37);  new base C (2/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

GTAGTCGCTATGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTC  >  NC_000913/1274775‑1275066
                                                                                                                                                 |                                                                                                                                                  
gTAGTCGCTATGACAATCCGCAGCGAAATTCGCTAGCGACTCTTGAGTACCGAGAGCGTAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAAACTGGCGAAACAGTGCCGCTTAAATTACAAATCCGCATTCGTTa                                                                                                                                                  <  6:216689/148‑1 (MQ=255)
  aGTCGCTATGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACgg                                                                                                                                               >  2:197584/1‑149 (MQ=255)
   gTCGCTATGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAATAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAACCCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACgg                                                                                                                                               <  3:101146/148‑1 (MQ=255)
   gTCGCTATGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACgg                                                                                                                                               <  5:362785/148‑1 (MQ=255)
        tATGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGGTAACGGTGTTTCTGGCGCCACCGCCGTGGGATCTAAAAACTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCg                                                                                                                                         >  8:392555/1‑149 (MQ=255)
         aTGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCgg                                                                                                                                        >  7:258618/1‑149 (MQ=255)
         aTGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCgg                                                                                                                                        >  4:88553/1‑149 (MQ=255)
         aTGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCCCAGTCGTTACGGTATTCgg                                                                                                                                        >  4:83632/1‑149 (MQ=255)
          tGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGc                                                                                                                                       <  4:86/149‑1 (MQ=255)
          tGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAAATACAAATCCGCAGTCGTTACGGTATTCgg                                                                                                                                        >  2:301921/1‑148 (MQ=255)
              aaTCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACt                                                                                                                                   >  8:264927/1‑149 (MQ=255)
                tCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGa                                                                                                                                 <  3:82071/149‑1 (MQ=255)
                   gCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGAttt                                                                                                                              <  7:175014/149‑1 (MQ=255)
                   gCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGAttt                                                                                                                              >  4:345981/1‑149 (MQ=255)
                     aGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTTGTATCTAAAACCTGGCGAAACAGTGCCGCTGATATTACAAATAAGCAGTCGTTACTGTATTCGGCAACTGATTTgg                                                                                                                            <  2:151167/149‑1 (MQ=255)
                     aGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTgg                                                                                                                            <  8:258618/149‑1 (MQ=255)
                     aGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTgg                                                                                                                            <  1:197584/149‑1 (MQ=255)
                       cGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGGGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCa                                                                                                                          <  4:247705/149‑1 (MQ=255)
                         aaaTAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAgg                                                                                                                        >  5:321321/1‑149 (MQ=255)
                              atCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGa                                                                                                                    <  7:335494/148‑1 (MQ=255)
                              atCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGa                                                                                                                    >  3:31976/1‑148 (MQ=255)
                              atCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGAt                                                                                                                   <  6:5651/149‑1 (MQ=255)
                                    agACTATTGAGTACCGACAGCGAAAAACGTTAACGGTTTTTCTGGCGAAACCGACGTGTGATATAAAACCTGGCGAAACAGTGCCTCTGAAATTACAAATACGCAGTCGTTACGGTCTTCGGCAACTGATTTGGCAGGGCGATACGCAg                                                                                                             <  6:273464/148‑1 (MQ=255)
                                      aCTCTTGAGTACCGCCAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGGGGGATCTAAAACCTGGCGAAACCGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGa                                                                                                            >  6:122479/1‑148 (MQ=255)
                                            gAGTACCGACAGCGAAAAACGTTAACGGTGTATCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAg                                                                                                     >  6:33858/1‑149 (MQ=255)
                                                 ccGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGa                                                                                                <  8:23760/149‑1 (MQ=255)
                                                  cGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGCTACGGTATTCgg                                                                                                                                        >  3:291680/1‑108 (MQ=255)
                                                  cGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGCTACGGTATTCgg                                                                                                                                        >  3:296685/1‑108 (MQ=255)
                                                  cGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGCTACGGTATTCgg                                                                                                                                        <  4:296685/108‑1 (MQ=255)
                                                  cGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGCTACGGTATTCgg                                                                                                                                        <  4:291680/108‑1 (MQ=255)
                                                       gCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCgg                                                                                          <  8:342885/149‑1 (MQ=255)
                                                        cGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATTCGCAGATATTAAGTTTGACGCCggg                                                                                         >  7:394312/1‑149 (MQ=255)
                                                        cGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCggg                                                                                         >  8:293678/1‑149 (MQ=255)
                                                         gAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGgc                                                                                        <  8:244552/149‑1 (MQ=255)
                                                           aaaaCGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATCCGGCAACTGAGTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGgcgc                                                                                      >  6:72358/1‑149 (MQ=255)
                                                            aaaCGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGt                                                                                                    >  1:175576/1‑134 (MQ=255)
                                                            aaaCGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGt                                                                                                    <  2:175576/134‑1 (MQ=255)
                                                              aCGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGcac                                                                                    >  6:187518/1‑148 (MQ=255)
                                                               cGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACaa                                                                                  <  4:31976/149‑1 (MQ=255)
                                                                 ttAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGAGATTAAGTTTGACGCCGGGCGCACAAGc                                                                                >  3:5929/1‑149 (MQ=255)
                                                                    aCGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCa                                                                              <  7:89107/148‑1 (MQ=255)
                                                                      ggTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAGACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCaa                                                                             >  7:143006/1‑147 (MQ=255)
                                                                      ggTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTATGGCATTc                                                                                                                                          <  7:102163/86‑1 (MQ=255)
                                                                      ggTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTAGGGAATTc                                                                                                                                          >  8:102163/1‑86 (MQ=255)
                                                                       gtgtTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAg                                                                          >  1:9843/1‑149 (MQ=255)
                                                                       gtgtTTCTGGCGACACCGCCGGGGGATCTAAAACCTGGCGAAACAGGGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAg                                                                          <  5:309806/149‑1 (MQ=255)
                                                                              tGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAg                                                                   >  8:108907/1‑149 (MQ=255)
                                                                               ggCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAgg                                                                  >  4:251387/1‑149 (MQ=255)
                                                                                gCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGa                                                                    <  4:115343/146‑1 (MQ=255)
                                                                                gCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGa                                                                    <  4:5929/146‑1 (MQ=255)
                                                                                gCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAggg                                                                 >  3:326877/1‑149 (MQ=255)
                                                                                  gACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCt                                                               <  7:293678/149‑1 (MQ=255)
                                                                                   acacCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTg                                                              >  8:84044/1‑149 (MQ=255)
                                                                                           gTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGa                                                      <  3:88553/149‑1 (MQ=255)
                                                                                           gTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGa                                                      <  3:83632/149‑1 (MQ=255)
                                                                                              ggATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATc                                                    <  1:301921/148‑1 (MQ=255)
                                                                                               gATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCAt                                                  >  4:356674/1‑149 (MQ=255)
                                                                                                aTCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCATCTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATg                                                 >  3:316733/1‑149 (MQ=255)
                                                                                                      aaCCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGAc                                           <  7:264927/149‑1 (MQ=255)
                                                                                                        ccTGGCGAGACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTg                                         <  8:143006/149‑1 (MQ=255)
                                                                                                         cTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTgg                                        <  3:345981/149‑1 (MQ=255)
                                                                                                         cTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTgg                                        >  7:369583/1‑149 (MQ=255)
                                                                                                         cTGGCGAAACAGTGCCGCGGAAATTACAAATCCGCAGTCGTTACGGTATTCCGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTgg                                        <  6:283170/149‑1 (MQ=255)
                                                                                                           ggCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTgg                                        >  8:58210/1‑147 (MQ=255)
                                                                                                              gAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGaa                                   <  7:58210/149‑1 (MQ=255)
                                                                                                               aaaCAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGAGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGAGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAAc                                  >  1:84457/1‑149 (MQ=255)
                                                                                                               aaaCAGTGCCGCTGAAATTACAAATCCGCAGGCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAAc                                  <  7:108907/149‑1 (MQ=255)
                                                                                                                     tGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGaa                            <  5:187518/149‑1 (MQ=255)
                                                                                                                       ccGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAg                           <  2:9843/148‑1 (MQ=255)
                                                                                                                        cGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAggg                         <  5:33858/149‑1 (MQ=255)
                                                                                                                                ttACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGgcaagcaa                 >  3:127832/1‑149 (MQ=255)
                                                                                                                                 tACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAAt                >  5:25794/1‑149 (MQ=255)
                                                                                                                                          gCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGCACGGGGAAGGGGCAAGCAATCACTGGCGa       >  6:269851/1‑149 (MQ=255)
                                                                                                                                          gCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGa       >  2:136958/1‑149 (MQ=255)
                                                                                                                                                gTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTc  <  6:321321/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GTAGTCGCTATGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTC  >  NC_000913/1274775‑1275066

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: