Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,138,988 G→A 6.3% G363S (GGC→AGC)  yegH → inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,138,9880GA 6.3% 197.3 / 2.6 63G363S (GGC→AGC) yegHinner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (26/33);  new base A (2/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

TATCGATCTCAACGCGCCGGAAGAAGAAATCCGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGT  >  NC_000913/2138843‑2139136
                                                                                                                                                 |                                                                                                                                                    
tatCGATCTCAACGCGCCGGAAGAAGAAATCCGCCAGTTGCTGGAGCGAAATCAGCATACGCGTCTGGTGGTTACCGAAGGTGAAGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCgcggc                                                                                                                                                    >  1:73351/1‑148 (MQ=255)
tatCGATCTCAACGCGCCGGAAGAAGAAATCCGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCgcgg                                                                                                                                                     >  5:385384/1‑147 (MQ=255)
tatCGATCTCAACGCGCCGGAAGAAGAAATCCGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCgcgg                                                                                                                                                     <  6:385384/147‑1 (MQ=255)
 atCGATCTCAACGCGCCGGAAGAAGAAATCCGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCgcggcg                                                                                                                                                   >  1:19254/1‑148 (MQ=255)
     aTCTCAACGCGCCGGAAGAAGAAATCCGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCg                                                                                                                                              >  7:231503/1‑149 (MQ=255)
          aaCGCGCCGGAAGAAGAAATCCGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCaa                                                                                                                                         >  7:201537/1‑149 (MQ=255)
           aCGCGCCGGAAGAAGAAATCCGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGTGTTTACCGACGGTGATGACGCAGACGATTTGCTCGGTGTTGTTCACGTAATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAAc                                                                                                                                        <  4:239799/149‑1 (MQ=255)
                 cGGAAGAAGAAATCCGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCATCCTGCgg                                                                                                                                  <  2:311354/149‑1 (MQ=255)
                 cGGAAGAAGAAATCCGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCgg                                                                                                                                  <  4:315199/149‑1 (MQ=255)
                   gaagaagaaATCCGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGgt                                                                                                                                <  3:197951/149‑1 (MQ=255)
                   gaagaagaaATCCGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGgt                                                                                                                                <  6:311272/149‑1 (MQ=255)
                     agaagaaATCCGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGgtgt                                                                                                                              <  2:28016/149‑1 (MQ=255)
                      gaagaaATCCGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGtt                                                                                                                             >  7:190362/1‑149 (MQ=255)
                          aaaTCCGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGTGGGTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGAtt                                                                                                                         <  3:250922/149‑1 (MQ=255)
                               cGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGCCCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGc                                                                                                                      >  3:75662/1‑147 (MQ=255)
                               cGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGGTGTTACCGACGGTGATGACGCATTAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCa                                                                                                                    <  2:73351/149‑1 (MQ=255)
                                     ttGCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCt                                                                                                              >  7:17925/1‑149 (MQ=255)
                                       gCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTg                                                                                                             <  8:231503/148‑1 (MQ=255)
                                       gCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTg                                                                                                             <  4:75662/148‑1 (MQ=255)
                                            aGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTc                                                                                                       >  7:343343/1‑149 (MQ=255)
                                            aGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTc                                                                                                       <  7:385904/149‑1 (MQ=255)
                                             gCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTcc                                                                                                      >  5:69246/1‑149 (MQ=255)
                                              cGAAATCAGCATACGCGGCTGGTTGATACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTccc                                                                                                     >  4:228788/1‑149 (MQ=255)
                                                  aTCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGa                                                                                                  >  7:379067/1‑148 (MQ=255)
                                                      gCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGgtg                                                                                                          <  7:246233/136‑1 (MQ=255)
                                                      gCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGgtg                                                                                                          >  8:246233/1‑136 (MQ=255)
                                                      gCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCt                                                                                             >  8:31296/1‑149 (MQ=255)
                                                           cgcgGCTGGGGGTTACCGACGTTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGTCTAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGcc                                                                                         <  6:77843/148‑1 (MQ=255)
                                                            gcgGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGGGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCgt                                                                                       <  2:19254/149‑1 (MQ=255)
                                                               gCTGGTGGTTACCGACGGTGATGACGCAGAAGCTTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCgttgt                                                                                    <  6:158048/149‑1 (MQ=255)
                                                               gCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCgttgt                                                                                    <  8:49702/149‑1 (MQ=255)
                                                                  ggtggtTACCGACGGTGATGACGCAGAATATTTGCTCGGTGATGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAAGCTCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTAc                                                                                 <  2:249343/149‑1 (MQ=255)
                                                                    tggtTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCt                                                                               >  1:186558/1‑149 (MQ=255)
                                                                       ttACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGCCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGcc                                                                            >  4:333803/1‑149 (MQ=255)
                                                                          ccGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCCCGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTg                                                                         >  1:214547/1‑149 (MQ=255)
                                                                            gACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTg                                                                         >  8:76511/1‑147 (MQ=255)
                                                                             aCGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAg                                                                      <  7:76511/149‑1 (MQ=255)
                                                                                    tgaCGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTcc                                                               >  4:256177/1‑149 (MQ=255)
                                                                                        gCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTaa                                                           >  5:102299/1‑149 (MQ=255)
                                                                                           gaagaTTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCAGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGc                                                        >  3:266514/1‑149 (MQ=255)
                                                                                           gaagaTTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCAGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGc                                                        >  3:271592/1‑149 (MQ=255)
                                                                                            aagaTTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGcc                                                       >  2:315804/1‑149 (MQ=255)
                                                                                               aTTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCcgc                                                    <  4:288934/149‑1 (MQ=255)
                                                                                                       ggTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACtt                                            <  4:116063/149‑1 (MQ=255)
                                                                                                       ggTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACtt                                            <  1:403309/149‑1 (MQ=255)
                                                                                                        gtgtTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACttt                                           >  3:19689/1‑149 (MQ=255)
                                                                                                           ttgttCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGc                                         <  6:102299/148‑1 (MQ=255)
                                                                                                           ttgttCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGc                                         <  5:133481/148‑1 (MQ=255)
                                                                                                            tgttCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTg                                                                                                             >  5:70395/1‑79 (MQ=255)
                                                                                                            tgttCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTg                                                                                                             <  6:70395/79‑1 (MQ=255)
                                                                                                             gttCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCttt                                      >  7:218754/1‑149 (MQ=255)
                                                                                                             gttCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCtt                                       <  8:201537/148‑1 (MQ=255)
                                                                                                                 aCGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTtgt                                   <  3:256177/148‑1 (MQ=255)
                                                                                                                       tCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGa                             <  8:190362/148‑1 (MQ=255)
                                                                                                                         gACCTGCTACAACAGTCTCTGCGCAGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGATGa                          <  4:271592/149‑1 (MQ=255)
                                                                                                                         gACCTGCTACAACAGTCTCTGCGCAGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGATGa                          <  4:266514/149‑1 (MQ=255)
                                                                                                                         gACCTGCTACAACAGGCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGATGa                          <  1:184463/149‑1 (MQ=255)
                                                                                                                         gACCTGCTACAACAGCCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGATGa                          <  3:333803/149‑1 (MQ=255)
                                                                                                                               cTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGATGAGTTTgg                    >  7:271044/1‑149 (MQ=255)
                                                                                                                                tACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGCGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGATGAGTTTGGc                   >  1:47132/1‑149 (MQ=255)
                                                                                                                                   aacaGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGATGAGTTTGGcac                 >  2:64206/1‑146 (MQ=255)
                                                                                                                                   aacaGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGATGAGTTTGGcac                 <  1:64206/148‑3 (MQ=255)
                                                                                                                                      aGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGATGAGTTTGGCTCGGTg             <  1:380313/149‑1 (MQ=255)
                                                                                                                                            tGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGATGAGTTTGGCTCGGTGGaa          <  8:17925/146‑1 (MQ=255)
                                                                                                                                            tGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGATGAGTTTGGCTCGGTGGaa          <  8:218754/146‑1 (MQ=255)
                                                                                                                                            tGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGATGAGTTTGGCTCGGTGGAAgg        >  7:81239/1‑148 (MQ=255)
                                                                                                                                             gcgcGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTACCCGGAAACCTTGCCGTTGTTACCTGCCCAGGAGCAGTTCCGTAAAGCCCGCACGCACTATCCTTTTGTGGTGGATGAGTTTGGCTCCGTGGAAGGGa      >  6:307708/1‑149 (MQ=255)
                                                                                                                                               gcggcgAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGAtt    <  8:343343/149‑1 (MQ=255)
                                                                                                                                               gcggcgAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGAtt    <  1:212726/149‑1 (MQ=255)
                                                                                                                                                cggcgAACCGCTCAAGCTGCGGGTGTTGATTCGCCAGCCGGTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCGCGCACGCACTTTGCTTTTGTGGTGGTTGAGTTTGGCTCGGTGGAAGGGATtg   <  5:307708/149‑1 (MQ=255)
                                                                                                                                                 ggcgAACCGCGCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGGTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGGGGAGGAGTTTGGCTCGGGGGAAGGGATtgt  <  5:350173/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TATCGATCTCAACGCGCCGGAAGAAGAAATCCGCCAGTTGCTGGAGCGAAATCAGCATACGCGGCTGGTGGTTACCGACGGTGATGACGCAGAAGATTTGCTCGGTGTTGTTCACGTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTTTTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGT  >  NC_000913/2138843‑2139136

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