Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,655,299 G→A 5.9% intergenic (‑158/+48) rspA ← / ← ynfA bifunctional D‑altronate/D‑mannonate dehydratase/UPF0060 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,655,2990GA 5.9% 215.1 / 2.8 68intergenic (‑158/+48)rspA/ynfAbifunctional D‑altronate/D‑mannonate dehydratase/UPF0060 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (28/36);  new base A (2/2);  total (30/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.45e-01

GCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGC  >  NC_000913/1655154‑1655447
                                                                                                                                                 |                                                                                                                                                    
gCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGt                                                                                                                                                     <  5:138722/147‑1 (MQ=255)
   tCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCt                                                                                                                                                <  3:96705/149‑1 (MQ=255)
   tCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCt                                                                                                                                                <  2:48058/149‑1 (MQ=255)
     gCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGa                                                                                                                                              <  5:44718/149‑1 (MQ=255)
       ttGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTATCAGCTGaca                                                                                                                                            >  6:368392/1‑149 (MQ=255)
         gTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGaca                                                                                                                                            <  5:9256/147‑1 (MQ=255)
            gTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTc                                                                                                                                        <  5:132464/148‑1 (MQ=255)
             tGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAACGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTATCAGCTGACACTTCct                                                                                                                                      >  6:79510/1‑149 (MQ=255)
              gATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCctc                                                                                                                                     >  1:245597/1‑149 (MQ=255)
               aTGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCctc                                                                                                                                     <  3:180883/148‑1 (MQ=255)
                 gCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCtctta                                                                                                                                  >  5:148383/1‑149 (MQ=255)
                  cATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCtcttat                                                                                                                                 >  3:67110/1‑149 (MQ=255)
                   aTGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCtcttatc                                                                                                                                >  2:186258/1‑149 (MQ=255)
                   aTGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCtcttatc                                                                                                                                >  1:88489/1‑149 (MQ=255)
                       aaTCTACGCAACTGCGCTACTACCATACAAGTATAAAGATCGAAACAAGCCGGAGGGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGCTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCtcttatctta                                                                                                                             <  6:5718/148‑1 (MQ=255)
                        aTCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTAtt                                                                                                                           >  2:191959/1‑149 (MQ=255)
                        aTCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTAtt                                                                                                                           <  4:67110/149‑1 (MQ=255)
                          cTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGa                                                                                                                         >  2:347908/1‑149 (MQ=255)
                           tACGCAACTGAGCTACTACCATACAAGTATTAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGa                                                                                                                         >  4:196039/1‑148 (MQ=255)
                              gCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTcttatcttat                                                                                                                            >  6:251485/1‑142 (MQ=255)
                              gCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTcttatcttat                                                                                                                            <  5:251485/142‑1 (MQ=255)
                               cAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATaaaa                                                                                                                    <  6:148383/149‑1 (MQ=255)
                               cAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATaaaa                                                                                                                    <  8:12005/149‑1 (MQ=255)
                                  cTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATgg                                                                                                                 <  7:238716/149‑1 (MQ=255)
                                   tGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGa                                                                                                                >  1:337361/1‑149 (MQ=255)
                                     aGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGTAGTGATCACAAAAAAAGGCGTATATTTGCGCTTTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGAt                                                                                                               <  3:6134/148‑1 (MQ=255)
                                          ctacCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTATCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTAt                                                                                                           <  5:368392/147‑1 (MQ=255)
                                          atacCATAAAATTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTAt                                                                                                           <  6:66825/146‑1 (MQ=255)
                                           tacCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGtt                                                                                                        >  6:332611/1‑149 (MQ=255)
                                             ccATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTATCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGtt                                                                                                        <  5:79510/147‑1 (MQ=255)
                                              cATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCt                                                                                                     <  7:189913/149‑1 (MQ=255)
                                              cATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGAAACTTCCTCTTATCTTATTGATAAAATGGATTTagatagg                                                                                                      >  3:350668/1‑142 (MQ=255)
                                              cATAAAAGTATAAAGATCGAAAAAAGCCGGATTGATCACAAAAAAAGGCGTATATTTGCGCTGTGATTGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGAtt                                                                                                              <  4:350668/140‑1 (MQ=255)
                                               aTACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTa                                                                                                    >  7:46998/1‑149 (MQ=255)
                                                tACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTAc                                                                                                   >  4:235669/1‑149 (MQ=255)
                                                  cAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTg                                                                                                                          <  8:98660/124‑1 (MQ=255)
                                                  cAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTg                                                                                                                          >  7:98660/1‑124 (MQ=255)
                                                  cAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGt                                                                                                 <  1:347908/149‑1 (MQ=255)
                                                   aaGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTg                                                                                                <  1:191959/149‑1 (MQ=255)
                                                        taAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGcc                                                                                            <  5:332611/148‑1 (MQ=255)
                                                          aaGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCctc                                                                                                                                     <  4:372583/105‑1 (MQ=255)
                                                          aaGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCctc                                                                                                                                     >  3:372583/1‑105 (MQ=255)
                                                                gAAAAAAGCCGGAGTGATGACAAAAAAAGGCGTATATTTGCGCTGTGATTGGTTGACAAATGATGAAATAGAATACCGTTTGTCTGCTGCCACTTCCTCTTATCGTATTGATAAGAGGGATTTATGTTCCTACGTGCGCCCcccgcccg                                                                                   >  6:379763/1‑149 (MQ=255)
                                                                 aaaaaaGCCGGAGTGATCACAAAAAAAGGGGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCATCTGACTCTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGc                                                                                  <  6:87479/149‑1 (MQ=255)
                                                                       gCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAAt                                                                            <  2:245597/149‑1 (MQ=255)
                                                                          ggAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGAt                                                                         <  2:353004/149‑1 (MQ=255)
                                                                            aGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATc                                                                        >  6:161322/1‑148 (MQ=255)
                                                                             gTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCaa                                                                      <  3:84211/149‑1 (MQ=255)
                                                                              tGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAAc                                                                     <  2:337361/149‑1 (MQ=255)
                                                                                aTCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACa                                                                    <  3:235669/148‑1 (MQ=255)
                                                                                 tCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATTCCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGCTTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATg                                                                  >  6:111455/1‑149 (MQ=255)
                                                                                   acaAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGc                                                                 <  1:186258/148‑1 (MQ=255)
                                                                                     aaaaaaaGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCg                                                               >  3:197919/1‑148 (MQ=255)
                                                                                             gCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAgcgc                                                       >  4:288106/1‑148 (MQ=255)
                                                                                                   atTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAAc                                                >  1:81155/1‑149 (MQ=255)
                                                                                                    tttGCGCTGTGAATGGTTGAGAAAAGATGAAATAGAATACCTTTTGTCAGCTGACCCTTCCTCTTATCTTATTGATAAAATGGATTTATGTACCTAAGTGCGCCCCCAGCCCGCAACAATGATCACCATGCCGCAAAGCGCAATCAACg                                               >  6:134132/1‑149 (MQ=255)
                                                                                                             tgAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTc                                       >  3:20365/1‑148 (MQ=255)
                                                                                                                      gACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTc                                  >  2:61887/1‑144 (MQ=255)
                                                                                                                      gACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTc                                  <  1:61887/144‑1 (MQ=255)
                                                                                                                        cAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAaga                           <  2:88489/149‑1 (MQ=255)
                                                                                                                           aaGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGt                         >  4:382054/1‑148 (MQ=255)
                                                                                                                              aTGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGt                     <  5:161322/149‑1 (MQ=255)
                                                                                                                                 aaaTAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTc                  <  2:81155/149‑1 (MQ=255)
                                                                                                                                   aTAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCAc                <  3:382054/149‑1 (MQ=255)
                                                                                                                                            cTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGt                     <  1:117201/135‑1 (MQ=255)
                                                                                                                                            cTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGt                     >  2:117201/1‑135 (MQ=255)
                                                                                                                                            cTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCcac       <  8:46998/149‑1 (MQ=255)
                                                                                                                                            cTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCca        >  4:87015/1‑148 (MQ=255)
                                                                                                                                                 gTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTACAGTCATAAAGAGTCAGTTTCACGCCATCCACAAcgc  >  8:286314/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGC  >  NC_000913/1655154‑1655447

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: ‑