Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 146,733 G→A 5.1% 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 5.1% 394.3 / 6.8 119intergenic (‑39/‑235)panD/yadDaspartate 1‑decarboxylase/transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (59/53);  new base A (4/2);  total (63/56)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 119 total aligned reads displayed.

CGGCTGCGTCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTGAAC  >  NC_000913/146586‑146877
                                                                                                                                                   |                                                                                                                                                 
cGGCTGCGTCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGa                                                                                                                                                  <  3:203097/149‑1 (MQ=255)
   cTGCGTCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTCCCCTGTCGTTAACGAAGCAGGCGATGc                                                                                                                                              <  4:324623/149‑1 (MQ=255)
   cTGCGTCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGc                                                                                                                                              <  1:550188/149‑1 (MQ=255)
          aaGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCttt                                                                                                                                       >  2:379094/1‑149 (MQ=255)
          aaGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCttt                                                                                                                                       >  2:454931/1‑149 (MQ=255)
          aaGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCATCGCCTGCttc                                                                                                                                       <  2:262955/149‑2 (MQ=255)
            gAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGa                                                                                                                                     <  1:119759/149‑1 (MQ=255)
            gAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGa                                                                                                                                     >  2:804742/1‑149 (MQ=255)
            gAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGa                                                                                                                                     >  3:809322/1‑149 (MQ=255)
            gAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGa                                                                                                                                     <  3:322423/149‑1 (MQ=255)
                aTCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACAAAGCAGGCATCGCCTGCttcgttaac                                                                                                                                 <  1:431627/149‑8 (MQ=255)
                  ccTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAAt                                                                                                                               >  1:180386/1‑149 (MQ=255)
                      gTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATtctc                                                                                                                           >  2:609065/1‑149 (MQ=255)
                      gTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATtctc                                                                                                                           <  1:664219/149‑1 (MQ=255)
                       tCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCa                                                                                                                          <  2:25236/149‑1 (MQ=255)
                       tCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCa                                                                                                                          >  3:735317/1‑149 (MQ=255)
                       tCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCa                                                                                                                          <  3:806139/149‑1 (MQ=255)
                       tCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCa                                                                                                                          <  4:679216/149‑1 (MQ=255)
                          aTGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGc                                                                                                                       >  1:191987/1‑149 (MQ=255)
                           tGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCa                                                                                                                      <  4:332546/149‑1 (MQ=255)
                              cgcAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTa                                                                                                                   >  2:618497/1‑149 (MQ=255)
                                 aaGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTg                                                                                                                >  3:477226/1‑149 (MQ=255)
                                  aGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGc                                                                                                               <  3:469365/149‑1 (MQ=255)
                                      ccTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGa                                                                                                           >  1:1435/1‑149 (MQ=255)
                                       cTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGAt                                                                                                          >  1:816842/1‑149 (MQ=255)
                                       cTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGAt                                                                                                          >  2:259527/1‑149 (MQ=255)
                                        ttCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGAtt                                                                                                         >  4:711179/1‑149 (MQ=255)
                                             aGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGt                                                                                                    <  4:477226/149‑1 (MQ=255)
                                              gTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGtt                                                                                                   <  1:91211/149‑1 (MQ=255)
                                                 cAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAg                                                                                                <  2:816842/149‑1 (MQ=255)
                                                 cAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAg                                                                                                <  1:721852/149‑1 (MQ=255)
                                                  aGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTAtt                                                                                                                 >  4:425272/1‑131 (MQ=255)
                                                  aGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTAtt                                                                                                                 <  3:425272/131‑1 (MQ=255)
                                                  aGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAgt                                                                                               >  1:660898/1‑149 (MQ=255)
                                                  aGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAgt                                                                                               <  1:779245/149‑1 (MQ=255)
                                                    gTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGAttttt                                                                                                      <  3:80785/140‑1 (MQ=255)
                                                    gTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGAttttt                                                                                                      >  4:80785/1‑140 (MQ=255)
                                                    gTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAgtgt                                                                                             <  4:380891/149‑1 (MQ=255)
                                                      ccGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCt                                                                                           <  1:692298/149‑1 (MQ=255)
                                                          aTGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTc                                                                                       >  3:141265/1‑149 (MQ=255)
                                                           tGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCa                                                                                      >  2:440341/1‑149 (MQ=255)
                                                            gAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGTGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCAt                                                                                     >  2:579433/1‑149 (MQ=255)
                                                            gAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCAt                                                                                     >  2:695542/1‑149 (MQ=255)
                                                              gTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCt                                                                                   >  1:25380/1‑149 (MQ=255)
                                                               tCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTg                                                                                  >  2:646264/1‑149 (MQ=255)
                                                               tCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTg                                                                                  >  2:667258/1‑149 (MQ=255)
                                                                 aCTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGAc                                                                                >  1:443530/1‑149 (MQ=255)
                                                                 aCTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGAc                                                                                >  3:173182/1‑149 (MQ=255)
                                                                 aCTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGAc                                                                                <  1:646264/149‑1 (MQ=255)
                                                                   tttCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACgg                                                                              >  2:9607/1‑149 (MQ=255)
                                                                    ttCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGc                                                                             <  1:609065/149‑1 (MQ=255)
                                                                       aCGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCAtt                                                                          <  2:452976/149‑1 (MQ=255)
                                                                          cgGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGc                                                                       <  2:191987/149‑1 (MQ=255)
                                                                             tGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTc                                                                    <  1:311199/149‑1 (MQ=255)
                                                                                 gTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTcagca                                                                >  3:695448/1‑149 (MQ=255)
                                                                                  tttGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTc                                                                                            <  2:817587/120‑1 (MQ=255)
                                                                                  tttGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTc                                                                                            >  1:817587/1‑120 (MQ=255)
                                                                                    tGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGtt                                                             <  1:259527/149‑1 (MQ=255)
                                                                                    tGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGtt                                                             <  2:480091/149‑1 (MQ=255)
                                                                                       ccTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGc                                                          >  1:250343/1‑149 (MQ=255)
                                                                                       ccTGCACCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGc                                                          >  1:412346/1‑149 (MQ=255)
                                                                                        cTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCg                                                         >  1:44249/1‑149 (MQ=255)
                                                                                          gcagcaTCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTa                                                       <  4:237434/149‑1 (MQ=255)
                                                                                               aTCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACcgcg                                                  >  4:463749/1‑149 (MQ=255)
                                                                                                 cGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGcc                                                <  2:443530/149‑1 (MQ=255)
                                                                                                        aTCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGt                                         <  2:841353/149‑1 (MQ=255)
                                                                                                        aTCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGt                                         <  2:660898/149‑1 (MQ=255)
                                                                                                         tCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTg                                        <  4:173182/149‑1 (MQ=255)
                                                                                                           aTAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCCGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGcc                                      >  1:604118/1‑149 (MQ=255)
                                                                                                           aTAACTTCTACCTTTATACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGcc                                      <  4:561602/149‑1 (MQ=255)
                                                                                                            tAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCa                                     >  4:213286/1‑149 (MQ=255)
                                                                                                               cTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTGTTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGt                                  >  1:572245/1‑149 (MQ=255)
                                                                                                               cTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGt                                  >  3:375935/1‑149 (MQ=255)
                                                                                                                   tACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTGACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCg                              >  2:802145/1‑149 (MQ=255)
                                                                                                                       tttCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTg                          >  1:392496/1‑149 (MQ=255)
                                                                                                                         tCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTa                        >  2:642425/1‑149 (MQ=255)
                                                                                                                         tCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTa                        <  3:668465/149‑1 (MQ=255)
                                                                                                                          cTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTaa                       <  1:804742/149‑1 (MQ=255)
                                                                                                                           tACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAAc                      >  2:153858/1‑149 (MQ=255)
                                                                                                                            aCCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGt                                         >  1:223589/1‑129 (MQ=255)
                                                                                                                             cccTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGa                    >  3:793878/1‑149 (MQ=255)
                                                                                                                              ccTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGaa                   >  1:301955/1‑149 (MQ=255)
                                                                                                                               cTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAAt                  >  2:489720/1‑149 (MQ=255)
                                                                                                                               cTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAAt                  >  2:429039/1‑149 (MQ=255)
                                                                                                                                tGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAAtt                 >  3:407854/1‑149 (MQ=255)
                                                                                                                                tGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAAtt                 <  2:563089/149‑1 (MQ=255)
                                                                                                                                    gTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGata             >  4:183172/1‑149 (MQ=255)
                                                                                                                                     ttAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGatat            >  1:652405/1‑149 (MQ=255)
                                                                                                                                     ttAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGatat            >  4:588638/1‑149 (MQ=255)
                                                                                                                                       aaCGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGTGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATAttt          <  4:674907/149‑1 (MQ=255)
                                                                                                                                       aaCGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATAttt          >  2:755625/1‑149 (MQ=255)
                                                                                                                                        tcaaaGCAGGCATCGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTtc         >  3:899435/4‑149 (MQ=255)
                                                                                                                                        tcaaaGCAGGCATCGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTtc         >  4:762229/4‑149 (MQ=255)
                                                                                                                                         cGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTtct        <  1:489720/149‑1 (MQ=255)
                                                                                                                                          gAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTtctc       <  1:864691/149‑1 (MQ=255)
                                                                                                                                           aaGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTctct      <  1:153858/149‑1 (MQ=255)
                                                                                                                                           aaGCAGGCATCGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTctct      >  2:312072/1‑149 (MQ=255)
                                                                                                                                            aGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTg     >  1:540455/1‑149 (MQ=255)
                                                                                                                                              cAGGCGATACCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTTCTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGTCCAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTGaa   <  4:569072/149‑1 (MQ=255)
                                                                                                                                               aGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTGAAc  <  1:9607/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGGCTGCGTCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTGAAC  >  NC_000913/146586‑146877

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