Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 958,773 C→G 5.6% intergenic (+32/‑39) serC → / → aroA 3‑phosphoserine/phosphohydroxythreonine aminotransferase/5‑enolpyruvylshikimate‑3‑phosphate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913958,7730CG 5.6% 552.2 / 17.4 159intergenic (+32/‑39)serC/aroA3‑phosphoserine/phosphohydroxythreonine aminotransferase/5‑enolpyruvylshikimate‑3‑phosphate synthetase
Reads supporting (aligned to +/- strand):  ref base C (72/78);  new base G (3/6);  total (75/84)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.08e-01

Only 100 of 164 total aligned reads displayed.

CTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACG  >  NC_000913/958625‑958920
                                                                                                                                                    |                                                                                                                                                   
cTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCAcagccagc                                                                                                                                                     >  1:249665/1‑149 (MQ=255)
 tGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACagccagcc                                                                                                                                                    >  2:889203/1‑149 (MQ=255)
     aGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCGATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACATACTTCATGGTTCAGTTCGAACGCCGTCACGTTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgtg                                                                                                                                                <  4:273737/149‑1 (MQ=255)
      ggTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTgg                                                                                                                                               <  3:817407/149‑1 (MQ=255)
      ggTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTgg                                                                                                                                               >  4:472311/1‑149 (MQ=255)
      ggTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTgg                                                                                                                                               <  1:613299/149‑1 (MQ=255)
        tCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTgggg                                                                                                                                             <  3:786126/149‑1 (MQ=255)
        tCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTgggg                                                                                                                                             <  3:456751/149‑1 (MQ=255)
         cACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGt                                                                                                                                            >  1:255140/1‑149 (MQ=255)
         cACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGt                                                                                                                                            >  3:334427/1‑149 (MQ=255)
            cGTGTGGTCGGCGGAACGCGCACTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAACGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGtttt                                                                                                                                         >  4:1077740/1‑149 (MQ=255)
             gtgtGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAAACCCCACAGGCTGGCTGTGGGGattaa                                                                                                                                        <  2:18240/149‑3 (MQ=255)
              tgtgGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCCGCCTGTGGGGTTTTTa                                                                                                                                       <  4:523229/149‑1 (MQ=255)
              tgtgGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTa                                                                                                                                       <  2:923567/149‑1 (MQ=255)
                 ggTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTAttt                                                                                                                                    <  1:881249/149‑1 (MQ=255)
                  gTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTa                                                                                                                                       <  2:472249/145‑1 (MQ=255)
                  gTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTa                                                                                                                                       >  1:472249/1‑145 (MQ=255)
                    cggcggAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtg                                                                                                                                 <  2:255140/149‑1 (MQ=255)
                    cggcggAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtg                                                                                                                                 <  1:607985/149‑1 (MQ=255)
                        ggAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtgttgt                                                                                                                             <  1:307607/149‑1 (MQ=255)
                           aTGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTaga                                                                                                                          <  4:36637/149‑1 (MQ=255)
                             gcgcgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTagaga                                                                                                                        <  3:554248/149‑1 (MQ=255)
                             gcgcgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTagaga                                                                                                                        <  3:288342/149‑1 (MQ=255)
                              cgcgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTagagag                                                                                                                       <  2:387143/149‑1 (MQ=255)
                                cgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagtt                                                                                                                     <  2:1087515/149‑1 (MQ=255)
                                 gcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAGATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttg                                                                                                                    >  4:878295/1‑149 (MQ=255)
                                 gcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttg                                                                                                                    >  3:48527/1‑149 (MQ=255)
                                  cTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttga                                                                                                                   <  3:1077740/149‑1 (MQ=255)
                                    tCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGtt                                                                                                                                           >  4:493715/1‑123 (MQ=255)
                                    tCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGtt                                                                                                                                           <  1:780201/123‑1 (MQ=255)
                                    tCTATTCATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGtt                                                                                                                                           >  2:780201/1‑123 (MQ=255)
                                      tATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTc                                                                                                               >  2:165173/1‑149 (MQ=255)
                                        tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCAt                                                                                                             <  4:815546/149‑1 (MQ=255)
                                        tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCAt                                                                                                             <  1:10186/149‑1 (MQ=255)
                                        tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCAt                                                                                                             >  2:256999/1‑149 (MQ=255)
                                            taACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGaa                                                                                                         <  1:1089236/149‑1 (MQ=255)
                                             aaCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAAt                                                                                                        >  1:80631/1‑149 (MQ=255)
                                                gCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTggggat                                                                                                                                           >  1:868534/1‑109 (MQ=255)
                                                gCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTggggat                                                                                                                                           <  2:868534/111‑3 (MQ=255)
                                                      ccGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGt                                                                                               >  2:698618/1‑149 (MQ=255)
                                                         cTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTagag                                                                                                                         >  2:32397/1‑120 (MQ=255)
                                                         cTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTAc                                                                                            <  4:334427/149‑1 (MQ=255)
                                                              aGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAAccc                                                                                       >  2:1029646/1‑149 (MQ=255)
                                                               ggCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCa                                                                                      >  4:390850/1‑149 (MQ=255)
                                                               ggCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCa                                                                                      >  3:87076/1‑149 (MQ=255)
                                                                  gTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGtt                                                                                                                                           >  4:47038/1‑93 (MQ=255)
                                                                   ttAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgt                                                                                                                                                 <  3:959916/86‑1 (MQ=255)
                                                                   ttAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgt                                                                                                                                                 >  4:959916/1‑86 (MQ=255)
                                                                            tGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGa                                                                         <  3:395743/149‑1 (MQ=255)
                                                                              aCAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgt                                                                                                                                                 <  3:533824/75‑1 (MQ=255)
                                                                              aCAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATg                                                                       >  4:10404/1‑149 (MQ=255)
                                                                                  aCTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCAc                                                                   <  1:939210/149‑1 (MQ=255)
                                                                                    ttCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTTGGGTTTTTATTTCTGTTGAAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACtc                                                                 >  1:787975/1‑148 (MQ=255)
                                                                                    ttCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAAc                                                                                         <  4:1054559/125‑1 (MQ=255)
                                                                                    ttCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAAc                                                                                         >  3:1054559/1‑125 (MQ=255)
                                                                                       aTGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTa                                                              <  2:824183/149‑1 (MQ=255)
                                                                                               gTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGccc                                                      >  4:325671/1‑149 (MQ=255)
                                                                                                ttCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCg                                                     >  1:685760/1‑149 (MQ=255)
                                                                                                ttCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCg                                                     >  2:802076/1‑149 (MQ=255)
                                                                                                   gAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGtt                                                  <  3:480530/149‑1 (MQ=255)
                                                                                                   gAACGCCGTCACGGTTAATGCCGAAACTTCGCTTAATCCCCACAGCCAGCCTGTGGGGTATTTATTTCTGTTGTAGAGAGATGAGTTCATGGAATCCCTGACGTTACAACCCAACGCTCGTGTCGATGGCACTATTAATCTGCCCGGtt                                                  <  4:670214/149‑1 (MQ=255)
                                                                                                    aaCGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTc                                                 >  3:432705/1‑149 (MQ=255)
                                                                                                           tCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGc                                          >  3:184243/1‑149 (MQ=255)
                                                                                                           tCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGc                                          >  1:117223/1‑149 (MQ=255)
                                                                                                             aCGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGt                                        <  3:220529/149‑1 (MQ=255)
                                                                                                                gTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTTTCGATGGCCCTATTAATCTGCCCGGTTCCAAGAGCGTTTc                                     >  1:475381/1‑149 (MQ=255)
                                                                                                                   aaTGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTaa                                  >  2:14124/1‑149 (MQ=255)
                                                                                                                      gCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCaa                                              <  1:444822/133‑1 (MQ=255)
                                                                                                                      gCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCaa                                              >  2:444822/1‑133 (MQ=255)
                                                                                                                        cGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACcgcg                             >  3:1060840/1‑149 (MQ=255)
                                                                                                                        cGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACcgcg                             >  1:359512/1‑149 (MQ=255)
                                                                                                                         gAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCgcgc                            <  4:682350/149‑1 (MQ=255)
                                                                                                                         gAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCgcgc                            >  1:633294/1‑149 (MQ=255)
                                                                                                                         gAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCgcgc                            <  4:87076/149‑1 (MQ=255)
                                                                                                                          aaaTTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCt                           <  1:913465/149‑1 (MQ=255)
                                                                                                                               ttGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGCCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTAtt                      <  1:361962/149‑1 (MQ=255)
                                                                                                                                   ttAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCaa                                              >  2:105862/1‑121 (MQ=255)
                                                                                                                                   ttAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTg                  >  4:453680/1‑149 (MQ=255)
                                                                                                                                   ttAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTg                  >  3:1050397/1‑149 (MQ=255)
                                                                                                                                   ttAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTg                  >  2:16328/1‑149 (MQ=255)
                                                                                                                                   ttAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTg                  >  2:368550/1‑149 (MQ=255)
                                                                                                                                   ttAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTg                  >  4:127384/1‑149 (MQ=255)
                                                                                                                                   ttAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGAACTCTTAATCTGCCCGGTTCCATGAGCGTTTCTAACCGCGCTTTATTGCTg                  >  3:678655/1‑149 (MQ=255)
                                                                                                                                     aaTCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGAt                                                                        >  3:614555/1‑93 (MQ=255)
                                                                                                                                     aaTCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGAt                                                                        <  4:614555/93‑1 (MQ=255)
                                                                                                                                     aaTCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTggc                <  1:707983/149‑1 (MQ=255)
                                                                                                                                      aTCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAAc                                 <  3:602995/131‑1 (MQ=255)
                                                                                                                                      aTCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAAc                                 >  4:602995/1‑131 (MQ=255)
                                                                                                                                       tCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTggcgg              >  3:88242/1‑149 (MQ=255)
                                                                                                                                       tCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTggcgg              >  1:23427/1‑149 (MQ=255)
                                                                                                                                        ccccACAGGCTGGCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTAtt                      >  3:544001/1‑140 (MQ=255)
                                                                                                                                        ccccACAGGCTGGCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTAtt                      <  4:544001/140‑1 (MQ=255)
                                                                                                                                        ccccACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGATAGTTGAGTGCATGGAATCCCTGACGTTACAACCCATCGCTCGGGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTggcggc             >  4:625559/1‑149 (MQ=255)
                                                                                                                                         cccACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCa            >  1:1029988/1‑149 (MQ=255)
                                                                                                                                          ccACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCa                                                                    >  2:740486/1‑92 (MQ=255)
                                                                                                                                          ccACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCAt           <  2:475381/149‑1 (MQ=255)
                                                                                                                                          ccACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCAt           <  4:415677/149‑1 (MQ=255)
                                                                                                                                           cacaGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCAtt          <  2:633294/149‑1 (MQ=255)
                                                                                                                                            acaGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTa         <  1:750432/149‑1 (MQ=255)
                                                                                                                                                  agccTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCacac   <  3:932504/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCATTAGCACACG  >  NC_000913/958625‑958920

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