Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,156,366 G→A 5.5% G338S (GGC→AGC)  yqhD → aldehyde reductase, NADPH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,156,3660GA 5.5% 234.7 / 2.4 73G338S (GGC→AGC) yqhDaldehyde reductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base G (33/36);  new base A (2/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

CTGCCTGCACTGTGGAATGAAAAACGCGATACCAAGCGCGCTAAGCTGCTGCAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCC  >  NC_000913/3156219‑3156509
                                                                                                                                                   |                                                                                                                                               
ctgcctgcACTGTGGAATGAAAAACGCGATACCAAGCGCGCTAAGCTGCTGCAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAgg                                                                                                                                                >  2:193552/1‑149 (MQ=255)
   cctgcACTGTGGAATGAAAAACGCGATACCAAGCGCGCTAAGCTGCTGCAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGt                                                                                                                                             <  2:26668/149‑1 (MQ=255)
      gcACTGTGGAATGAAAAACGCGATACCAAGCGCGCTAAGCTGCTGCAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGcc                                                                                                                                          >  7:284881/1‑149 (MQ=255)
        aCTGTGGAATGAAAAACGCGATACCAAGCGCGCTAAGCTGCTGCAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGa                                                                                                                                        <  8:129681/149‑1 (MQ=255)
          tgtgGAATGAAAAACGCGATACCAAGCGCGCTAAGCTGCTGCAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGAcc                                                                                                                                      >  1:102856/1‑149 (MQ=255)
            tgGAATGAAAAACGCGATACCAAGCGCGCTAAGCTGCTGCAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGAcc                                                                                                                                      <  7:10596/147‑1 (MQ=255)
              gAATGAAAAACGCGATACCAAGCGCGCTAAGCTGCTGCAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCAc                                                                                                                                   <  7:135278/148‑1 (MQ=255)
               aaTGAAAAACGCGATACCAAGCGCGCTAAGCTGCTGCAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACct                                                                                                                                 >  6:318112/1‑149 (MQ=255)
                 tGAAAAACGCGATACCAAGCGCGCTAAGCTGCTGCAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACctct                                                                                                                               <  8:302557/149‑1 (MQ=255)
                               cctaGCGCGCTAAGCTGCTGCAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGGAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTg                                                                                                                 <  5:245211/146‑1 (MQ=255)
                               ccAAGCGCGCTAAGCTGCTGCAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTg                                                                                                                 <  1:35678/149‑1 (MQ=255)
                                       gcTAAGCTGCTGCAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTg                                                                                                                 >  7:373835/1‑141 (MQ=255)
                                       gcTAAGCTGCTGCAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTg                                                                                                                 <  8:373835/141‑1 (MQ=255)
                                              tgctgcAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGAAGATCCATc                                                                                                  >  5:318437/1‑149 (MQ=255)
                                                   cAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGc                                                                                             >  2:243996/1‑149 (MQ=255)
                                                      tatGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACGACGGTCTGGACGGCAGCTCCATCCCGGCttt                                                                                          <  2:242194/149‑1 (MQ=255)
                                                        tGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTg                                                                                         <  1:72731/148‑1 (MQ=255)
                                                        tGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCAACCCGGCTTTGc                                                                                        >  5:340334/1‑149 (MQ=255)
                                                        tGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGAGGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGc                                                                                        <  2:87290/149‑1 (MQ=255)
                                                         gCTGAACGCGTCTGGAACATCACTGAAGGTTACGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGAGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCt                                                                                       >  5:301981/1‑149 (MQ=255)
                                                          cTGAAGGCGTCTTGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCTCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCt                                                                                       <  6:301981/148‑1 (MQ=255)
                                                           tGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGCCCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTg                                                                                      >  3:217504/1‑148 (MQ=255)
                                                           tGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTg                                                                                      >  7:29496/1‑148 (MQ=255)
                                                              aCGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGaaaa                                                                                  <  5:318112/149‑1 (MQ=255)
                                                                gcgTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGaaaaaa                                                                                <  1:398314/149‑1 (MQ=255)
                                                                 cgTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAAc                                                                               <  1:243996/149‑1 (MQ=255)
                                                                  gTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACt                                                                              >  8:54939/1‑149 (MQ=255)
                                                                   tCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTg                                                                             <  7:64285/149‑1 (MQ=255)
                                                                    cTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTg                                                                             >  6:18687/1‑148 (MQ=255)
                                                                     tGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTgg                                                                            >  3:66156/1‑148 (MQ=255)
                                                                      ggAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTgg                                                                            >  3:348930/1‑147 (MQ=255)
                                                                       gAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTgg                                                                            >  5:57552/1‑146 (MQ=255)
                                                                          cATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAAc                                                                               >  4:326463/1‑140 (MQ=255)
                                                                          cATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAAc                                                                               <  3:326463/140‑1 (MQ=255)
                                                                          cATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGc                                                                      <  5:78811/149‑1 (MQ=255)
                                                                           aTCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCa                                                                     <  6:248162/149‑1 (MQ=255)
                                                                             cACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAAGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACg                                                                   >  2:80095/1‑149 (MQ=255)
                                                                             cACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAAGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACg                                                                   >  2:74858/1‑149 (MQ=255)
                                                                              aCTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCa                                                                                                                                    >  2:68932/1‑83 (MQ=255)
                                                                              aCTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCGTTGCTGAAAAAACTGGAAGAGCACgg                                                                  <  4:237612/149‑1 (MQ=255)
                                                                               cTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAAGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGc                                                                 <  1:74858/149‑1 (MQ=255)
                                                                               cTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAAGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGc                                                                 <  1:80095/149‑1 (MQ=255)
                                                                                tGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCa                                                                >  6:315529/1‑149 (MQ=255)
                                                                                  aaGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCAt                                                               <  2:102856/148‑1 (MQ=255)
                                                                                  aaGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCAGCAACCCGCAATTTCTTTGATCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCAt                                                               <  6:100655/148‑1 (MQ=255)
                                                                                   aGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGa                                                             <  5:315529/149‑1 (MQ=255)
                                                                                     gTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGAcc                                                           <  6:57552/149‑1 (MQ=255)
                                                                                     gTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGAcc                                                           >  5:369774/1‑149 (MQ=255)
                                                                                         agatgatGAGCGTATTGACGCCGATATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGAAACACCTCTCCGACTACGGAATGGAAGGAAGCTACATACCCTCTTTCCTGAAAAAACTGGAAGAGCACGGCATCCCCCAAc                                                       >  3:153401/2‑149 (MQ=255)
                                                                                               tgaGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCCGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGc                                                  >  5:211569/1‑148 (MQ=255)
                                                                                                gaGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCg                                                 >  7:69390/1‑148 (MQ=255)
                                                                                                  gCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACg                                                                   <  2:5655/128‑1 (MQ=255)
                                                                                                  gCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACg                                                                   >  1:5655/1‑128 (MQ=255)
                                                                                                  gCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCACGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGaa                                               >  4:147777/1‑148 (MQ=255)
                                                                                                     ttgttaCGCCGGTATTGCCGCATCCCGCAATTTGTTTGAGTTATTAGTCTTGCCGACGCACCTCTCCTACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGGTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAACATc                                           <  6:117863/144‑1 (MQ=255)
                                                                                                       ttGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCAt                                         <  5:258619/149‑1 (MQ=255)
                                                                                                         gACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGa                                       >  7:30442/1‑149 (MQ=255)
                                                                                                          aCGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGa                                       >  1:175047/1‑148 (MQ=255)
                                                                                                             ccgcGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACAtt                                   <  3:290753/149‑1 (MQ=255)
                                                                                                               gcgATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTAc                                 <  4:366201/149‑1 (MQ=255)
                                                                                                                   ttGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTg                             <  4:150653/149‑1 (MQ=255)
                                                                                                                      ccGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGa                           <  4:348930/148‑1 (MQ=255)
                                                                                                                      ccGCAACCCGCAATTTCTTTGAGAAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTATGTTGGa                           <  1:289744/148‑1 (MQ=255)
                                                                                                                               gCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAgccgc                 >  3:292016/1‑149 (MQ=255)
                                                                                                                                 aaTTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGccgccg               <  1:193552/149‑1 (MQ=255)
                                                                                                                                    ttCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGta             >  8:37988/1‑148 (MQ=255)
                                                                                                                                      cTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGtata           >  2:36857/1‑148 (MQ=255)
                                                                                                                                       tttGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCCGCCGCCGtatata         >  4:65811/1‑149 (MQ=255)
                                                                                                                                       tttGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGtatat          >  6:85519/1‑148 (MQ=255)
                                                                                                                                        ttGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGtatata         >  6:364015/1‑148 (MQ=255)
                                                                                                                                         tGAGCATTTAGGCGTGCCGACCCACCTCTCCGCCTACGTTCTGGACGGCAGCTCCATCCCGGCTTTTCTGAAAAAACTGTAAGAGCACGTCATGACCCAACTGTGCGAAAATCATGACATTACGTTTGATGTCAGCAGCCGtatata         <  6:318437/147‑1 (MQ=255)
                                                                                                                                            gCAATTAGGCGTGCCGACCCCCCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATAc        >  2:91254/1‑145 (MQ=255)
                                                                                                                                            gCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATAc        <  1:91254/145‑1 (MQ=255)
                                                                                                                                              aaTTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAgcc  <  5:18687/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CTGCCTGCACTGTGGAATGAAAAACGCGATACCAAGCGCGCTAAGCTGCTGCAATATGCTGAACGCGTCTGGAACATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCC  >  NC_000913/3156219‑3156509

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: