Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 958,769 C→G 13.3% intergenic (+28/‑43) 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,7690CG13.3% 408.8 / 49.2 128intergenic (+28/‑43)serC/aroA3‑phosphoserine/phosphohydroxythreonine aminotransferase/5‑enolpyruvylshikimate‑3‑phosphate synthetase
Reads supporting (aligned to +/- strand):  ref base C (52/59);  new base G (8/9);  total (60/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

Only 100 of 135 total aligned reads displayed.

TGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCAT  >  NC_000913/958621‑958911
                                                                                                                                                    |                                                                                                                                              
tGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCAcagc                                                                                                                                                >  4:266126/1‑149 (MQ=255)
  cACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCAcagcca                                                                                                                                              <  1:381967/149‑1 (MQ=255)
   aCTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCAcagccag                                                                                                                                             >  2:1445532/1‑149 (MQ=255)
     tGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACagccagcc                                                                                                                                           >  3:874395/1‑149 (MQ=255)
        aaGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgt                                                                                                                                        >  1:869158/1‑149 (MQ=255)
            tCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTgggg                                                                                                                                    <  1:1445532/149‑1 (MQ=255)
            tCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTgggg                                                                                                                                    >  4:55570/1‑149 (MQ=255)
               ccGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTggg                                                                                                                                     >  2:525522/1‑145 (MQ=255)
               ccGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTggg                                                                                                                                     <  1:525522/145‑1 (MQ=255)
                  tgtgGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTa                                                                                                                              <  1:1053398/149‑1 (MQ=255)
                   gtgGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTAt                                                                                                                             <  2:869158/149‑1 (MQ=255)
                   gtgGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTAt                                                                                                                             <  3:55570/149‑1 (MQ=255)
                    tgGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTggg                                                                                                                                     >  3:1606968/1‑140 (MQ=255)
                    tgGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTggg                                                                                                                                     <  4:1606968/140‑1 (MQ=255)
                     ggTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTAttt                                                                                                                           <  4:874395/149‑1 (MQ=255)
                     ggTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTAttt                                                                                                                           <  3:534879/149‑1 (MQ=255)
                              aaTGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTag                                                                                                                  <  2:218621/149‑1 (MQ=255)
                                 gcgcgcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTagaga                                                                                                               <  3:266126/149‑1 (MQ=255)
                                     gcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTagaga                                                                                                               >  4:699791/1‑145 (MQ=255)
                                     gcTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTagaga                                                                                                               <  3:699791/145‑1 (MQ=255)
                                       ttCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTa                                                                                                                              >  2:1823220/1‑128 (MQ=255)
                                       ttCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTa                                                                                                                              <  1:1823220/128‑1 (MQ=255)
                                       ttCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTagaga                                                                                                               >  2:753813/1‑143 (MQ=255)
                                            tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTgggga                                                                                                                                   >  1:1654637/1‑117 (MQ=255)
                                            tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttga                                                                                                          >  2:1809009/1‑143 (MQ=255)
                                            tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttga                                                                                                          <  1:1809009/143‑1 (MQ=255)
                                            tttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCAt                                                                                                    <  4:422771/149‑1 (MQ=255)
                                             ttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTggg                                                                                                                                     <  2:751949/115‑1 (MQ=255)
                                             ttATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTggg                                                                                                                                     >  1:751949/1‑115 (MQ=255)
                                                taACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtgtt                                                                                                                      <  1:757197/127‑1 (MQ=255)
                                                      cATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGtttt                                                                                                                                >  2:17962/1‑111 (MQ=255)
                                                       aTGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGa                                                                                         <  1:1184095/149‑1 (MQ=255)
                                                                  aGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTaga                                                                                                                 >  4:888571/1‑114 (MQ=255)
                                                                  aGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTaga                                                                                                                 <  3:888571/114‑1 (MQ=255)
                                                                    gCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttga                                                                                                          <  4:927442/119‑1 (MQ=255)
                                                                    gCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttga                                                                                                          >  3:927442/1‑119 (MQ=255)
                                                                     cGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTggg                                                                                                                                     >  4:831568/1‑91 (MQ=255)
                                                                     cGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTggg                                                                                                                                     <  3:831568/91‑1 (MQ=255)
                                                                       ttAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCtgtg                                                                                                                                       >  1:284475/1‑87 (MQ=255)
                                                                          aaGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATc                                                                                              <  4:412144/125‑1 (MQ=255)
                                                                          aaGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATc                                                                                              >  3:412144/1‑125 (MQ=255)
                                                                            gcgcTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtg                                                                                                                        <  1:385595/97‑1 (MQ=255)
                                                                                tGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATc                                                                                              >  2:533537/1‑119 (MQ=255)
                                                                                tGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATc                                                                                              <  1:533537/119‑1 (MQ=255)
                                                                                  aCAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATg                                                                                                   <  2:228996/112‑1 (MQ=255)
                                                                                  aCAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATg                                                                                                   >  1:228996/1‑112 (MQ=255)
                                                                                    aGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttga                                                                                                          <  3:1554982/103‑1 (MQ=255)
                                                                                    aGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttga                                                                                                          >  4:1554982/1‑103 (MQ=255)
                                                                                       cTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtgttgt                                                                                                                    >  1:1851953/1‑90 (MQ=255)
                                                                                       cTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtgttgt                                                                                                                    <  2:1851953/90‑1 (MQ=255)
                                                                                        ttCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTgg                                                                                                                                      <  4:296499/71‑1 (MQ=255)
                                                                                        ttCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTgg                                                                                                                                      >  3:296499/1‑71 (MQ=255)
                                                                                        ttCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtgtt                                                                                                                      <  2:1507274/87‑1 (MQ=255)
                                                                                        ttCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtgtt                                                                                                                      >  1:1507274/1‑87 (MQ=255)
                                                                                        ttCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagtt                                                                                                            <  4:90689/97‑1 (MQ=255)
                                                                                        ttCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagtt                                                                                                            >  3:90689/1‑97 (MQ=255)
                                                                                        ttCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGa                                                                                         <  3:1284634/116‑1 (MQ=255)
                                                                                        ttCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTa                                                        >  3:1232629/1‑149 (MQ=255)
                                                                                          cATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTAtt                                                      >  1:1711564/1‑149 (MQ=255)
                                                                                             ggTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTaga                                                                                                                 >  3:1119959/1‑87 (MQ=255)
                                                                                             ggTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAAt                                                   <  2:1711564/149‑1 (MQ=255)
                                                                                                tGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACaa                                                                                 <  3:396925/116‑1 (MQ=255)
                                                                                                tGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACaa                                                                                 >  4:396925/1‑116 (MQ=255)
                                                                                                   gTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATgg                                                             >  3:395982/1‑133 (MQ=255)
                                                                                                   gTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATgg                                                             <  4:395982/133‑1 (MQ=255)
                                                                                                   gTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATc                                                  <  2:377164/144‑1 (MQ=255)
                                                                                                   gTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATc                                                  >  1:377164/1‑144 (MQ=255)
                                                                                                    ttCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGa                                                                <  2:361104/129‑1 (MQ=255)
                                                                                                    ttCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCg                                            >  3:882164/1‑149 (MQ=255)
                                                                                                         aCGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagtt                                                                                                            >  1:1019279/1‑80 (MQ=255)
                                                                                                         aCGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagtt                                                                                                            <  2:1019279/80‑1 (MQ=255)
                                                                                                                        aTGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAAc                        <  3:1594937/149‑1 (MQ=255)
                                                                                                                         tGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAAcc                       <  4:882164/149‑1 (MQ=255)
                                                                                                                           ccGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGtt                                                                                                                                  <  4:1280929/40‑1 (MQ=255)
                                                                                                                           ccGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGtt                                                                                                                                  >  3:1280929/1‑40 (MQ=255)
                                                                                                                             gAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagtt                                                                                                            <  2:208955/60‑1 (MQ=255)
                                                                                                                             gAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagtt                                                                                                            >  1:208955/1‑60 (MQ=255)
                                                                                                                             gAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGtt                              >  3:252234/1‑138 (MQ=255)
                                                                                                                             gAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGtt                              <  4:252234/138‑1 (MQ=255)
                                                                                                                             gAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCgcgc                   >  1:1488854/1‑149 (MQ=255)
                                                                                                                                 ttttGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAAc                        <  4:1434884/140‑1 (MQ=255)
                                                                                                                                   ttGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTAtt             <  3:1693383/149‑1 (MQ=255)
                                                                                                                                       ttAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttgagt                                                                                                        <  1:551219/54‑1 (MQ=255)
                                                                                                                                       ttAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTcc                                       <  4:1627216/119‑1 (MQ=255)
                                                                                                                                       ttAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTcc                                       >  3:1627216/1‑119 (MQ=255)
                                                                                                                                       ttAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTAtt             <  2:692144/145‑1 (MQ=255)
                                                                                                                                       ttAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTAtt             >  1:692144/1‑145 (MQ=255)
                                                                                                                                       ttAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGGAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCGATCGATCGTGTCGATGGCACTGTTAATCTGCCCGGTTCCAAGGGCGTTTCTAACCGCGCGTTATTGCtt         >  3:1376291/1‑148 (MQ=255)
                                                                                                                                        tAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTgg        <  2:1488854/149‑1 (MQ=255)
                                                                                                                                        tAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTgg        <  4:1232629/149‑1 (MQ=255)
                                                                                                                                        aaaaccccACAGGCTGGCTGTGGGGATTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGc                                 >  3:976064/5‑124 (MQ=255)
                                                                                                                                         aaTCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTaa                         >  2:1708557/1‑131 (MQ=255)
                                                                                                                                         aaTCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTaa                         <  1:1708557/131‑1 (MQ=255)
                                                                                                                                           accccACAGGCTGGCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATg                                                              <  4:1616350/91‑1 (MQ=255)
                                                                                                                                           accccACAGGCTGGCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATg                                                              <  4:1264962/91‑1 (MQ=255)
                                                                                                                                             cccACAGGCTGGCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTcc                                       <  4:702863/113‑1 (MQ=255)
                                                                                                                                             cccACAGGCTGGCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTcc                                       >  3:702863/1‑113 (MQ=255)
                                                                                                                                              ccACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACACCCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTc                                        <  2:536296/111‑1 (MQ=255)
                                                                                                                                                    ccagccTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTAt              <  4:1884376/131‑1 (MQ=255)
                                                                                                                                                    ccagccTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTAt              >  3:1884376/1‑131 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TGCACTGAAAGGTCACCGTGTGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGCTGGCGGCAT  >  NC_000913/958621‑958911

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