Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 146,733 G→A 6.2% intergenic (‑39/‑235) panD ← / → yadD aspartate 1‑decarboxylase/transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913146,7330GA 6.2% 197.5 / 2.4 66intergenic (‑39/‑235)panD/yadDaspartate 1‑decarboxylase/transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (33/27);  new base A (2/2);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

CGGCTGCGTCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTGAACGC  >  NC_000913/146586‑146879
                                                                                                                                                   |                                                                                                                                                  
cGGCTGCGTCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGa                                                                                                                                                   >  5:153120/1‑149 (MQ=255)
   cTGCGTCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCACGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGc                                                                                                                                                <  6:125108/149‑1 (MQ=255)
       gTCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGc                                                                                                                                            >  7:28806/1‑149 (MQ=255)
        tCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCt                                                                                                                                           >  1:246889/1‑149 (MQ=255)
          aaGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCATCGCCTGCtt                                                                                                                                          <  8:47458/148‑1 (MQ=255)
            gAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCttt                                                                                                                                         >  7:149971/1‑147 (MQ=255)
            gAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCttt                                                                                                                                         >  7:209219/1‑147 (MQ=255)
            gAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCATCGCCTGCttcgt                                                                                                                                       <  7:39730/149‑4 (MQ=255)
              aaaTCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAgg                                                                                                                                     <  4:45219/149‑1 (MQ=255)
                 tCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGaaa                                                                                                                                  <  5:188169/149‑1 (MQ=255)
                         aaTGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCAc                                                                                                                           <  7:149711/148‑1 (MQ=255)
                          aTGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACg                                                                                                                          >  6:151356/1‑148 (MQ=255)
                               gcAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTAt                                                                                                                    >  7:58316/1‑149 (MQ=255)
                                      ccTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGa                                                                                                             >  3:90292/1‑149 (MQ=255)
                                       cTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACACTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGAt                                                                                                            >  6:165450/1‑149 (MQ=255)
                                       cTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCAGAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGAt                                                                                                            >  6:165451/1‑149 (MQ=255)
                                              gTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGtt                                                                                                     <  5:208598/149‑1 (MQ=255)
                                              gTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCACTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCAGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGACTATTGtt                                                                                                     >  6:193346/1‑149 (MQ=255)
                                                 cAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTa                                                                                                   <  8:230798/148‑1 (MQ=255)
                                                       cGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTa                                                                                            <  5:151356/149‑1 (MQ=255)
                                                       cGCATGAGCCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCTATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTa                                                                                            >  6:84345/1‑149 (MQ=255)
                                                                 aCTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGAc                                                                                  >  2:96494/1‑149 (MQ=255)
                                                                   tttCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACgg                                                                                >  3:67334/1‑149 (MQ=255)
                                                                       aCGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCAtt                                                                            <  1:96494/149‑1 (MQ=255)
                                                                        cgcgGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCAttt                                                                           >  4:66354/1‑149 (MQ=255)
                                                                          cgGTGGAGTTTTCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGc                                                                         <  2:96867/149‑1 (MQ=255)
                                                                            gTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCTATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGt                                                                       <  5:84345/149‑1 (MQ=255)
                                                                            gTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGt                                                                       >  4:179975/1‑149 (MQ=255)
                                                                             tGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTc                                                                      <  3:66354/149‑1 (MQ=255)
                                                                              ggAGTTTGCCCTGCAGCATCGGGCGAATCATAACTTCTACCTTTCTACCCTGGCGTTAACGAAGCAGGCGAGGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTca                                                                     <  2:185236/149‑1 (MQ=255)
                                                                                  tttGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTcagcag                                                                 >  2:17002/1‑149 (MQ=255)
                                                                                       ccTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGc                                                            >  4:229793/1‑149 (MQ=255)
                                                                                         tGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGttt                                                              <  2:2449/145‑1 (MQ=255)
                                                                                         tGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGttt                                                              >  1:2449/1‑145 (MQ=255)
                                                                                         tGCAGCATCGTGCGAATCATAACTTCTACCTTTCTAAACTGTCGTTAACGAAGCAGGAGCTCCCTGCTTTGAGGAAAAACTCCCGCAGTATTGCCCGATTTTTGTATAGTGTCTACTCATCTGAAGGCATTTGCGTCAGAAGTTTGCg                                                           >  3:194299/1‑148 (MQ=255)
                                                                                          gcagcaTCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCAttt                                                                           <  7:160125/131‑1 (MQ=255)
                                                                                          gcagcaTCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACCAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCAttt                                                                           >  8:160125/1‑131 (MQ=255)
                                                                                              caTCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGTGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACcgc                                                     >  6:184854/1‑149 (MQ=255)
                                                                                                   tGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCg                                                 >  7:123512/1‑148 (MQ=255)
                                                                                                           atatCTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTTTTGCCCGTTTTTTGTTTAGTGTCTACTCACCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGc                                         <  5:57534/148‑1 (MQ=255)
                                                                                                           aTAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGc                                         <  1:76565/148‑1 (MQ=255)
                                                                                                              aCTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCAt                                      <  3:179975/148‑1 (MQ=255)
                                                                                                               cTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCg                                                           <  7:52067/126‑1 (MQ=255)
                                                                                                               cTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCg                                                           >  8:52067/1‑126 (MQ=255)
                                                                                                                ttCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGtt                                   >  2:223447/1‑148 (MQ=255)
                                                                                                                  cTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCa                                       >  8:39069/1‑143 (MQ=255)
                                                                                                                  cTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCa                                       <  7:39069/143‑1 (MQ=255)
                                                                                                                   tACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCg                                >  2:86848/1‑149 (MQ=255)
                                                                                                                   tACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCg                                >  6:71032/1‑149 (MQ=255)
                                                                                                                      cTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTTCCTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCt                             >  7:109417/1‑149 (MQ=255)
                                                                                                                             cccTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTaa                         <  8:28806/146‑1 (MQ=255)
                                                                                                                             cccTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTaa                         <  3:229793/146‑1 (MQ=255)
                                                                                                                             cccTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTaa                         <  8:209219/146‑1 (MQ=255)
                                                                                                                             cccTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTaa                         <  8:149971/146‑1 (MQ=255)
                                                                                                                             cccTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACg                       >  8:20693/1‑148 (MQ=255)
                                                                                                                               cTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAAt                    <  5:165450/149‑1 (MQ=255)
                                                                                                                               cTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAAt                    <  5:165451/149‑1 (MQ=255)
                                                                                                                                  tCGTTAACGAAGCATGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGa                 <  6:224328/149‑1 (MQ=255)
                                                                                                                                        tcaaaGCAGGCATCGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTtc           >  5:140802/4‑149 (MQ=255)
                                                                                                                                        tcaaaGCAGGCATCGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTtc           >  2:214067/4‑149 (MQ=255)
                                                                                                                                        aCGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTtc           >  8:100938/1‑149 (MQ=255)
                                                                                                                                          gAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTGATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTtct          <  3:97428/148‑1 (MQ=255)
                                                                                                                                          gAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCACTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTtct          >  4:97428/1‑148 (MQ=255)
                                                                                                                                              cAGGCGATGCCTGCTTTGCGGAAATTCTCACGCAGAATTGCCCGATATTTGTTAAGTGTCTACTCATCTGACGGCATATGCTTCAGCAGTTTGCGTACAGCGCCGAAGCGAGCCATGTATCGCCTGTAACGAATAGATATTTCTCTCaa     >  1:26732/1‑149 (MQ=255)
                                                                                                                                                 gCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTGAACgc  <  8:123512/149‑1 (MQ=255)
                                                                                                                                                 gCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTGAACgc  >  4:128933/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGGCTGCGTCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTGAACGC  >  NC_000913/146586‑146879

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: